lagrange.m

来自「一些数值计算用的MATLAB代码,适合初学MATLAB和数字计算的人使用,根据M」· M 代码 · 共 20 行

n=input('n=');         %输入要做的次数
n=n+2;
 x=input('x[]=');      %输入矩阵
 y=input('y[]=');      %输入矩阵
u=input('要求值：');
v=0;
j=1;
while j~=n
    l=1;
    i=1;
    while i~=n
        if i~=j
           l=l*(u-x(i))/(x(j)-x(i));
        end
        i=i+1;
    end
    v=v+l*y(j);
    j=j+1;
end
u,v