# set terminal pngcairo  background "#ffffff" enhanced font "times" fontscale 1.0 size 640, 480 
# set output 'complex_trig.3.png'
set bar 1.000000 front
set style circle radius graph 0.02, first 0.00000, 0.00000 
set style ellipse size graph 0.05, 0.03, first 0.00000 angle 0 units xy
unset key
set style textbox transparent margins  1.0,  1.0 border
unset logscale
set view map scale 1
set isosamples 100, 100
set size ratio 1 1,1
set style data lines
set xtics  norangelimit 
set xtics   ("-π/2" -1.57080, "-π/4" -0.785398, "0" 0.00000, "π/4" 0.785398, "π/2" 1.57080)
set ytics  norangelimit 
set ytics   ("-π/2" -1.57080, "-π/4" -0.785398, "0" 0.00000, "π/4" 0.785398, "π/2" 1.57080)
set cbtics  norangelimit 
set cbtics   ("0" -3.14159, "2π" 3.14159)
unset paxis 1 tics
unset paxis 2 tics
unset paxis 3 tics
unset paxis 4 tics
unset paxis 5 tics
unset paxis 6 tics
unset paxis 7 tics
set title "acos( x + iy )" 
set xrange [ -1.57080 : 1.57080 ] noreverse nowriteback
set yrange [ -1.57080 : 1.57080 ] noreverse nowriteback
set cblabel "Phase Angle" 
set cblabel  offset character -2, 0, 0 font "" textcolor lt -1 rotate by -270
set cbrange [ -3.14159 : 3.14159 ] noreverse nowriteback
set paxis 1 range [ * : * ] noreverse nowriteback
set paxis 2 range [ * : * ] noreverse nowriteback
set paxis 3 range [ * : * ] noreverse nowriteback
set paxis 4 range [ * : * ] noreverse nowriteback
set paxis 5 range [ * : * ] noreverse nowriteback
set paxis 6 range [ * : * ] noreverse nowriteback
set paxis 7 range [ * : * ] noreverse nowriteback
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, 0 size screen 0.05, 0.6, 0 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) = acos(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. )
save_encoding = "utf8"
GPFUN_Hue = "Hue(x,y) = (pi + atan2(-y,-x)) / (2*pi)"
GPFUN_phase = "phase(x,y) = hsv2rgb( Hue(x,y), sqrt(x**2+y**2), 1. )"
x = 0.0
GPFUN_rp = "rp(x,y) = real(f(x,y))"
GPFUN_ip = "ip(x,y) = imag(f(x,y))"
GPFUN_color = "color(x,y) = hsv2rgb( Hue( rp(x,y), ip(x,y) ), abs(f(x,y)), 1. )"
GPFUN_f = "f(x,y) = acos(x + y*{0,1})"
## Last datafile plotted: "++"
splot '++' using 1:2:(color($1,$2)) with pm3d lc rgb variable