要求
代码
运行结果
#要求
#代码
下面展示一些 内联代码片
。
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';
#运行结果