solve.m

来自「Lagrange插值、Newton插值函数等插值函数」· M 代码 · 共 14 行

function m=solve(u,v,d)

xs=zeros(size(u,2)+2);
U=[u,0];
V=[0,v];
for i=1:size(xs,2)
    xs(i,i)=2;
end
for i=1:size(xs,2)-1
    xs(i,i+1)=V(i);
    xs(i+1,i)=U(i);
end
D=[0,d,0]';
m=xs\D;