No BSD License  

Highlights from
Mandel121

image thumbnail
from Mandel121 by Dr. Jens Koopmann
Cool tool to zoom into the Mandebrot Set. 2/3D images are supported.

[x,y,Zvalues]=dummy2(lowerR,higherR,lowerI,higherI,stepsR,MaxIter,colorMap,viewPoint,x_view,y_view,myLight,axis_value)
function [x,y,Zvalues]=dummy2(lowerR,higherR,lowerI,higherI,stepsR,MaxIter,colorMap,viewPoint,x_view,y_view,myLight,axis_value)


%This is the short routine that calculates the Mandelbrot set
% Parameter Usage
% This function calculates the mandelbrot set within the limits of 
% lowerR higher R lowerI and higherI
% stepsR stepI maxiter define the image resolution, higher settings mean
% better resolution and more CPU time.
% Colormap ranges from 1-18 and defines the look of the result
% viewPoint: 2 or 3 -> 2 or 3D output
% x_view, y_view: 0-360 define the viewing angle of 3D graphs
% myLight: 1-4; defines the lighting effects used for the graph

% Example 
% [x, y, Zvalues ]=dummy2 (-2,2,-2,2,500,700,25,15,3,29,72,3);
% Example increase MaxIter
% [x, y, Zvalues ]=dummy2 (-2,2,-2,2,500,700,100,15,3,29,72,3);
% Example zoom into mandelbrot Set
% [x, y, Zvalues ]=dummy2 (-1.9,0.8,-1.3,1.28,500,700,100,15,3,29,72,3);

handles.col=colorMap;
handles.view=viewPoint;
stepsI=stepsR;


     %compute other constants  
    Rwidth=higherR-lowerR;
    Iwidth=higherI-lowerI;
    slR=Rwidth/(stepsR-1);
    slI=Iwidth/(stepsI-1);
    [x,y]=meshgrid((0:stepsR-1)*slR+lowerR,(0:stepsI-1)*slI+lowerI);
    Zvalues=ones(size(x));
    initZ=zeros(size(x));
    c=(x+i*y);
    obergrenze=stepsR*stepsI;
    z=initZ;
    ddx=1/MaxIter;
    h = waitbar(0,'Please wait...');
    h_z=1:(obergrenze);
    for counter=1:MaxIter 
      ddx=ddx+ddx;
      waitbar(ddx)  
      z(h_z)=z(h_z).^2+c(h_z);
      h_z= h_z(logical(abs(z(h_z))<2));
      Zvalues(h_z)=Zvalues(h_z)+1;
    end
    close(h);
    color_change(handles.col,handles.view,x,y,Zvalues,x_view,y_view,MaxIter,myLight,axis_value);
 
end

Contact us at files@mathworks.com