QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2833|回复: 4
打印 上一主题 下一主题

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

[复制链接]
字体大小: 正常 放大
lxy444        

712

主题

16

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2025-10-12 00:09
  • 签到天数: 3988 天

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    跳转到指定楼层
    1#
    发表于 2014-8-31 15:19 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    我想用mathematica编写一个求根函数,用牛顿切线法,我编的函数如下:; `8 u1 _+ I3 x* f) q! @
    Newton[f_,{x_,x0_},e_]:=
    " X" F! c+ g: r) W$ m6 s# Y6 \3 F! eModule[{df,x1,y0,y1,k},# J5 w4 Z7 g8 f, l( t5 @% ^" }
    df=D[f,x];
    & P5 b' _: @% b) py0=f/.x->x0;y1=df/.x->x0;, s( {. C+ Y/ _( b' C" o3 I8 [
    x1=x0-y0/y1;6 n" h* r, u# y& t! |
    k=1;4 [8 r: }* c) ?6 n- f# g
    While[Abs[x1-x0]>e,. j9 |1 [  e: I( T# \1 c" \8 ?
    x0=x1;
    ; v, D- F: a; t) H* ^/ r& Jy0=f/.x->x0;y1=df/.x->x0;
    : z' O0 _- F1 R, y  t# Qx1=x0-y0/y1;" G1 H$ ?6 U- P) E- d  ^+ ], L
    k=k+1;
    1 [% h1 v. k( c9 @7 i7 v]! [. z3 U+ R* z* y0 a2 a
    Print["迭代次数为:",k];' f3 Y8 ?, I# H
    Print["近似根为:",x1];
    4 j% r! n5 V  D: I% \. C]
    8 I9 e* Y+ O& z* K
    ) }3 J& }8 v: x1 e但是当我调用它时,会出现这样的提示,如下:9 u- G. r& e( e# u4 i' O0 \1 W
    ' G+ _5 U, a% _5 a. j, d
    这是为什么啊?求高手解答,应该怎么修改啊
    4 D- a- j7 _0 |
    * z0 b: b8 o# D# {  z# }+ q0 o7 y7 i" v5 Q9 }, e

    1.PNG (19.82 KB, 下载次数: 225)

    1.PNG

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    lxy444        

    712

    主题

    16

    听众

    1万

    积分

    升级  0%

  • TA的每日心情
    开心
    2025-10-12 00:09
  • 签到天数: 3988 天

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    回复

    使用道具 举报

    madio        

    3万

    主题

    1311

    听众

    5万

    积分

  • TA的每日心情
    奋斗
    2024-7-1 22:21
  • 签到天数: 2014 天

    [LV.Master]伴坛终老

    自我介绍
    数学中国站长

    社区QQ达人 邮箱绑定达人 优秀斑竹奖 发帖功臣 风雨历程奖 新人进步奖 最具活力勋章

    群组数学建模培训课堂1

    群组数学中国美赛辅助报名

    群组Matlab讨论组

    群组2013认证赛A题讨论群组

    群组2013认证赛C题讨论群组

    1. Newton[f_, {x_, x0_}, e_] := Module[{df, x1, y0, y1, k}, df = D[f, x];
      0 N, w0 `4 g$ Y& T1 e7 _  O
    2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
      * l0 d' t2 Y9 i% d2 Z2 f\" J3 R- D( h
    3.   x1 = x0 - y0/y1;
      : j# p% [2 E: [
    4.   k = 1;
      $ A& |0 w( e7 Y) o& u5 P
    5.   While[Abs[x1 - x0] > e, x0 = x1;
      $ ]+ T! N! A* z+ f
    6.     y0 = f /. x -> x0; y1 = df /. x -> x0;+ W$ J9 H( {9 ?5 D$ |8 X
    7.     x1 = x0 - y0/y1;
      \" Q- _- V, n. S' I- a5 I1 T
    8.     k = k + 1;] Print["迭代次数为:", k];
      ! d+ U* [- ^( Z$ E' T7 e4 m, d
    9.   Print["近似根为:", x1];]
    复制代码
    可选参数应该放在一个表里面
    数学建模社会化
    回复

    使用道具 举报

    lxy444        

    712

    主题

    16

    听众

    1万

    积分

    升级  0%

  • TA的每日心情
    开心
    2025-10-12 00:09
  • 签到天数: 3988 天

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    madio 发表于 2014-8-31 21:59 ; ~( U6 c+ |5 z7 `/ r* x" A5 \5 W6 y
    可选参数应该放在一个表里面

    / ?5 I) Y: q7 `+ h# ?% i哦,明白了,十分感谢!
    回复

    使用道具 举报

    0

    主题

    11

    听众

    102

    积分

    升级  1%

  • TA的每日心情
    开心
    2014-10-17 07:39
  • 签到天数: 22 天

    [LV.4]偶尔看看III

    自我介绍
    建模苦孩子。
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-10-23 20:05 , Processed in 0.671539 second(s), 76 queries .

    回顶部