QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2834|回复: 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编写一个求根函数,用牛顿切线法,我编的函数如下:
    2 U9 j( ^! m$ y6 P0 tNewton[f_,{x_,x0_},e_]:=
    ( ~" i4 e2 N% c# Y" {Module[{df,x1,y0,y1,k},
    : n* d& E8 f: P6 odf=D[f,x];
    + n% k* z) S3 P9 u7 Qy0=f/.x->x0;y1=df/.x->x0;
    ! ^0 ?) s% O* T' c' a. g+ rx1=x0-y0/y1;
    5 p7 @& b( b- V& {k=1;
    ( S) B/ }2 n) d* s/ H5 y7 K8 C3 B5 VWhile[Abs[x1-x0]>e,2 F2 O& u, l3 c9 \& \1 \
    x0=x1;0 _8 {0 X, Q$ }( x1 B, l
    y0=f/.x->x0;y1=df/.x->x0;
    6 T! J4 G1 k8 m, C+ q& d- M0 G; _x1=x0-y0/y1;% g( E# M6 ]5 X! f
    k=k+1;! Y4 \2 S2 [! u1 r$ p, k5 |
    ]
    5 A& p* }0 j2 G( V8 ?& O( }0 PPrint["迭代次数为:",k];+ {3 G; k6 v! o: `4 M
    Print["近似根为:",x1];! F' H! n. _# L- _) ?0 O; \
    ]
    ! E) H7 K+ P: ~! {/ S, D* Z0 k" `+ L9 A! h" H! F; r
    但是当我调用它时,会出现这样的提示,如下:( M  x2 Z" ?. b7 K

    2 Y! C: u; A/ X) w* h9 h: d这是为什么啊?求高手解答,应该怎么修改啊) x$ ?) X& [$ }% b
    : ]# b+ x+ N' s: v1 {

    4 u/ O! C0 ]0 L) ^

    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];$ y: B1 k% z% u) O
    2.   y0 = f /. {x -> x0; y1 = df} /. x -> x0;
      . I6 \0 M1 B0 u( h7 k4 g
    3.   x1 = x0 - y0/y1;
      ( ]9 h  B) ]/ v# e. ]7 s
    4.   k = 1;7 [% m. z0 }/ J0 J! v
    5.   While[Abs[x1 - x0] > e, x0 = x1;% `7 n8 w4 `, T  W9 G
    6.     y0 = f /. x -> x0; y1 = df /. x -> x0;
      ! M( P5 c4 y% |
    7.     x1 = x0 - y0/y1;; N* O6 ]: A$ F1 t6 g\" t2 |( M
    8.     k = k + 1;] Print["迭代次数为:", k];
      % P8 P* B- }2 u3 W( L: ?* B
    9.   Print["近似根为:", x1];]
    复制代码
    可选参数应该放在一个表里面
    数学建模社会化
    回复

    使用道具 举报

    lxy444        

    712

    主题

    16

    听众

    1万

    积分

    升级  0%

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

    [LV.Master]伴坛终老

    邮箱绑定达人

    群组2013年数学建模国赛备

    madio 发表于 2014-8-31 21:59
    5 G! H# U5 b' n$ b( h! d9 |% e0 g可选参数应该放在一个表里面

    & A4 P$ P4 H0 p, U3 N哦,明白了,十分感谢!
    回复

    使用道具 举报

    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:16 , Processed in 0.580267 second(s), 75 queries .

    回顶部