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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
Module[{df,x1,y0,y1,k},
! o& Y' s- w2 I0 l/ n% b+ _- w; Kdf=D[f,x];
y0=f/.x->x0;y1=df/.x->x0;
' R, e$ [# O, ^0 w3 kx1=x0-y0/y1;
k=1;
5 |! M, C" u3 ]# _" \) C8 X5 k) oWhile[Abs[x1-x0]>e,
0 ^! \: G8 B* a% ?, ~x0=x1;
: o( `  `. O8 j" y& Dy0=f/.x->x0;y1=df/.x->x0;
: d/ V# l4 @0 y0 Vx1=x0-y0/y1;
k=k+1;
9 o) y' C1 t+ }5 }8 A; z5 `% P]
Print["迭代次数为：",k];
Print["近似根为：",x1];
7 `1 k; q; }0 Y* t7 @/ u]
# p' z1 v' a+ x/ K/ _* _
但是当我调用它时，会出现这样的提示，如下：
% [# c1 u5 H  [! i% j, G
这是为什么啊？求高手解答，应该怎么修改啊


1 [0 Z6 k/ J1 \

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];4 c% h, c/ k; ?\" U% r+ U
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;& S4 [3 |: n9 [
   
3.   x1 = x0 - y0/y1;4 T- X6 x5 w: o5 m4 r7 E' t8 F6 B
   
4.   k = 1;$ Z+ P+ q5 w( I  J6 L
   
5.   While[Abs[x1 - x0] > e, x0 = x1;# c3 Z& s6 Q  h) |& T2 Z7 E
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   2 O* D. Q  w6 B& ]% ]\" u1 I
7.     x1 = x0 - y0/y1;
   ) G2 X4 c+ R\" _- a$ S. Q9 x
8.     k = k + 1;] Print["迭代次数为：", k];
   0 t& R0 ~% S3 w4 T
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
, g9 n! `; S0 h3 j% t可选参数应该放在一个表里面

哦，明白了，十分感谢！

Agent-Zero

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