F = @(z) z^3-1;
Fprime = @(z) 2*z^2;
z0 = 0;
d = 1;
n = 2^11;
m = 10^8;
maxiter = 2^7;

newton_fractal(F,Fprime,z0,d,n,maxiter,m)