# original from TANAKA Masaki  
# masaki@isea.is.titech.ac.jp(Tokyo Institute of technology)
compl(a,b)=a*{1,0}+b*{0,1}
mand(z,a,n) = n<=0 || abs(z)>100 ? 1:mand(z*z+a,a,n-1)+1
#
set pm3d depthorder
set palette rgb 7,5,15 
unset colorbox 
unset key
set ticslevel 0
set mapp cart
set view 50,30,1.2
set isosamples 100
splot [-2:1][-1.5:1.5] mand({0,0},compl(x,y),50) w pm3d
pause 3
set term pngcairo transparent size 480,480
set out "mandelbrot.png"
replot