teaching:progappchim:polynomes-8

Différences

Ci-dessous, les différences entre deux révisions de la page.

Lien vers cette vue comparative

Les deux révisions précédentes Révision précédente
Prochaine révision
Révision précédente
Dernière révisionLes deux révisions suivantes
teaching:progappchim:polynomes-8 [2012/11/16 07:09] villersdteaching:progappchim:polynomes-8 [2016/02/23 13:22] villersd
Ligne 13: Ligne 13:
     cf. http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Ruffini-Horner     cf. http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Ruffini-Horner
     """     """
-    n=len(a)-1 # n = ordre du polynome +    n = len(a)-1 # n = ordre du polynome 
-    p=a[n] +    p = 0. 
-    for i in range(n-1,-1,-1): +    for i in range(n,-1,-1): 
-        p=p*x+a[i]+        p = p*x + a[i]
     return p     return p
  
-def polyscal(a,s): +def polyscal(s,a): 
-    """polynome multiplié par un scalaire s """ +    """polynôme multiplié par un scalaire s """ 
-    b=[]+    b = []
     for coef in a:     for coef in a:
         b.append(coef*s)         b.append(coef*s)
Ligne 29: Ligne 29:
     """ Addition de deux polynomes de coefficients a et b     """ Addition de deux polynomes de coefficients a et b
     """     """
-    r=a[:] # on travaille sur une copie de a pour ne pas le modifier +    r = a[:] # on travaille sur une copie de a pour ne pas le modifier 
-    t=b[:] # idem pour b  +    t = b[:] # idem pour b  
-    g=[]   # polynome somme +    g = []   # polynome somme 
-    n1=len(r) # ordre du premier polynome +    n1 = len(r) # ordre du premier polynome 
-    n2=len(t) # ordre du second polynome +    n2 = len(t) # ordre du second polynome 
-    if n1>n2: # premier polynome de plus haut degré que le second+    if n1 > n2: # premier polynome de plus haut degré que le second
         for i in range (n1-n2):         for i in range (n1-n2):
             t.append(0)             t.append(0)
-    elif n1<n2: # second polynome de plus haut degré que le premier+    elif n1 < n2: # second polynome de plus haut degré que le premier
          for i in range (n2-n1):          for i in range (n2-n1):
              r.append(0)              r.append(0)
Ligne 47: Ligne 47:
 # différents tests : # différents tests :
  
-absc=[] +absc = [] 
-ordo=[] +ordo = [] 
-coef=[0,9,0,-120,0,432,0,-576,0,256]   # un polynome un peu particulier +coef = [0,9,0,-120,0,432,0,-576,0,256]   # un polynome un peu particulier 
-xa=-1.04+xa = -1.04
 while xa < 1.04: while xa < 1.04:
     absc.append(xa)     absc.append(xa)
     ordo.append(polyeval(xa,coef))     ordo.append(polyeval(xa,coef))
-    xa=xa+0.01+    xa = xa + 0.01
  
 plot(absc,ordo) plot(absc,ordo)
  • teaching/progappchim/polynomes-8.txt
  • Dernière modification : 2017/02/24 11:43
  • de villersd