埃尔米特
function y=hermite(x0,y0,y1,x)n=length(x0),m=length(x);
for k=1:m
yy=0.0;
for i=1:n
h=1.0;
a=0.0;
for j=1:n
if j~= i
h=h*((x(k)-x0(j))/(x0(i)-x0(j)))^2;
a=a+1/(x0(i)-x0(j));
end
end
yy=yy+h*((x0(i)-x(k))*(2*a*y0(i)-y1(i))+y0(i));
end
y(k)=yy;
end
请问这个程序到底出了啥问题吗？如何调用它在命令窗口 使用？
页:
[1]