function z=lagr(x,y)
// z=lagr(x,y)
// 
// Cria polinomio de lagrange
// Entradas:  x = pontos a interpolar x = [x0 x1 ... xn] 
//            y = [f(x0) f(x1) ... f(xn)]
// Saida:     
//            z = [ a0 a1 ...] tais que 
//            P(x) = a0 + a1*x + ...  eh o pol interpolador de Lagange
//
// Eduardo. 14/03/10


c=size(x,2);
M=[];
for k=1:c,

 l = 1;
 for j=1:c,
 if j~=k,
   l = inv(x(k)-x(j))*convol([-x(j) 1],l);
 end
 end
M = [ M; y(k)*l];

end

z = sum(M,1);
endfunction




