# set terminal pngcairo  transparent enhanced font "arial,10" fontscale 1.0 size 600, 400 
# set output 'cerf.3.png'
unset key
set view map scale 1
set samples 200, 200
set isosamples 200, 200
set contour base
set cntrlabel onecolor format '%8.3g' font '' start 5 interval 20
set cntrparam levels 12
set cntrparam levels discrete 0.1,0.2 ,0.5 ,1 ,2 ,5 ,10 ,20 ,50 ,100 ,200 ,500 
set size ratio 1 1,1
set xtics border in scale 0,0 mirror norotate  autojustify
set ytics border in scale 0,0 mirror norotate  autojustify
set ztics border in scale 0,0 nomirror norotate  autojustify
set cbtics border in scale 0,0 mirror norotate  autojustify
set cbtics  norangelimit 
set cbtics   ("0" -3.14159, "2π" 3.14159)
set rtics border in scale 0,0 nomirror norotate  autojustify
set title "Complex error function cerf( x + iy )" 
set trange [ * : * ] noreverse nowriteback
set urange [ -3.00000 : 3.00000 ] noreverse nowriteback
set vrange [ -3.00000 : 3.00000 ] noreverse nowriteback
set xrange [ * : * ] noreverse writeback
set x2range [ * : * ] noreverse writeback
set yrange [ * : * ] noreverse writeback
set y2range [ * : * ] noreverse writeback
set zrange [ * : * ] noreverse writeback
set cblabel "Phase Angle" 
set cblabel  offset character -2, 0, 0 font "" textcolor lt -1 rotate
set cbrange [ -3.14159 : 3.14159 ] noreverse nowriteback
set rrange [ * : * ] noreverse writeback
set bmargin at screen 0.1
set tmargin at screen 0.9
set palette positive nops_allcF maxcolors 0 gamma 1.5 color model HSV 
set palette defined ( 0 0 1 1, 1 1 1 1 )
set colorbox vertical origin screen 0.9, 0.2 size screen 0.05, 0.6 front  noinvert bdefault
Hue(x,y) = (pi + atan2(-y,-x)) / (2*pi)
phase(x,y) = hsv2rgb( Hue(x,y), sqrt(x**2+y**2), 1. )
rp(x,y) = real(f(x,y))
f(x,y) = cerf(x+y*{0,1})
ip(x,y) = imag(f(x,y))
color(x,y) = hsv2rgb( Hue( rp(x,y), ip(x,y) ), abs(f(x,y)), 1. )
NO_ANIMATION = 1
save_encoding = "utf8"
fwhm1 = 3.60287466889735
fwhm2 = 3.63018986329991
fwhm3 = 3.60113567720316
fwhm4 = 3.6
## Last datafile plotted: "++"
splot '++' using 1:2:(color($1,$2)) with pm3d lc rgb variable nocontour,       '++' using 1:2:(abs(cerf($1+$2*{0,1}))) with lines nosurf lt -1