function mandel(z0,z1,maxsteps,hres,fla) %построение множества Мандельброта %z0 - координата левой нижней точки %z1 - координата правой верхней точки %maxsteps - число шагов для вычисления %hres - дискретность по оси x x0 = real(z0); y0 = imag(z0); x1 = real(z1); y1 = imag(z1); vres = ceil( abs( hres * (y1-y0)/(x1-x0)) ); dx = (x1-x0)/hres; dy = (y1-y0)/vres; [X,Y] = MESHGRID(x0:dx:x1, y0:dy:y1); c = X + Y*i; z = c; pict = 0*z; clear X Y for k=1:maxsteps bigs = find(abs(z)>=2); pict(bigs) = k*ones(size(bigs)); z(bigs) = zeros(size(bigs)); c(bigs) = zeros(size(bigs)); z = z.^2 + c; end; if (fla==1), figure;end; image(pict) title('Множество Мандельброта') colormap(hsv(maxsteps)); axis('image') axis('xy') axis('off')