拉格朗日插值函数

function [g n]=Lagrange(A)
n=size(A,1)-1
syms x;
for j=1:n+1
Li(j)=sym('1')
for i=[1:j-1 j+1:n+1]
Li(j)=Li(j)*(x-A(i,1))/(A(j,1)-A(i,1))
end
end
g=Li*A(:,2);
g=expand(g)
end