Différences
Ci-dessous, les différences entre deux révisions de la page.
Les deux révisions précédentes Révision précédente Prochaine révision | Révision précédenteDernière révisionLes deux révisions suivantes | ||
teaching:progappchim:polynomes-8 [2012/11/16 07:09] – villersd | teaching:progappchim:polynomes-8 [2016/02/23 13:22] – villersd | ||
---|---|---|---|
Ligne 13: | Ligne 13: | ||
cf. http:// | cf. http:// | ||
""" | """ | ||
- | 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, |
- | p=p*x+a[i] | + | p = p*x + a[i] |
return p | return p | ||
- | def polyscal(a,s): | + | def polyscal(s,a): |
- | """ | + | """ |
- | b=[] | + | b = [] |
for coef in a: | for coef in a: | ||
b.append(coef*s) | b.append(coef*s) | ||
Ligne 29: | Ligne 29: | ||
""" | """ | ||
""" | """ | ||
- | 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): | ||
| | ||
Ligne 47: | Ligne 47: | ||
# différents tests : | # différents tests : | ||
- | absc=[] | + | absc = [] |
- | ordo=[] | + | ordo = [] |
- | coef=[0, | + | coef = [0, |
- | xa=-1.04 | + | xa = -1.04 |
while xa < 1.04: | while xa < 1.04: | ||
absc.append(xa) | absc.append(xa) | ||
ordo.append(polyeval(xa, | ordo.append(polyeval(xa, | ||
- | xa=xa+0.01 | + | xa = xa + 0.01 |
plot(absc, | plot(absc, |