数学建模社区-数学中国

标题: 求助，用牛顿切线法编写的一个求根程序，恳请各位大侠指点一下 [打印本页]

---

作者: lxy444    时间: 2014-8-31 15:19
标题: 求助，用牛顿切线法编写的一个求根程序，恳请各位大侠指点一下
我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
2 Y" v$ p6 }3 OModule[{df,x1,y0,y1,k},
df=D[f,x];
/ N8 `, u' L  q' q, P! ^$ X; Py0=f/.x->x0;y1=df/.x->x0;
. W  v8 f" B" Q. a5 C- t5 rx1=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;
2 v. u, q9 H) E! r: O3 }, i]
! M3 V; z- b$ n' f, `Print["迭代次数为：",k];
& y0 `' _9 [! h( z5 [. dPrint["近似根为：",x1];
]

  ]9 a  M( O2 f4 M但是当我调用它时，会出现这样的提示，如下：

5 G% j! c# |# D2 N/ r这是为什么啊？求高手解答，应该怎么修改啊

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

2014-8-31 15:17 上传

点击文件名下载附件

---

作者: lxy444    时间: 2014-8-31 15:22
恳请大神解答！@madio  

---

作者: madio    时间: 2014-8-31 21:59

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
   , g7 T# U9 D2 I! F2 \- ?. f
3.   x1 = x0 - y0/y1;
   % k8 M' ~* o0 Z3 h0 l& f4 d1 o& N3 S
4.   k = 1;
   
5.   While[Abs[x1 - x0] > e, x0 = x1;
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   
7.     x1 = x0 - y0/y1;
   
8.     k = k + 1;] Print["迭代次数为：", k];
   
9.   Print["近似根为：", x1];]

复制代码

可选参数应该放在一个表里面

---

作者: lxy444    时间: 2014-9-1 00:08

madio 发表于 2014-8-31 21:59
5 V- ?) |2 t4 v4 s* D可选参数应该放在一个表里面

; K4 s# |6 N7 v哦，明白了，十分感谢！

---

作者: Agent-Zero    时间: 2014-9-1 10:23
不明觉厉。就是看看。。。

---

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

Powered by Discuz! X2.5