\(\(ni = 100; \  (*\ the\ number\ of\ iteration\ *) 

    it[z_, c_] := 
      Module[{r, j, x}, 
        \((r = \((1 + Sqrt[1 + 4  Abs[c]])\)/2; x = z; j = 0; 
          While[Abs[x] < r\  && \ j < ni + 1, x = x*x + c; \(j++\)]; 
          Return[j])\)]; 
    
 (*\ function\ giving\ the\ escape\ time\ of\ the\ point\ *) \)\)

\(c0 = 1; 

    DensityPlot[it[x + I\ y, c0], {x, \(-1\), 1}, {y, \(-1.5\), 1.5}, 
      PlotPoints -> 200, Mesh -> False, 
      AspectRatio -> Automatic]\)

\(c0 = \(- .13\) +  .75  I; 

    DensityPlot[it[x + I\ y, c0], {x, \(-1.5\), 1.5}, {y, \(-1.5\), 1.5}, 
      PlotPoints -> 200, Mesh -> False, 
      AspectRatio -> Automatic]\)

\(c0 = \(- .13\) +  .75  I; 

    Plot3D[it[x + I\ y, c0], {x, \(-1.5\), 1.5}, {y, \(-1.5\), 1.5}, 
      PlotPoints -> 100, Mesh -> False, 
      AspectRatio -> Automatic]\)