lagrange Algoritmo

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function r=lagran(A, Y,L,)
 x=sym('x');
 format long
 Z=A;
 p=length(A);
 J=1;
 Q=1;
 R=0;
 for k=1:p
 for n=1:p
 J=J*(x-(A(n)));
 Q=Q*(Z(k)-A(n));
 end
 end
 L(k)=J/Q;
 R=R+(Y(k)*L(k));
 J=1;
 Q=1;
end