拉格朗日方法的多项式插值

Question

这是我的Lagrange插值代码，我做了3个lagrange函数以及多项式和插值，在多项式函数中进行了lagrande调用，在插值函数中进行了多项式。 它没有解决，我看不到我的近似图。能帮我吗

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