这程序设c的模小于2,迭代次数小于50次
用了个Block函数将z作为局部分量和全局变量区分清楚
通过画密度函数画出c的图像
实际上,这个程序可简化成
M1[x_,y_]:=
Block[{z,k=0},
z=x+y*I
While[(Abs[z]<2.0)&&(k<50),++k
z=z^2+(x+y*I)]
Return[k]
]
M2[px_,py_,p_]:=
Block[{t},
t=DensityPlot[M1[xx,yy],{xx,-1.5,0.5},{yy,-1.2,1.2},p,
Mesh->False,ColorFunction->Hue]
Return[t]
]
Mandelbrot=M2[x,y,{PlotPoints->120,PlotLabel->"Mandelbrot集"}]
你看看行不行
这次我们用循环实现了一个Mandelbrot图形的计算。在《xialulee大战pythonchallenge—...下面是程序运行的结果: 本回答由网友推荐 评论 xiangjuan314 采纳率:94% ...欢迎分享,转载请注明来源:内存溢出
评论列表(0条)