# Author: Ika, 2015-10-17
# Anyone can copy and revise this code without any copyright concern.

# This program generates the imaginary part of g*exp(bz^2).

set term svg size 800,600 enhanced font 'Times, 12'
set output 'gexp_im.svg'

set multiplot

set xlabel "x"
set ylabel "iy"
set zlabel "Gexp"
set xrange [-0.5:0.5]
set yrange [-0.5:0.5]
unset colorbox
set view 50,20

set border lw 0.5 
set key on left top

set hidden3d front
set ticslevel 0
set isosamples 30,30

# Draw the surfaces of f(x,y) = g*exp(b*(x^2-y^2))sin(2bxy) and y=1.
splot	exp(0.5*(x*x-y*y))*sin(2*0.5*x*y) lc rgb "blue" t "g=1, b=0.5",\
		exp(1*(x*x-y*y))*sin(2*1*x*y) lc rgb "green" t "g=1, b=1",\
		exp(-2*0.5*x*y)*sin(0.5*(x*x-y*y)) lc rgb "red" t "g=1, b=0.5i",\
		exp(-2*1*x*y)*sin(1*(x*x-y*y)) lc rgb "black" t "g=1, b=i",\
		exp(x*x-y*y-2*x*y)*sin(2*x*y+x*x-y*y) lc rgb "purple" t "g=1, b=1+i"

unset multiplot
set term x11
set output