Résultat de défaut de codage d'interpolation de Lagrange

Question

J'ai essayé de calculer la méthode d'interpolation de Lagrange dans MATLAB pour y=x^2+xe^(x). J'ai écrit le code suivant:

clc 
clear 
close all 

x0=4.7; 
n=10; 
x=linspace(0,5,n); 
y=x.^2+x.*exp(x); 

syms t 
L=sym(ones(1,n)); 
P_x=sym(0); 
for i=1:n 
    for j=1:n 
     L_improcess=(t-x(j))/(x(i)-x(j)); 
     if(i==j) 
      continue 
     end 
     L(i)=L(i)*L_improcess; 
     P_x=y(i)*L(i)+P_x; 
    end 
end 
P=double(subs(P_x,t,x0)); 
disp(['Lagrange interpolation: P= ',num2str(P)]) 
disp(['the real value from original function is:' num2str(x0^2+x0*exp(x0))]) 

donc le résultat pour x0=4.7 était:

Lagrange interpolation: P= 20195.8626 
the real value from original function is:538.8417 

Je me demande comment la différence entre deux résultats (les deux de la même doit être presque) Lagrange procédé d'interpolation pour f(x) est comme ceci:

plus informat ion sur l'interpolation de Lagrange est disponible here.

Answer 1

Cette ligne:

P_x=y(i)*L(i)+P_x; 

devrait être dans la boucle sur i, non j.