MATLAB 牛顿插值代码程序
我的程序一直出错，我就差列不出牛顿插值多项式了，麻烦各位高手出手相助！！
function f=agui_newtonInsert(x0,y0,x)
%x0为节点向量，y0为节点上的函数值，x为插值点，f返回插值
n=length(x0);
%均差对称性 （6.3.6）
s=0.0;
g=0.0;
for k=1:n
p=1.0;
for i=1:k
for j=1:k
if j~=i
p=p*(x0(i)-x0(j));
end
s=s+(y0(i)/p);
q=1.0;
q=q*(x-x0(k));
g=s*q+g;
end
end
end
f=g+y0(1);
f=vpa(f,6);
end
[解决办法]
vpa是什么函数？
[解决办法]
你这差商求得不对吧，怎么只有x的没y的
[解决办法]
还有函数结束的end会导致函数没法调用
[解决办法]
擅长matlab。
vpa， 是 matlab， 高精度显示
如： vpa(pi,40)
3.141592653589793238462643383279502884197
matlab的 end，有多种用途， 一种是用来划分代码块， function end， if end， for end..
