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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
2 U9 j( ^! m$ y6 P0 tNewton[f_,{x_,x0_},e_]:=
( ~" i4 e2 N% c# Y" {Module[{df,x1,y0,y1,k},
: n* d& E8 f: P6 odf=D[f,x];
+ n% k* z) S3 P9 u7 Qy0=f/.x->x0;y1=df/.x->x0;
! ^0 ?) s% O* T' c' a. g+ rx1=x0-y0/y1;
5 p7 @& b( b- V& {k=1;
( S) B/ }2 n) d* s/ H5 y7 K8 C3 B5 VWhile[Abs[x1-x0]>e,
x0=x1;
y0=f/.x->x0;y1=df/.x->x0;
6 T! J4 G1 k8 m, C+ q& d- M0 G; _x1=x0-y0/y1;
k=k+1;
]
5 A& p* }0 j2 G( V8 ?& O( }0 PPrint["迭代次数为：",k];
Print["近似根为：",x1];
]
! E) H7 K+ P: ~! {/ S, D* Z0 k" `
但是当我调用它时，会出现这样的提示，如下：

2 Y! C: u; A/ X) w* h9 h: d这是为什么啊？求高手解答，应该怎么修改啊


4 u/ O! C0 ]0 L) ^

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];$ y: B1 k% z% u) O
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
   . I6 \0 M1 B0 u( h7 k4 g
3.   x1 = x0 - y0/y1;
   ( ]9 h  B) ]/ v# e. ]7 s
4.   k = 1;7 [% m. z0 }/ J0 J! v
   
5.   While[Abs[x1 - x0] > e, x0 = x1;% `7 n8 w4 `, T  W9 G
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   ! M( P5 c4 y% |
7.     x1 = x0 - y0/y1;; N* O6 ]: A$ F1 t6 g\" t2 |( M
   
8.     k = k + 1;] Print["迭代次数为：", k];
   % P8 P* B- }2 u3 W( L: ?* B
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
5 G! H# U5 b' n$ b( h! d9 |% e0 g可选参数应该放在一个表里面

& A4 P$ P4 H0 p, U3 N哦，明白了，十分感谢！

Agent-Zero

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