
// On teste la fontion Newton sur l'exemple de l'exo 1 du TD1. On doit obtenir le même résultat qu'avec Lagrange (car il y a unicité du polynome d'interpolation)

T=[-10:0.1:10]
R1=newton([0,1,2],[1,2,0],T)
R2=lagrangeP([0,1,2],[1,2,0],T)

scf(0)
clf()
plot2d(T,R1)

scf(1)
clf()
plot2d(T,R2)

// Un autre test

T= linspace(0,2*%pi,100)

X=[1,2,3]

scf(2)
clf()
plot2d(T,sin(T),2)
plot2d(T,lagrangef(X,sin,T))

scf(3)
clf()
plot2d(T,sin(T),2)
plot2d(T,newton(X,sin(X),T))


// Comme on a deux méthodes différentes pour faire un même calcul on peut comparer leur vitesse d'execution

T= linspace(0,2*%pi,100)

X=[1:10]

timer();lagrangeP(X,sin(X),T);t1=timer()

disp(t1)


timer();newton(X,sin(X),T);t2=timer()

disp(t2)

// La première méthode est bien plus rapide ce qui n'est pas étonant car la seonde méthode n'a pas été optimisée (cf remarque fin du fichier exo3.sci)
