±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 698  |  »Ø¸´: 1

ÎÒÒª·É

Ìú³æ (ÕýʽдÊÖ)

[ÇóÖú] ÇóÖúmatlabÀïÃæ»­³öx^3-1=0µÄNewtonµü´úµÄÎüÒýÓòµÄ³ÌÐò

ÇóÖúmatlabÀïÃæ»­³öx^3-1=0µÄNewtonµü´úµÄÎüÒýÓòµÄ³ÌÐò
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢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Â¥2011-10-16 10:04:54
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ÎÒÒª·É µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Çóµ÷¼Á +8 xzghyuj 2026-04-04 8/400 2026-04-06 06:47 by houyaoxu
[¿¼ÑÐ] 0855Çóµ÷¼Á²ÄÁÏ +11 ºìÌÒׯׯ 2026-04-04 11/550 2026-04-05 22:13 by ×íÎÌwl
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϺ½£¬ÊýһӢһѧ˶317Çóµ÷¼Á£¡£¡ +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by ²«»÷518
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜ ³õÊÔ316·Ö Çóµ÷¼Á +6 ²ÐÐÇ·÷Êï 2026-03-31 6/300 2026-04-05 12:15 by rainbow11
[¿¼ÑÐ] 0835ѧ˶299Çóµ÷¼Á 08´óÀà¿É½ÓÊÜ +5 useryy 2026-04-03 5/250 2026-04-04 20:07 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 298Çóµ÷¼Á +5 zzz£¬£¬r 2026-04-02 8/400 2026-04-04 19:55 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 359Çóµ÷¼Á +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÄÏʦ·¶361·Ö£¬»¯Ñ§Çóµ÷¼Á +7 Nicole88888 2026-04-01 7/350 2026-04-04 18:28 by macy2011
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄÏÅ©090401£¬268£¬Çóµ÷¼Á +5 һľÄñÈ» 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[¿¼ÑÐ] 325Çóµ÷¼Á +4 ´º·ç²»½èÒâ 2026-04-04 4/200 2026-04-04 14:46 by ÏæÅ©´¢ÄܲÄÁÏ
[¿¼ÑÐ] ÇóÉúÎïѧרҵµ÷¼Á-332·Ö +5 ÔÆ¶äåÞÍäÖ¸ÄÏ 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
[¿¼ÑÐ] Çóµ÷¼Á +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[¿¼ÑÐ] 320Çóµ÷¼Á +3 ũҵ¹¤³ÌÓëÐÅÏ¢¼ 2026-04-03 3/150 2026-04-03 11:40 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 081200-11408-276ѧ˶Çóµ÷¼Á +6 ´Þwj 2026-04-02 6/300 2026-04-03 10:19 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Ò»Ö¾Ô¸Éî´ó085601²ÄÁϹ¤³Ìרҵ£¨×¨Ë¶£©300·Ö¿ÉÒÔµ÷¼ÁÈ¥ÄÄ +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[¿¼ÑÐ] 293Çóµ÷¼Á +4 çæçæÀÖ 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[¿¼ÑÐ] 318Çóµ÷¼Á£¬¼ÆËã²ÄÁÏ·½Ïò +10 Îüß÷Óк¦óÏÃü 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +11 Ò»ÑùYWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[¿¼ÑÐ] 262Çóµ÷¼Á +9 ÀøÖ¾Ò»¶¨·¢ÎÄÕ 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[¿¼ÑÐ] 301Çóµ÷¼Á +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by ÎҵĴ¬Îҵĺ£
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û