lx.m

来自「Lagrange插值、Newton插值函数等插值函数」· M 代码 · 共 19 行

function [lx]=lx(x,X,y)



lx=zeros(1,size(x,2));
for e=1:size(x,2)
 for i=1:size(X,2)
     j=1;
     k=1;
   for l=1:size(X,2)
        if i==l
        else 
           j=j*(x(1,e)-X(1,l));
           k=k*(X(1,i)-X(1,l));
        end
   end
   lx(1,e)=lx(1,e)+j/k*y(1,i);
 end
end