; input
;   x - array of [n,2] values. x[*,0] are x values, x[*,1] are y values
;   p - gaussian curve parameters
;    p[0] - mean in x direction
;    p[1] - st dev in x direction
;    p[2] - "area" more closely related to volume under surface
;    p[3] - mean in y direction
;    p[4] - st dev in y direction
function gauss2,x,p
  p1=p[0:2] ;mean, stdev, area of curve 1
  p2=[p[3],p[4],p[2]] ;mean, stdev, area of curve 2
  a=gauss1(x[*,0],p1)
  b=gauss1(x[*,1],p2)

  return,a*b
end