LAYOUDE.MesTravauxEns

Pour plus d'information sur l'approximation de Lagrange faites une petite recherche sur Google

<Vous ne serez ignorant que pour quelques minutes>

function y=lagrange(X,Y,x0)
     
  % calcul des fonctions de lagrange
     j=length(X);
     result=0;
      for i=1:j
          num=1;
          denum=1;
          k=1;
           while k<j+1
               if k~=i
                   num=num*(x0-X(k));
                   denum=denum*(X(i)-X(k));
               end  
               k=k+1;
           end
           L(i)=num/denum;                 
      end
      % calcul de la fonction en x0
      for i=1:j
          result=result+(L(i)*Y(i));
      end
      y=result;
     
    
     return

merçi de toujours critiquer: aucune oeuvre de l'esprit n'est parfaite