|
unset key
set view map
set bmargin at screen 0.1
set sample 200
set isosamples 200, 200
set size ratio 1 1,1
set cbtics ("0" -pi, "2π" pi)
set title "Complex error function cerf( x + iy )"
set xrange [ -3 : 3 ]
set yrange [ -3 : 3 ]
set tics scale 0
set cblabel "Phase Angle"
set cblabel offset character -2, 0, 0 font "" textcolor lt -1 rotate by -270
set cbrange [ -pi : pi ] noreverse nowriteback
set contour
set cntrparam levels discrete 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50, 100, 200, 500
set cntrlabel onecolor
set palette positive nops_allcF maxcolors 0 gamma 1.5 color model HSV
set palette defined ( 0 0 1 1, 1 1 1 1 )
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))
ip(x,y) = imag(f(x,y))
color(x,y) = hsv2rgb( Hue( rp(x,y), ip(x,y) ), abs(f(x,y)), 1. )
f(x,y) = cerf(x+y*{0,1})
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
Click here for minimal script to generate this plot
|