conj(z) = z - {0,2}*imag(z)
Fz(z,r,n) = 0 < n ? r%2 == 0 ? Fz(zb*conj(z), r/2, n-1): \
                   Fz(zc*(z - 1) + 1, r/2, n-1) : z
unset border
unset colorbox
unset tics
unset key
set pm3d nohidden3d
set pal rgb 15,7,5
N = 14
Range = 2**N
zb = {0.4616, 0.4616}
zc = {0.6200, -0.197}
set samples Range
XL = 1.1
set xrange [0:XL]
set view 0,0,1.6,
set label sprintf("b = %.4f %+.4fi, c = %.4f %+.4fi", \
                  real(zb), imag(zb), real(zc), imag(zc)) \
                  at 0.2,-0.05 font ",12" 
set term pngc enh font "Times,10" size 480,310
set out "fern_complex.png"
Zo = {0,0}
splot '+' u (z = Fz(Zo,int($1*Range/XL),N)):(imag(z)):\
            (int($1*Range)%64) w d lt palette