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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
/ k' R9 P8 ^1 Y2 }! M5 s9 ~1 mNewton[f_,{x_,x0_},e_]:=
4 G! l  N9 f- `$ p8 Y$ l$ [  cModule[{df,x1,y0,y1,k},
% Q% s2 b: x3 s/ T( z4 Sdf=D[f,x];
' `) k" e( }3 m! c% Z/ ^y0=f/.x->x0;y1=df/.x->x0;
" y8 M$ V: \( L, I7 W9 k& vx1=x0-y0/y1;
* a* F, m9 g' w; ^) nk=1;
7 B+ U% c' j. j8 lWhile[Abs[x1-x0]>e,
3 P9 w, B; q9 L* z8 ax0=x1;
y0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
k=k+1;
( v8 G% B; b9 @9 i]
7 A- y- t4 a" g' J) bPrint["迭代次数为：",k];
Print["近似根为：",x1];
" M/ H9 p( Q6 ~]
/ p3 N! a2 g* ]- M6 a8 {0 }0 v
$ `! s/ Q4 Q- D  S但是当我调用它时，会出现这样的提示，如下：

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

, c/ l- C8 p. P4 s# y

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];7 i0 J1 _5 s/ a( g5 T7 u9 K% n
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;0 ~/ b0 Q' {2 L2 D
   
3.   x1 = x0 - y0/y1;1 y; t0 i$ [* m1 N1 H- F
   
4.   k = 1;
   ! ~9 R# {5 V  [# U, a0 P
5.   While[Abs[x1 - x0] > e, x0 = x1;) ^0 d6 {6 B: [/ p, D1 ?' ?1 e
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
     b1 a3 i5 y8 k5 x' H& ?
7.     x1 = x0 - y0/y1;; s4 C, p7 f8 ^. n! x7 x
   
8.     k = k + 1;] Print["迭代次数为：", k];/ ]) D8 A4 S, H7 X
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
' t! P" Z3 s/ \" Y# D1 {- Y' M可选参数应该放在一个表里面

1 w; o; ^3 u! [7 r  t哦，明白了，十分感谢！

Agent-Zero

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