function y=lagr(x0,y0,x)
%x0,y0為節(jié)點(diǎn)
%x是插值點(diǎn)
n=length(x0);
m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end
標(biāo)簽：
lagr
上傳時(shí)間：
2020-06-09
上傳用戶：shiyc2020
