la.asv

来自「fit programme for a serials of data」· ASV 代码 · 共 15 行

function y1=L(x,y,x1)
% input x is a vector that contains a list of abscissas
%       y is a vector that cntains a list of ordinates
%       x1 is the number at where the interpolation is to be obained
% output X is the la interpolation at x1
y=0;
for k=0:n
   p=1;
   for j=0:n
      if j~=k
         p=p*(x1-x(j))/(x(k)-x(j)); 
      end
   end
   
end