#-------------------------------------------------------------
#
#  Sample code
#              The Julia set for map     z->z**2 + C
#              with                      C = (0.32,0.043)
#
#------------------------------------------------------------

function list {
    creal = 0.32;
    cimag = 0.043;
    J(x,y,z)=((z>100)?(0.1):((x*x+y*y>10000.0)?(0.001*z): \
                  (J(x*x-y*y+creal,2.0*x*y+cimag,z+1))) );
    };
julia = surface{[J(x,y,0)][x=-1.0:1.0][y=-1.2:1.2]\
                        [sample=100:100]};

#-----------------------------------------------------------
plot2d julia;
#-----------------------------------------------------------

