#-------------------------------------------------------------
#
#  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))) );
    };
graph list {
       j_half = surface{[J(x,y,0)][x=-1.0:1.0][y=0.0:1.2]\
                        [sample=200:100]};
    };
Julia_set = object{
              j_half;
              pushmatrix
                scale -1.0,-1.0,1.0;
                 j_half;
              popmatrix;
    };
#-----------------------------------------------------------
center;
plot2d Julia_set;
#-----------------------------------------------------------

