newton[z_] := z - f[z]/f'[z]
plotplot[r_] :=
ListDensityPlot[Arg@FixedPoint[newton, #, 50] & /@
Table[i + j I, {j, -2.r, 2.r, 02 r/365.011}, {i, -2.r, 2.r, 02 r/365.011}],
ColorFunction -> "Rainbow"]"Rainbow",
DataRange -> {{-r, r}, {-r, r}}]
f[z_] := z^3 - 1; plotplot[2.0]
f[z_] := z^5 - 1; plotplot[2.0]