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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
" I8 i1 B' t5 u* C6 G/ K+ \Module[{df,x1,y0,y1,k},
' K# e8 Q4 U9 r2 Gdf=D[f,x];
: H7 e4 Y6 w9 D; y9 x1 w! d5 Ty0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
k=1;
1 f, K# n' @# Z1 IWhile[Abs[x1-x0]>e,
x0=x1;
" t/ R4 _1 _% c7 sy0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
, m# ?4 c! P! c! S- x3 ek=k+1;
6 y" s+ V1 Q- @) Q1 N/ ]9 O]
% [5 T: O! l6 q( |( u3 }! ~Print["迭代次数为：",k];
Print["近似根为：",x1];
  T% [" P. H& r* [6 C! i  a! Q]

但是当我调用它时，会出现这样的提示，如下：
% ^8 J# p* N* ]5 G
这是为什么啊？求高手解答，应该怎么修改啊
; ^: |7 E$ j. S$ ?7 k2 P& u. u) G

: _& s: Y$ k7 _8 A1 K9 I. e# E& r" j

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];& C  w7 J, \\" @
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;# D1 `, k4 A; E8 M+ q; U
   
3.   x1 = x0 - y0/y1;9 R, r. a* r( J\" [) l
   
4.   k = 1;
   9 g5 H+ k( m; Q2 b- a$ V; b9 Q
5.   While[Abs[x1 - x0] > e, x0 = x1;
   ) z. W( a! i: K) b# i9 M
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   ) F. j: v2 X$ b& |4 M
7.     x1 = x0 - y0/y1;' m0 E4 q' n+ u8 f' w
   
8.     k = k + 1;] Print["迭代次数为：", k];  M/ J' v) A9 `# Q
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
) l1 S" S& i/ g& k) b4 r; g可选参数应该放在一个表里面

哦，明白了，十分感谢！

Agent-Zero

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