# set terminal pngcairo  transparent enhanced font "arial,10" fontscale 1.0 size 600, 400 
# set output 'BesselK.4.png'
set border -1 front lt black linewidth 1.000 dashtype solid
unset key
set view 42, 33, 1, 1
set samples 300, 300
set isosamples 300, 300
set hidden3d back offset 1 trianglepattern 3 undefined 1 altdiagonal bentover
set style data lines
set xyplane relative 0
set xtics  norangelimit 1 
set ytics  norangelimit 1 
set ztics  norangelimit 1 
set cbtics  norangelimit 
set cbtics   ("-π" -3.14159, "0" 0.00000, "π" 3.14159)
set rtics  norangelimit 1 
set title "{/Times=16 complex Bessel function} {/Times:Italic=16 |K_{1/2} (x + iy)|}" 
set xlabel "x" 
set xrange [ -1.57080 : 1.57080 ] noreverse nowriteback
set x2range [ * : * ] noreverse writeback
set ylabel "iy" 
set yrange [ -1.57080 : 1.57080 ] noreverse nowriteback
set y2range [ * : * ] noreverse writeback
set zrange [ -2.00000 : 8.00000 ] noreverse writeback
set cbrange [ -3.14159 : 3.14159 ] noreverse writeback
set rrange [ * : * ] noreverse writeback
set pm3d lighting primary 0.5 specular 0.2 spec2 0.6
set pm3d spotlight rgb 0xff0044 rot_x 50 rot_z 90 Phong 4.0
set palette cubehelix start 0.5 cycles -1.5 saturation 1
set colorbox vertical origin screen 0.9, 0.2 size screen 0.05, 0.6 front  noinvert bdefault
NO_ANIMATION = 1
save_encoding = "utf8"
## Last datafile plotted: "++"
splot '++' using 1:2:( abs(BesselK( 0.5, x+y*I )) ):( arg(BesselK( 0.5, x+y*I )) ) with pm3d