полиномиальная интерполяция методом Лагранжа

Question

Вот мой код для интерполяции Лагранжа, я сделал 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