simple.m

来自「牛顿插值法」· M 代码 · 共 11 行

n=10;
A1=sparse(1:n,1:n,-2,n,n);
A2=sparse(1:n-1,2:n,1,n,n);
A=A1+A2+A2';
B1=ones(1,8);
B2=-1.5*B1;
B3=[-0.5];
B4=cat(2,B3,B2,B3);
B=(B4)';
x=A\B;
x