// Exemplo do polinomio interpolador usando o metodo de Lagrange.
getf("interpoly.sci");
getf("lagrangep.sci");
// Vamos aproximar o cos(x) no intervalo 0..2pi
tab=[0,%pi/2,%pi,3*%pi/2,2*%pi;
1,0,-1,0,1]
g=interpoly(tab)
// Para construir a funca devemos definir
deff("y=p(x)","y=horner(g,x)");
// plotagem da funcao p com cosseno
x=linspace(0,2*%pi,100)';
plot2d([x,x],[p(x),cos(x)]);
