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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
Module[{df,x1,y0,y1,k},
5 y' B* H6 u' hdf=D[f,x];
: C- n% b4 |# R  P( V  e" Cy0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
k=1;
" V8 T9 q0 ^1 z/ fWhile[Abs[x1-x0]>e,
x0=x1;
y0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
9 K$ a0 M. ?$ Z! X- Q, Nk=k+1;
]
Print["迭代次数为：",k];
2 |$ F6 t! W' V* fPrint["近似根为：",x1];
9 _$ Z' ~, [- L) Q]
9 u, @( y, Q9 G- B' Z, n! ~
; t! j% S. k  k! e0 c但是当我调用它时，会出现这样的提示，如下：

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

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
   3 r6 U4 m# v0 R- L
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
   : H9 [+ g4 G, j  g. {
3.   x1 = x0 - y0/y1;. k* O0 D8 k% E: V- Y
   
4.   k = 1;
   9 m1 }; W: n0 |. K
5.   While[Abs[x1 - x0] > e, x0 = x1;) }5 g1 J/ R\" M& [) k
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;0 J( R1 C7 l2 R& v# L
   
7.     x1 = x0 - y0/y1;
   ) N- Q\" c3 n  |* J, y
8.     k = k + 1;] Print["迭代次数为：", k];
   , J' [# F6 o7 j: G! R# c: Z: `
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
可选参数应该放在一个表里面

哦，明白了，十分感谢！

Agent-Zero

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