# set terminal pngcairo  transparent enhanced font "arial,10" fontscale 1.0 size 600, 400 
# set output 'expint.2.png'
set border -1 front lt black linewidth 1.000 dashtype solid
set grid nopolar
set grid xtics nomxtics ytics nomytics ztics nomztics nortics nomrtics \
 nox2tics nomx2tics noy2tics nomy2tics nocbtics nomcbtics
set grid vertical layerdefault   lt 0 linecolor 0 linewidth 0.500,  lt 0 linecolor 0 linewidth 0.500
unset key
set view 35, 43, 0.84, 1.4
set samples 101, 101
set isosamples 101, 101
set xyplane relative 0
set xtics border in scale 1,0.5 mirror norotate  offset character 0, -0.5, 0 autojustify
set ytics border in scale 1,0.5 mirror norotate  offset character 0, -0.5, 0 autojustify
set ztics  norangelimit 2 
set cbtics  norangelimit 
set cbtics   ("-π" -3.14159, "+π" 3.14159)
set title "Complex Exponential Integral\n\n {/:Italic {E_n(z) = }{/*2 ∫@_{/*0.5 1}^{/*.5 ∞}}{t^ne^{-tz}}}       for {/:Italic n=2}" 
set title  offset character -25, 0, 0 font "" textcolor lt -1 norotate
set urange [ -2.00000 : 2.00000 ] noreverse nowriteback
set vrange [ -2.00000 : 2.00000 ] noreverse nowriteback
set xlabel "Real(z)" 
set xrange [ -2.00000 : 2.00000 ] noreverse nowriteback
set x2range [ * : * ] noreverse writeback
set ylabel "Imag(z)" 
set yrange [ -2.00000 : 2.00000 ] noreverse nowriteback
set y2range [ * : * ] noreverse writeback
set zlabel "Abs(E_2(z))" 
set zlabel  font "" textcolor lt -1 rotate
set zrange [ 0.00000 : 8.00000 ] noreverse nowriteback
set cblabel "Phase Angle" 
set cblabel  offset character -1, 0, 0 font "" textcolor lt -1 rotate
set cbrange [ -3.14159 : 3.14159 ] noreverse writeback
set rrange [ * : * ] noreverse writeback
set pm3d interpolate 1,1 flush begin noftriangles noborder corners2color c1
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
NO_ANIMATION = 1
save_encoding = "utf8"
## Last datafile plotted: "++"
splot '++' using 1:2:(abs( expint(2,x+y*I))):(arg( expint(2,x+y*I))) with pm3d