function buddhabrot();
% needs tweaking to get the desired image. I superimposed separate output
% from points in c and points outside c to get a mr ebalanced effect.
viewsize=640;
IFSimage = zeros(viewsize,viewsize,3);
xmin=-2;
xmax=1;
ymin=-1.5;
ymax=1.5;
%its=6120;
its=100;
wb = waitbar(0,'Please wait...');
for xx=1:viewsize
    for yy=1:viewsize
        waitbar(xx/viewsize,wb)
        cr=xmin+(xx-1)*(xmax-xmin)/viewsize;
        ci=ymin+(yy-1)*(ymax-ymin)/viewsize; 
        x=0;
        y=0;
        myflag=0;
        for i=1:its
            tmp=x*x-y*y+cr;
            y=2*x*y+ci;
            x=tmp;
            if x*x+y*y>100
                myflag=1; %c is not in M
                break
            end
        end
            x=0;
            y=0;
            for i=1:its
                tmp=x*x-y*y+cr;
                y=2*x*y+ci;
                x=tmp;
                if (x>=xmin)&&(x<=xmax)&&(y>=ymin)&&(y<=ymax)
                    if myflag==1 %if c is not in M add 5 to the value
                        IFSimage(int16((viewsize-1)*(1-(y-ymin)/(ymax-ymin))+1),int16((viewsize-1)*(x-xmin)/(xmax-xmin)+1),1:2) = IFSimage(int16((viewsize-1)*(1-(y-ymin)/(ymax-ymin))+1),int16((viewsize-1)*(x-xmin)/(xmax-xmin)+1),1:2)+5;
                    else %if c is in M add 1 to the value
                        IFSimage(int16((viewsize-1)*(1-(y-ymin)/(ymax-ymin))+1),int16((viewsize-1)*(x-xmin)/(xmax-xmin)+1),3) = IFSimage(int16((viewsize-1)*(1-(y-ymin)/(ymax-ymin))+1),int16((viewsize-1)*(x-xmin)/(xmax-xmin)+1),3)+1;
                    end
                end
                if x*x+y*y>100
                  break
                end

            end   
             end 
    end
close(wb);
save budd2e IFSimage
for i=1:viewsize
    for j=1:viewsize
        for k=1:3
            IFSimage(i,j,k)=(IFSimage(i,j,k)^(1/3));
        end
    end
end
for i=1:3
    IFSimage(:,:,i)=IFSimage(:,:,i)/max(max(IFSimage(:,:,i)));
end
image(IFSimage);
imwrite(IFSimage,'jul.jpg','jpg','Quality',100);