Différences

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

--- teaching:progappchim:polynomes-9 [2012/11/16 07:15] – créée villersd
+++ teaching:progappchim:polynomes-9 [2017/02/24 11:45] (Version actuelle) – villersd
@@ Ligne 1: / Ligne 1: @@
 ====== Polynômes : graphe multiple fonctions polynomiales ======
-<sxh python; title : poly09-plot-numpy-chebby.py>
+<code python poly09-plot-numpy-chebby.py>
 # -*- coding: utf-8 -*-
-# graphe multiple de polynômes de Tchebyshev
+"""
-# cf. http://fr.wikipedia.org/wiki/Polyn%C3%B4me_de_Tchebychev
+graphe multiple de polynômes de Tchebyshev
+cf. http://fr.wikipedia.org/wiki/Polyn%C3%B4me_de_Tchebychev
+"""
 from pylab import *   # librairie graphique (Matplotlib)
 def polyeval(x,a):
-    """application de l'agorithme de Horner
+    """
+    application de l'agorithme de Horner
     cf. http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Ruffini-Horner
-    variante de l'agorithme utilisant un while...
     """
-    n=len(a)
+    n = len(a)-1 # n = ordre du polynôme
-    result=0
+    p = 0.
-    while n>=1:
+    for i in range(n,-1,-1):
-        result=result*x+a[n-1]
+        p = p*x + a[i]
-        n-=1
+    return p
-    return result
-# utilisation des objets numpy
+# utilisation des objets de la librairie numpy
-x=arange(-1.,1.0001,0.01)   # avec des tableaux numpy, on peut évaluer les polynome
+x = arange(-1.,1.0001,0.01)   # avec des tableaux numpy, on peut évaluer les polynome
 # en une fois sur tous les éléments du tableau
-cheby=[[1],[0,1],[-1,0,2],[0,-3,0,4]]  # liste de polynomes
+cheby = [[1], [0,1], [-1,0,2], [0,-3,0,4]]  # liste de polynomes
 # création des graphes de tous les polynomes
@@ Ligne 33: / Ligne 34: @@
 legend()
 show()
-</sxh>
+</code>
 [[polynomes-10|Des fonctionnalités supplémentaires, en vue de visualiser un nombre quelconque de polynômes...]]