# -*- coding: utf-8 -*-
"""
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
    cf. http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Ruffini-Horner
    """
    n = len(a)-1 # n = ordre du polynôme
    p = 0.
    for i in range(n,-1,-1):
        p = p*x + a[i]
    return p

# utilisation des objets de la librairie numpy
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

# création des graphes de tous les polynomes

for pol in cheby:
    plot(x,polyeval(x,pol))

axis([-1,1,-1,1]) # xmin, xmax, ymin, ymax   
title('Polynomes de Tchebyshev')
legend()
show()