应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 汇编语言 » 求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序
北京石油化工学院2026年研究生招生接收调剂公告
21/1返回列表
查看: 697  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

我要飞

铁虫 (正式写手)

[求助] 求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序

求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-15 21:36:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

★ ★
jjdg(金币+2): 感谢支持 2011-10-16 23:35:49
jjdg: 图真漂亮！牛！ 2011-10-16 23:36:10
CODE:
import numpy as np
import pylab as pl
from matplotlib import cm

def iter_point(c):
    z=c
    tol=1.0e-2
    for i in range(1,1000):
        err=-(z**3-1)/(3*z**2)
        if abs(err)             break
        z=z-(z**3-1)/(3*z**2)
    return i

def draw_newton(cx,cy,d):
    x0,x1,y0,y1=cx-d,cx+d,cy-d,cy+d
    y,x=np.ogrid[y0:y1:200j,x0:x1:200j]
    c=x+y*1j
    newton=np.frompyfunc(iter_point,1,1)(c).astype(np.float)
    pl.imshow(newton,cmap=cm.Blues_r,extent=[x0,x1,y0,y1])
    pl.gca().set_axis_off()
    return

if __name__=='__main__':
    x,y=0.0,0.0
    draw_newton(x,y,0.5)
    pl.show()
赞一下(2人) 回复此楼
2楼2011-10-16 10:04:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 我要飞 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 307分材料专业求调剂 +7 Hll胡 2026-04-05 7/350 2026-04-05 18:47 by 无际的草原
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 化学357分，考研调剂 +10 .Starry. 2026-04-04 11/550 2026-04-05 10:57 by cql1109
[考研] 070300化学学硕311分求调剂 +10 梁富贵险中求 2026-04-04 12/600 2026-04-05 09:37 by guoweigw
[考研] 材料与化工363求推荐 +7 zh096 2026-04-04 7/350 2026-04-05 09:11 by 陌秋26
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 293分求调剂，外语为俄语 +6 加一一九 2026-03-31 6/300 2026-04-04 14:57 by 聪明的大松鼠
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 建环，能源，土木老师路过看一看！！！ +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 化学070300-总分378-求调剂 +5 挪椅子的泡泡糖 2026-04-02 5/250 2026-04-02 22:20 by ZXlzxl0425
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 材料化工340求调剂 +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
信息提示
关闭
请填处理意见
关闭