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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
Module[{df,x1,y0,y1,k},
df=D[f,x];
y0=f/.x->x0;y1=df/.x->x0;
2 D) I% Y9 n8 k) R" A# G! v) nx1=x0-y0/y1;
) {2 O. o0 l! a! J6 ~k=1;
1 o3 R- M; F1 O, E- \. s, oWhile[Abs[x1-x0]>e,
x0=x1;
; w3 B, O7 _/ h, E3 X" oy0=f/.x->x0;y1=df/.x->x0;
0 k  q+ ]6 M$ y6 |& J1 Jx1=x0-y0/y1;
k=k+1;
]
" }  g3 V7 r! y4 ]) KPrint["迭代次数为：",k];
; q& r' |# b/ RPrint["近似根为：",x1];
$ k/ i- w- I+ R4 a" }3 U  v. c, V8 e]

但是当我调用它时，会出现这样的提示，如下：
! m, B* V4 F* q0 B  s% C
这是为什么啊？求高手解答，应该怎么修改啊
/ X2 g2 J6 k0 e/ [- x
/ L" ^8 a4 N* f6 T3 ]! F

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
   0 K+ w) \: A$ ~. Q
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
   \" i5 Z4 R4 n+ W, {6 A  S
3.   x1 = x0 - y0/y1;
   7 a1 |3 z) o- V  d) H
4.   k = 1;0 C# t6 S5 z# Q: L+ j; U+ w
   
5.   While[Abs[x1 - x0] > e, x0 = x1;
   ' @0 q+ y7 G# _# s4 v) m' n) M
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;( \& F+ Z- @9 U) K9 ~
   
7.     x1 = x0 - y0/y1;
   9 N. S$ @$ r2 i- o7 q$ B' w
8.     k = k + 1;] Print["迭代次数为：", k];
   ; V9 j0 C. N+ g
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

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

哦，明白了，十分感谢！

Agent-Zero

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