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!