This is the formula I've been using lately and have really enjoyed it. I'd be curious to see what others can do with it

All I did was add funtions to all terms of dmj.ufm(Nova (Mandelbrot)) and clear out some of the stuff I didn't understand

JABMOD {
;        a modification of Nova (Mandelbrot) from dmj.ufm
;
init:
  complex zsquared = (0,0)
  complex zcubed = (0,0)
  complex zold = (0,0)
  
  z = @start
  
loop:
  IF (@power == (3,0)); special optimized routine for power 3
    zsquared = sqr(z)
    zcubed = zsquared * z
    zold = z
    z = z - @relax * @func(zcubed-1) / @func2(3*zsquared) + @func3(#pixel)
  ELSE
    zold = z
    z = z - @relax * @func(z^@power-1) / @func2(@power * z^(@power-1)) + @func3(#pixel)
  ENDIF
  
bailout:
  |z - zold| > @bailout
  
default:
  title = "Nova (Mandelbrot) JABMOD"
  maxiter = 10000
  periodicity = 0
  center = (-0.5,0)
  magn = 1.5
  
  param start
    caption = "Start Value"
    default = (2,0)
    hint = "Starting value for each point.  You can use this to
            'perturb' the fractal."
  endparam
  param power
    caption = "Exponent"
    default = (5,0)
    hint = "Overall exponent for the equation. The value (3,0) gives
            the classic Nova Mandelbrot-type fractal."
  endparam
  param bailout
    caption = "Bailout"
    default = 0.000001
$IFDEF VER40
    exponential = true
$ENDIF
    hint = "Bailout value; smaller values will cause more
            iterations to be done for each point."
  endparam
  param relax
    caption = "Relaxation"
    default = (1,0)
    hint = "This can be used to slow down the convergence of
            the formula."
  endparam


}

If you use it I'd love to see the result! Please bring it to my attention!