# generic commands for twographbox output size
load "twograph3Dsquare.in"



set data style lines

set zlabel "Density" 0,-0.3
set xrange [0:1]
set yrange [0:1]
set nokey
set samples 200,200
set isosamples 40,40
set hidden3d

## axis-parallel 2D Gaussian
## n(x,y,mx,my,sx,sy) = (1/(sx*sy*2*pi))*exp(-(1/2)*( ((x-mx)**2)/(sx**2) + ((y-my)**2)/(sy**2)))

## Standardized normal coordinates:
sn(x,mx,sx) = (x - mx)/sx
n(nx,ny,sx,sy) = (1/(sx*sy*2*pi))*exp((-0.5)*(nx**2 + ny**2))

w1 = 0.2
m1x = 0.2
m1y = 0.75
s1x = 0.1
s1y = 0.1

w2 = 0.3
m2x = 0.5
m2y = 0.5
s2x = 0.1
s2y = 0.05

w3 = 0.5
m3x = 0.75
m3y = 0.75
s3x = 0.07
s3y = 0.07



set output "mixture-3D.eps"
#splot exp((-0.5)*(sn(x,m1x,s1x)**2 + sn(y,m1y,s1y)**2))
splot w1*n(sn(x,m1x,s1x),sn(y,m1y,s1y),s1x,s1y) + w2*n(sn(x,m2x,s2x),sn(y,m2y,s2y),s2x,s2y) +  w3*n(sn(x,m3x,s3x),sn(y,m3y,s3y),s3x,s3y)