查看: 3022|回复: 4

求助，用牛顿切线法编写的一个求根程序，恳请各位大侠指点一下

字体大小: 正常放大

lxy444

712 主题	16 听众	1万积分

升级 0%

TA的每日心情

	开心 2026-6-3 23:37

签到天数: 3994 天

[LV.Master]伴坛终老

群组: 2013年数学建模国赛备

电梯直达

1^#

发表于 2014-8-31 15:19 |只看该作者 |正序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
Module[{df,x1,y0,y1,k},
df=D[f,x];
y0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
k=1;
While[Abs[x1-x0]>e,
x0=x1;
y0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
k=k+1;
]
Print["迭代次数为：",k];
Print["近似根为：",x1];
]

但是当我调用它时，会出现这样的提示，如下：

这是为什么啊？求高手解答，应该怎么修改啊

1.PNG (19.82 KB, 下载次数: 247)

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

Agent-Zero

0 主题	11 听众	102 积分

升级 1%

TA的每日心情

	开心 2014-10-17 07:39

签到天数: 22 天

[LV.4]偶尔看看III

自我介绍: 建模苦孩子。

5^#

发表于 2014-9-1 10:23 |只看该作者 |招呼Ta 关注Ta

不明觉厉。就是看看。。。

我也说一句

发表

使用道具举报

lxy444

712 主题	16 听众	1万积分

升级 0%

TA的每日心情

	开心 2026-6-3 23:37

签到天数: 3994 天

[LV.Master]伴坛终老

群组: 2013年数学建模国赛备

4^#

发表于 2014-9-1 00:08 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

madio 发表于 2014-8-31 21:59 ) i) A- G" f+ {
可选参数应该放在一个表里面

哦，明白了，十分感谢！

我也说一句

发表

使用道具举报

madio

3万主题	1312 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

3^#

发表于 2014-8-31 21:59 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
1 N7 L7 \6 y\" o6 z3 D2 c* S
y0 = f /. {x -> x0; y1 = df} /. x -> x0;
+ p- B1 L3 _$ _
x1 = x0 - y0/y1;
; ]: i' c$ \7 C6 p1 [! l: x* G( e
k = 1;
- v0 E/ c/ O4 P2 _/ t0 F
While[Abs[x1 - x0] > e, x0 = x1;- D0 A( e8 ~5 B0 V; v! b2 R& f0 k
y0 = f /. x -> x0; y1 = df /. x -> x0;& t; X# A& Z5 B2 Q% b
x1 = x0 - y0/y1;
% M* O- g! a7 q+ }6 u
k = k + 1;] Print["迭代次数为：", k];* s+ P6 s/ o\" _0 `5 }
Print["近似根为：", x1];]