LAYOUDE.MesTravauxEns

Il s'agit d'une méthode d'approximation

Pour plus d'information sur Newton et ses différents travaux faites un petit tour sur Google

<Vous ne serez ignorant que pour quelques instants>

function y=newton(X,Y,x0)
 
  % calcul des Cn
   j=length(X);
   C(1)=Y(1);
   for i=2:j
       k=i;
       denum=1;
       while k>1
           denum= denum*(X(k)-X(k-1));
           k=k-1;
       end
       C(i)=(Y(i)-Y(i-1))/denum;
   end
    
   % calcul de la foction au point x0
    
   U(1)=C(j);
   for i=2:j
       U(i)=C(j-i+1)+U(i-1)*(x0-X(j-i+1));
   end
   y=U(j);
   return

Merçi dde toujours critiquer. "Aucune oeuvre de l'esprit n'est parfaite"