1 %
   2 %	Filename: mandelplot.m
   3 %
   4 function mandelplot(run,anz,filename)
   5 my_res = 512;
   6 my_x = 0.87485;
   7 my_y = 0.75145;
   8 my_log_r0 = 1.5;
   9 my_log_r1 = -12;
  10 my_log_r = my_log_r0 + my_log_r1 - my_log_r0*(run-1)/anz;
  11 my_r = exp(my_log_r);
  12 x1 = my_x - my_r;
  13 x2 = my_x + my_r;
  14 y1 = my_y - my_r;
  15 y2 = my_y + my_r;
  16 x=[x1:(x2-x1)/my_res:x2];
  17 y=[y1:(y2-y1)/my_res:y2];
  18 m=zeros(my_res);
  19 m=mandelbrot(my_res,x1,x2,y1,y2);
  20 
  21 %clf;surf(x,y,m)
  22 %axis tight;
  23 %shading interp; map=colormap;
  24 %map(1,3)=0;
  25 %colormap(map);
  26 %view(2)
  27 
  28 imwrite(m,filename);
  29 end