Formula:
init:
z = #pixel
loop:
zold = z
nwt1 = (@p1 - 1) * z^@p1 + (@p2 - 1)*z^@p2 + @r
nwt2 = @p1 * z ^ (@p1 - 1)+ @p2 * z ^ (@p2 - 1)
zn = nwt1/nwt2
nwtm1 = zn^3 - 1
nwtm2 = 3*zn^2
z = z -(nwtm1/nwtm2)
bailout:
|z - zold|>= @bailout
The last one used @p1 = 3,@p2 = 0 and Two Layer(1 - Exponent Smoothing,2 - Smothing Iteration)