incr = 0.5;
xMin = -2;
xMax = 2;
yMin = -2;
yMax = 2;

x = -2:incr:2;
y = -2:incr:2;
[X,Y] = meshgrid(x,y);
counter = 0;
maxIterations = 25;

z = [];
a = 0;
for(i=xMin:incr:xMax)
   a = a + 1;
   b = 0;
    for(j=yMin:incr:yMax)
        b = b + 1;
        z(a,b) = complex(i,j);

        while(counter < maxIterations)
            z(a,b) = z(a,b) - f(z(a,b))/fderiv(z(a,b));
            counter = counter + 1;
        end
        counter = 0;
    end
end