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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
  ^8 A) f$ t! f  n6 b- |9 Q% ?6 aModule[{df,x1,y0,y1,k},
df=D[f,x];
y0=f/.x->x0;y1=df/.x->x0;
! s' t) ]4 Y4 r+ C9 u3 H$ jx1=x0-y0/y1;
! R" t0 }, b7 Y) [0 f. w; h8 Vk=1;
' q" E, d4 R2 O8 [While[Abs[x1-x0]>e,
5 q9 ^, h- d% \9 Xx0=x1;
y0=f/.x->x0;y1=df/.x->x0;
8 p) K2 d4 H( Q  g' y* m& Z# }x1=x0-y0/y1;
  U: Y4 s% P7 J. ik=k+1;
1 N2 P4 M  H0 g]
Print["迭代次数为：",k];
7 i1 a9 o6 X( }0 D4 H4 J+ \0 {Print["近似根为：",x1];
]
' z; i0 K# \  h) l3 Y0 ?
但是当我调用它时，会出现这样的提示，如下：
8 y# ^% M7 k4 f3 Y9 i# y0 a# R& G
5 p! F- F: f; T) {这是为什么啊？求高手解答，应该怎么修改啊
) e: \, ]6 F5 ^  z; @

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];) Z7 X) _: J4 O/ r1 g8 u: N
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;1 W3 N: g# j; |* O- A
   
3.   x1 = x0 - y0/y1;
   ( y. a6 K5 f* [& E7 A  F
4.   k = 1;\" t6 _$ X\" j2 ^/ ~  L& h+ w
   
5.   While[Abs[x1 - x0] > e, x0 = x1;
   * [( P) Y# M\" C: t9 Z0 E
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;3 R; \$ M) _2 y: {: t\" K/ u2 R
   
7.     x1 = x0 - y0/y1;1 @, f% ]9 P% F3 h
   
8.     k = k + 1;] Print["迭代次数为：", k];3 s4 ^* q  K7 O
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
* F4 w! N! W7 b* B/ k5 y可选参数应该放在一个表里面

哦，明白了，十分感谢！

Agent-Zero

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