M集是:
- Z(0) = 0 + 0i
- Z(n+1) = Z(n)^2 + C n=1,2,3,… 当|Z|<2时C的集合
M.m
function J = M(a,b,c,d,J,ite)
% Explanation:
% a : Z的实部虚部取值下界
% b : Z的实部虚部取值上界
% c : C的虚部
% d : C的实部
% J : Z的初始取值
% ite: 迭代次数N
N = a:0.001:b;
[Re,Im] = meshgrid(N);
Z = Re+Im*1i;
C = c*sqrt(-1)-d;
for n = 1:ite
t = abs(Z)<=2;
J = J+t;
Z = Z.^2+C;
end
main.m
clear all,clc;
J = M(-2,2,0.156,0.9,0,50);
imagesc(J);
colormap('jet');
axis equal
axis off