王指导的博客

牛顿插值法求差值 %求P(x) for i=1:1:m a=1; b=f(1,1); for j=2:1:n a=a*(xx(i)-x(j-1)); b=b+a*f(j,j); end yy(i)=b; end

spline 第一边界条件 其中用列主元消去法求解线性方程组