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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
! W1 M; D3 N3 |& ~, X' xModule[{df,x1,y0,y1,k},
df=D[f,x];
" J( z7 K1 ]7 U& z) M  ky0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
k=1;
# d* w& j" K, @: k* eWhile[Abs[x1-x0]>e,
x0=x1;
6 P+ Q  r5 w9 uy0=f/.x->x0;y1=df/.x->x0;
4 y4 G0 b% r, O% n# L; q( Yx1=x0-y0/y1;
7 i+ a, U+ q# o, B4 h1 Nk=k+1;
. A, s( K1 B  h( H# p$ O3 Q]
3 Z# Y# Z3 r5 z& f+ qPrint["迭代次数为：",k];
0 f% Y  L7 F+ @7 a: oPrint["近似根为：",x1];
]
0 K2 ?$ H  q/ H
8 M; H3 B) z' w. @" B1 s但是当我调用它时，会出现这样的提示，如下：
5 y% v' ^3 @1 z
# o$ s8 q' g4 I- O7 k9 z这是为什么啊？求高手解答，应该怎么修改啊
( s2 ~9 l  d4 T! A6 I8 x' ^

' h' r3 C1 h9 d8 D2 B

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];; c7 K, K, E: Q! R( d2 W
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;% l& j  U( n- t\" `
   
3.   x1 = x0 - y0/y1;
   # ~4 x5 g+ `; V6 V2 N* I+ R; O
4.   k = 1;: K8 R. J3 f3 ?0 L6 [
   
5.   While[Abs[x1 - x0] > e, x0 = x1;$ B/ h: m1 Z# U2 u$ ]
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   0 b3 T# W) `) `8 j
7.     x1 = x0 - y0/y1;
   5 Y, M, o5 h' R7 h
8.     k = k + 1;] Print["迭代次数为：", k];
   \" l* x) ^9 T' b, J7 `( n' b
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

madio 发表于 2014-8-31 21:59
" j# w* \; J3 B5 B可选参数应该放在一个表里面

4 q( V( C! x' h- _; y哦，明白了，十分感谢！

Agent-Zero

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