数学建模社区-数学中国

标题: 埃尔米特 [打印本页]

作者: suntangshi 时间: 2012-11-28 11:46
标题: 埃尔米特
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

请问这个程序到底出了啥问题吗？如何调用它在命令窗口使用？

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)