clear all
n=input('numero nodi interpolazione=');
a=input('estremo a='); b=input('estremo b=');
x=linspace(a,b,n); y=sin(2*x);
z=a:0.01:b; 
for i=1:length(z), 
    pz(i)=lagr(x,y,z(i));
end
xdisegno=z; ydisegno=sin(2*xdisegno);
plot(xdisegno,ydisegno,'r:',z,pz,'b-.',x,y,'r*')
xlabel('x');ylabel('y');title('interpolazione lagrange sin(2x)');
legend('funzione','pol. interp.','nodi')
errore=max(abs(pz-ydisegno));
fprintf('numero nodi=%3d,   errore=%12.6e\n',n,errore);   