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

lxy444

我想用mathematica编写一个求根函数，用牛顿切线法，我编的函数如下：
Newton[f_,{x_,x0_},e_]:=
2 p* c0 r9 m6 Z3 |3 r3 c/ mModule[{df,x1,y0,y1,k},
df=D[f,x];
% V9 B' n/ \, {2 ]" _& ]0 B* Zy0=f/.x->x0;y1=df/.x->x0;
: G+ _) F  b3 D7 Mx1=x0-y0/y1;
3 b& @1 Q* i  D! i1 ]k=1;
While[Abs[x1-x0]>e,
. }6 N; ?4 `4 x/ D: E. z2 S, C! J' nx0=x1;
% _% m# {$ C. a/ W3 H  Qy0=f/.x->x0;y1=df/.x->x0;
x1=x0-y0/y1;
8 ~% w) h( k; I# Rk=k+1;
]
. U7 e; p' d. U$ hPrint["迭代次数为：",k];
$ B4 m% k- ]* j; ]9 HPrint["近似根为：",x1];
% j" w! o' O3 G# u, b" |]
7 k9 c1 V! v6 Y- G
. ~* |: n5 Z# c4 P9 O, |; T( Q& z+ D但是当我调用它时，会出现这样的提示，如下：

- W3 a1 w8 h! I. [0 }这是为什么啊？求高手解答，应该怎么修改啊

' }4 c% X0 q; Q

lxy444

恳请大神解答！@madio  

madio

1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];' k5 S6 D7 `+ J! p  q3 A
   
2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;/ r3 p( y+ t6 S' T
   
3.   x1 = x0 - y0/y1;
   * q' F1 `0 y, ?$ I, J5 j
4.   k = 1;+ \5 Z+ H( n4 |5 G
   
5.   While[Abs[x1 - x0] > e, x0 = x1;2 W0 l; t' `7 G4 |2 R
   
6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
   + x( c- W+ c# B$ D: y
7.     x1 = x0 - y0/y1;
   ' w9 e/ @1 Z6 |4 `% r' ~$ U
8.     k = k + 1;] Print["迭代次数为：", k];# |8 h0 \) K( ~# A
   
9.   Print["近似根为：", x1];]

复制代码

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

lxy444

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

8 d8 }  [7 S# b- n哦，明白了，十分感谢！

Agent-Zero

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