# mandel.ode
#  drawing the mandelbrot set
#
x'=x^2-y^2+cx
y'=2*x*y+cy
cx'=cx
cy'=cy
init x=0,y=0
#
aux amp=x^2+y^2-4
#
#
# Set up Poincare map with SOS and amp=0 as section
# then color code by t (0,10). Integrate 2D range/array with
# -1.5 < cx < .5 and -1 < cy < 1  with 100 as the grid size
# and No to reset storage
 
@ xp=cx,yp=cy,xlo=-1.5,xhi=.5,ylo=-1,yhi=1,lt=-1
@ poimap=section,poivar=amp,poipln=0,poisgn=1,poistop=1
@ maxstor=40000,meth=disc,total=50
done