版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (2879)
>虫友互识 (340)
>导师招生 (149)
>考博 (143)
>文献求助 (116)
>论文投稿 (74)
>基金申请 (72)
>硕博家园 (64)
>找工作 (51)
>博后之家 (50)
>仿真模拟 (49)
>休闲灌水 (40)
>公派出国 (33)
>教师之家 (32)
>考研 (31)
>招聘信息布告栏 (28)

返回列表

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

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

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

回复此楼

» 猜你喜欢

河北省自然科学基金已经有8人回复
西安交大新媒学院副院长用撤稿论文结题已经有5人回复
论文撤稿了已经有5人回复
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！已经有5人回复
26/27申博自荐已经有9人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有7人回复
揭秘青基评审内幕：几个A才能顺利中标已经有4人回复
青B发送上会通知了吗已经有7人回复
博士申请已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！已经有11人回复
求用牛顿迭代法解非线性方程组的matlab程序已经有4人回复
请教matlab如何增加迭代次数已经有11人回复
求用牛顿迭代法求非线性方程组已经有4人回复
matlab数值计算时迭代项分母为零的时候怎么办？已经有15人回复
matlab菜鸟问3个问题~~~~~~~~~~~ 已经有2人回复
求助matlab里面画出x^3-1=0的Newton迭代的吸引域的程序已经有5人回复
求助matlab画出螺旋线的程序已经有6人回复
【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】已经有15人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼 2011-10-15 21:36:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
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的回帖

相关版块跳转我要订阅楼主我要飞的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 河北省自然科学基金 +5	Peterchao 2026-05-18	8/400	2026-05-24 11:58 by 晓晓爱翠翠
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3	bjvtcliu 2026-05-24	5/250	2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3	bjvtcliu 2026-05-24	5/250	2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！ +4	zju2000 2026-05-21	5/250	2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4	ZXW0202 2026-05-22	9/450	2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	pmo95bazuy 2026-05-23	7/350	2026-05-24 06:35 by fpo5ljpv91
[基金申请] 揭秘青基评审内幕：几个A才能顺利中标 +3	国自然国社科中 2026-05-23	4/200	2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5	chemBioBro 2026-05-22	7/350	2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3	焦晓明 2026-05-21	3/150	2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助，期刊 +4	希冀，有书读 2026-05-20	8/400	2026-05-22 10:16 by 希冀，有书读
[文学芳草园] 献血感触 +7	呀呀好傻 2026-05-19	13/650	2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页，违规吗？会被低分嘛？ +14	1234567wang 2026-05-17	16/800	2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4	无名者登山 2026-05-20	5/250	2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7	无名者登山 2026-05-18	11/550	2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9	fummck 2026-05-20	10/500	2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15	popular289 2026-05-18	26/1300	2026-05-21 10:35 by 西葫芦炒鸡蛋
[有机交流] 反应很差，大量原料没有反应 5+3	Mr.Zot 2026-05-19	8/400	2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3	铁达火车 2026-05-18	5/250	2026-05-20 09:33 by tfang
[考博] 博士申请 +5	星…… 2026-05-18	6/300	2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3	Love_MH 2026-05-17	3/150	2026-05-18 02:22 by 竹林孤影