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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
$ |' S6 x# t; t# KNewton[f_,{x_,x0_},e_]:=
3 W3 ?( _' C- uModule[{df,x1,y0,y1,k},
$ a$ ?+ E2 h" Z3 m0 |df=D[f,x];
y0=f/.x->x0;y1=df/.x->x0;
8 |+ y0 g3 |0 z0 y# Y$ Xx1=x0-y0/y1;
% f+ r( A9 a6 b6 c- w  t; mk=1;
While[Abs[x1-x0]>e,
x0=x1;
# G: S& j4 i' |5 ty0=f/.x->x0;y1=df/.x->x0;
! Y& C9 S* \6 x; d6 h( G; \* Jx1=x0-y0/y1;
k=k+1;
, U5 p6 g4 _- y9 e/ G/ r' X: k' q]
Print["迭代次数为：",k];
Print["近似根为：",x1];
]

但是当我调用它时，会出现这样的提示，如下：
# I: W8 B/ I) b
这是为什么啊？求高手解答，应该怎么修改啊
$ P$ H% k0 A: i, y. ?5 k. g, V& h
2 ~/ V% R! ?  r, g) m

Agent-Zero

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

lxy444

madio 发表于 2014-8-31 21:59
; g! K8 z( J- C7 a7 a可选参数应该放在一个表里面

6 B$ |" H) v$ |! g3 y6 {  S' q哦，明白了，十分感谢！

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];& P* k( K' q: A2 c# K' e
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
   2 i! B; S1 W) A8 Z
3.   x1 = x0 - y0/y1;
   . U- A1 \% \3 P2 T. X9 M( ?
4.   k = 1;% Y  x1 f+ `$ S
   
5.   While[Abs[x1 - x0] > e, x0 = x1;
   / V- i% \! X- B6 t5 r
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;0 A* ]( T& c# z$ r( M; |\" v5 Y
   
7.     x1 = x0 - y0/y1;
   2 k, u' P. {+ m
8.     k = k + 1;] Print["迭代次数为：", k];  n9 I/ t\" J; ?\" i# [7 {8 @2 `* Z
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

恳请大神解答！@madio