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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
1 r: }( J! j! l) pNewton[f_,{x_,x0_},e_]:=
Module[{df,x1,y0,y1,k},
df=D[f,x];
y0=f/.x->x0;y1=df/.x->x0;
0 E9 b- X8 {% r  d! xx1=x0-y0/y1;
+ s* i+ f5 P" ?* Q2 G, w" W# Z: ek=1;
2 |) _) E' @: n3 M7 {  QWhile[Abs[x1-x0]>e,
+ T' n4 B# d4 S4 l: bx0=x1;
2 e7 t* ]0 h0 _* k9 Cy0=f/.x->x0;y1=df/.x->x0;
1 e6 i, R, b' ix1=x0-y0/y1;
k=k+1;
1 L5 d$ |& G. Y# Y1 R; O' l4 N]
Print["迭代次数为：",k];
Print["近似根为：",x1];
  @. L' |. r0 R, B5 K$ Q]

但是当我调用它时，会出现这样的提示，如下：
" Y# d+ ~) [' |7 z% j
这是为什么啊？求高手解答，应该怎么修改啊
4 L1 Q% ?5 ]* E- K5 O

2 E" Z1 Q5 w! Q! w3 U+ H) Q

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
   % k. ]0 l6 n2 ^3 T0 A4 R$ U
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
   # \3 f3 z( ?7 q& ^\" u3 m* C8 U+ X
3.   x1 = x0 - y0/y1;
   ' h( h8 M4 G1 c6 t
4.   k = 1;' R2 d: J7 n4 }) S! S9 u% ?! |
   
5.   While[Abs[x1 - x0] > e, x0 = x1;9 _! a\" f$ a9 j7 K
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   & @: E. j' Q$ j/ d) s) E- B
7.     x1 = x0 - y0/y1;. M  B8 a/ y/ L
   
8.     k = k + 1;] Print["迭代次数为：", k];- s5 l\" Y! T9 M. f  e  ^
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
8 U: J3 o! B/ w" [/ c2 N* _可选参数应该放在一个表里面

$ c* b! p9 Q* f" k哦，明白了，十分感谢！

Agent-Zero

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