24小时热门版块排行榜

返回列表

tzh0724

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 194.2
帖子: 89
在线: 13.9小时
虫号: 1026238
注册: 2010-05-22

[交流] 【求助】求助程序已有2人参与

lagrange.m
%lagrange insert
function y=lagrange (x0, y0, x)
n=length (x0);m=length (x);
for i=1:m
x=x (i);
s=0.0;
for k=1:n
      p=1.0;
      for j=1:n
         if j~=k
            p=p*(x-x0(j))/(x0(k)-x0(j));
         end
      end
      s=p*y0 (k)+s;
end
y(i)=s;
end
这个程序不知道为什么运行不了，求各位帮帮忙，先谢谢了。

回复此楼

» 猜你喜欢

1楼 2010-07-22 16:07:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

仿真EPI: 3
应助: 157 (高中生)
贵宾: 5.735
金币: 40609
散金: 2394
红花: 63
沙发: 4
帖子: 4046
在线: 1330.2小时
虫号: 653607
注册: 2008-11-13
性别: GG
专业: 机械动力学
管辖: 仿真模拟

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+1):鼓励应助。 2010-07-23 18:47:46

引用回帖:

Originally posted by tzh0724 at 2010-07-22 16:07:32:
lagrange.m
%lagrange insert
function y=lagrange (x0, y0, x)
n=length (x0);m=length (x);
for i=1:m
x=x (i);
s=0.0;
for k=1:n
      p=1.0;
      for j=1:n
         if j~= ...

程序错在第5行
x=x (i);

是不是该 x=x0(i); ?