# 
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
# 
#-----

m1 = 1
sigma1 = 1./2

coeff_lbeta  = -3./(8*sigma1**2)
coeff_lalpha =  7.*m1/(8*sigma1**2)

l_alpha = seq( -3, +3, length.out=100 )
l_beta = l_alpha

z = outer( coeff_lalpha*l_alpha, coeff_lbeta*l_beta, FUN="+" )

gamma = 4

#postscript("../../WriteUp/Graphics/Chapter2/prob_2.2.19.eps", onefile=FALSE, horizontal=FALSE)

contour(l_alpha,l_beta,z, levels=c(gamma), xlab="L_ALPHA", ylab="L_BETA", col="black")

#dev.off()