这是拉格朗日插值的代码，我做了3个拉格朗日函数和多项式和插值，我在多项式函数和插值函数中的多项式中进行了拉格朗日调用 . 它没有用，我看不到我的近似图 . 可以帮助我

function interpolation
xx =0:pi/4:pi;
F =[0,0.078,0.57,1.64,3.14];
n=size(xx)
x1=0:pi/12:pi;
y=@(x) x-sin(x)
for i=1:n
 y=polynome(n,y,x1(i),xx)

 plot(xx,y,'r');

end
end


function R=polynome(n,F,x,xx)
R=0;
for i=1:n

R=R+F(i)*Lagrange(n,x,xx,i)
end
end

function L=Lagrange(n,x,xx,i)
  p=1;
  for j=1:n
  if j~=i 
  p = p*(x-xx(j))/(xx(i)-xx(j)) 
  end
  end
  L=p;
 end
matlab