mandel_plot.m

% compile the mex function
mex mandelbrot.c

close all

% number of iterations
iter=10;
% number of points to use
points=1001;
% square of the complex plane 4; width=4 means [-2,2, -2, 2]
width=4;

% call mex function
W=mandelbrot(iter, points, width);

% choose nice color and plot
colormap vga;
pcolor(W);
shading flat;
axis('square','equal','off');