www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著，每章的MATLAB源代码程序 > 第21章/代码21-9.txt

    
iter=27; %最大迭代次数
pixel=400; %迭代步长，即像素点
% 按分辨率划分x轴和y轴
x = linspace(-2.0,1.0,pixel); %x的范围从-2.0到1.0
y = linspace(-1.5,1.5,pixel)'; %y的范围从-1.5到1.5
% 以下定义矩阵C值，包含所有像素点
[Re,Im] = meshgrid(x,y);
C = Re + i * Im;
% 定义零矩阵B，标记发散速度
B = zeros(pixel,pixel);
Cn = B; % C0 = 0+0i，初值始终固定为0
for l = 1:iter %开始迭代
Cn = Cn.*Cn + C;
B = B + (abs(Cn)<2); %当Cn的绝对值大于2时,认为已经逃逸.
end;
% 绘图设置
imagesc(B);          % 画图
colormap(jet);       % 设置色彩索引图
axis equal           % 设置显示比例
axis off             % 不显示坐标轴