Limulus {
;
; adapted Newton fractal
; with z -> sin(z) transform
init:
  z = pixel
loop:
  z1 = sin(z)-p1
  z2 = z1^p2-1
  z3 = p2*(z1^(p2-1))
  z=z-(z2/z3)
bailout:
  |z2| >= 0.00001
default:
  title = "Limulus"
  helpfile = "Uf3.chm"
  helptopic = "Html\formulas\standard\newton.html"
  maxiter = 150
  param p1
    caption = "relaxation factor"
    default = (0.1,0)
    endparam
  param p2
    caption = "Exponent"
    default = (3,0)
    hint = "Specifies the exponent of the equation that is solved by \
            Newton's method. Use real numbers (set the imaginary component \
            to zero) to obtain classic Newton fractals."
  endparam
}