init:z = #pixelsum = 0.0loop:zold = znwt1 = (@p1 - 1)*z^@p1 + (@c)*sum + @rnwt2 = @p1 * z ^ (@p1 - 1)z = nwt1/nwt2sum = sum + exp(-cabs(z))bailout:|z - zold|>= @bailout