load "gnuplot_beta_defs"

w(x) = omega_n1_q2(x)
f(x,y) = (w(y)  + w(y-3) + w(y+3))*( (w(x) + w(x-3) + w(x+3))*1.2 + (w(x-1) + w(x-4) + w(x+2))*0.3 - (w(x-2) + w(x-5) + w(x+1))*0.2) \
       + (w(y-1)+ w(y-4) + w(y+2))*(-(w(x) + w(x-3) + w(x+3))*0.1 + (w(x-1) + w(x-4) + w(x+2))*0.5 + (w(x-2) + w(x-5) + w(x+1))*0.1) \
	   + (w(y-2)+ w(y-5) + w(y+1))*( (w(x) + w(x-3) + w(x+3))*0.5 - (w(x-1) + w(x-4) + w(x+2))*0.3 + (w(x-2) + w(x-5) + w(x+1))*0.6)
set contour
set cntrparam levels 28
set cntrparam levels incremental -0.25, 0.05, 1.15
set contour base
set nosurface
set view 0,0
set isosample 250
set sample 250, 250
set key at graph 1.25, graph 0.95
unset ztics
#set noclabel
#load "~/term_eps_bw"
set term postscript eps font ",14" size 6,6
set output "n1_q2.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n3_q4(x)
set output "n3_q4.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n5_q4(x)
set output "n5_q4.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n3_q6(x)
set output "n3_q6.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n5_q6(x)
set output "n5_q6.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n7_q6(x)
set output "n7_q6.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n9_q6(x)
set output "n9_q6.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n3_q8(x)
set output "n3_q8.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n5_q8(x)
set output "n5_q8.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n7_q8(x)
set output "n7_q8.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n9_q8(x)
set output "n9_q8.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n11_q8(x)
set output "n11_q8.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output
w(x) = omega_n13_q8(x)
set output "n13_q8.eps"
splot [-0.5:3.5][-0.5:3.5][] f(x,y) notitle
set output