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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
' O0 E8 z4 s4 J3 y) R) oNewton[f_,{x_,x0_},e_]:=
Module[{df,x1,y0,y1,k},
! t2 @+ J" C9 `4 H' o2 l6 {df=D[f,x];
, V- U5 b0 _5 a, Z/ g' yy0=f/.x->x0;y1=df/.x->x0;
  F; o) z9 I3 s3 P% l. ~* K9 hx1=x0-y0/y1;
k=1;
While[Abs[x1-x0]>e,
/ x0 T  f) h/ `% i( Nx0=x1;
$ ?+ Y* m8 F, e- w9 l  B5 O0 r# ]- ^y0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
, D9 S5 o3 T( }/ xk=k+1;
1 N6 A9 P: g6 n+ H]
( t; f! [% H: l4 H) \8 r: bPrint["迭代次数为：",k];
# @# o* m) T% H# hPrint["近似根为：",x1];
0 _1 o7 F( D2 r: U8 P: i/ r]

2 b6 x0 h7 H: ?6 C( m8 u5 Z. f8 i但是当我调用它时，会出现这样的提示，如下：

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

0 i. ~' B% G: g' Q& E. Z; i

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
   ! g8 ~) `' \9 W; u3 R
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
     H7 e! d/ G* Q; w0 n
3.   x1 = x0 - y0/y1;2 ^9 r, W. s. |; c9 w8 v\" Z2 _
   
4.   k = 1;
   \" O' e- |) Y- U0 D: ?: v0 E2 ?4 U
5.   While[Abs[x1 - x0] > e, x0 = x1;5 Z\" w, Z8 z) K) q  z
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   + R2 N& V. ~\" L% [+ t* h+ p( p& Y& {
7.     x1 = x0 - y0/y1;: C0 A6 f; V( j: U; y
   
8.     k = k + 1;] Print["迭代次数为：", k];, y! f8 P# m% q, L9 A3 S
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
+ \( t/ P, y! Q1 a) Y5 u4 u: n; w可选参数应该放在一个表里面

哦，明白了，十分感谢！

Agent-Zero

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