QQ登录

只需要一步,快速开始

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

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

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

712

主题

16

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2026-7-5 10:02
  • 签到天数: 3998 天

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    跳转到指定楼层
    1#
    发表于 2014-8-31 15:19 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    我想用mathematica编写一个求根函数,用牛顿切线法,我编的函数如下:* C2 z; T/ J% e6 a% d
    Newton[f_,{x_,x0_},e_]:=- Y. `( p7 s$ k/ C' g
    Module[{df,x1,y0,y1,k},; @0 v8 _# r. ^/ u
    df=D[f,x];
    8 n& `" k% E/ R2 m) Ay0=f/.x->x0;y1=df/.x->x0;
    * h/ R9 Y* l/ k0 ?# D  j- R7 rx1=x0-y0/y1;
    $ _- m8 Z0 t  b- Ck=1;
    - ^+ ?( g5 W% D+ W$ YWhile[Abs[x1-x0]>e,! y8 a, u: G$ M, D6 H5 Y
    x0=x1;
      a4 l5 `/ P* m- b8 W& I1 x$ h9 I. O) X: Gy0=f/.x->x0;y1=df/.x->x0;
    $ k1 H. ]' b6 v1 |9 I+ s8 px1=x0-y0/y1;; q! @7 x) P# ~1 D
    k=k+1;' J9 X' F( e. W  y
    ]2 N. F& F! |. V/ Y& ^/ n, y
    Print["迭代次数为:",k];
    / Q2 I  q' X5 J" \( |' uPrint["近似根为:",x1];
    ! o3 Q" U! x5 ^! q9 c], _) p" O) ^: ~8 |2 T

    : N! w# ^0 S2 X  ^% N. z但是当我调用它时,会出现这样的提示,如下:+ P7 L1 l# {/ G
    , n1 b1 w  K- s! R+ O0 ~' }' E
    这是为什么啊?求高手解答,应该怎么修改啊
    & f6 p! t2 z5 f7 q
    7 m' A+ e8 B* {; z, o1 |4 L& @
    / W% q) N' o: o

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

    1.PNG

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

    712

    主题

    16

    听众

    1万

    积分

    升级  0%

  • TA的每日心情
    开心
    2026-7-5 10:02
  • 签到天数: 3998 天

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    回复

    使用道具 举报

    madio        

    3万

    主题

    1312

    听众

    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];
      ) T' F8 J. S5 P
    2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;# v# u0 F# A# I' }
    3.   x1 = x0 - y0/y1;
      ) H5 i5 b: W- T1 [6 k0 |4 A
    4.   k = 1;- Z5 Y  x+ h# |; W1 G
    5.   While[Abs[x1 - x0] > e, x0 = x1;
      6 n! ?1 U9 }, s: U9 g
    6.     y0 = f /. x -> x0; y1 = df /. x -> x0;' n' [( k- d: v6 M
    7.     x1 = x0 - y0/y1;\" [$ b9 x8 t& F$ K1 b. _2 g
    8.     k = k + 1;] Print["迭代次数为:", k];
      . u8 G0 S8 g1 ~\" G( ]! M- u
    9.   Print["近似根为:", x1];]
    复制代码
    可选参数应该放在一个表里面
    数学建模社会化
    回复

    使用道具 举报

    lxy444        

    712

    主题

    16

    听众

    1万

    积分

    升级  0%

  • TA的每日心情
    开心
    2026-7-5 10:02
  • 签到天数: 3998 天

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    madio 发表于 2014-8-31 21:59
    ' a, M; q5 U( t5 q9 }3 Z8 Y可选参数应该放在一个表里面

    ( M6 G5 Z0 _. P" x8 }哦,明白了,十分感谢!
    回复

    使用道具 举报

    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, 2026-8-7 16:01 , Processed in 0.522005 second(s), 76 queries .

    回顶部