Matlab实现插值算法（二）

要求

代码

运行结果

#要求

#代码
下面展示一些 内联代码片。

function f=lagrange(x0,y0,x)
%x0为节点向量，y0为节点上的函数值，x为插值点，f返回插值
n=length(x0);
m=length(x);
format long
s=0.0;
for k=1:n
    p=1.0;
    for j=1:n
        if j~=k
            p=p*(x-x0(j))/(x0(k)-x0(j));
        end
    end
    s=p*y0(k)+s;
end
f=s;
end

需在命令行键入要求参数

x0=[0.10 0.15 0.25 0.40 0.50 0.57 0.70 0.85 0.93 1.00];
y0=[0.904837 0.860708 0.778801 0.670320 0.606531 0.565525 0.496585 0.427415 0.394554 0.367879];
lagrange(x0,y0,0.63) An highlighted block
var foo = 'bar';

#运行结果