fixed bug in bivariate normal

jdh2358 · jdh2358 · commit 18fb2085746c · 2006-07-09T20:57:27.000Z
svn path=/trunk/matplotlib/; revision=2554
diff --git a/lib/matplotlib/_pylab_helpers.py b/lib/matplotlib/_pylab_helpers.py
@@ -14,6 +14,7 @@ def get_fig_manager(num):
     get_fig_manager = staticmethod(get_fig_manager)
 
     def destroy(num):
+        
         if not Gcf.has_fignum(num): return
         figManager = Gcf.figs[num]
 
diff --git a/lib/matplotlib/mlab.py b/lib/matplotlib/mlab.py
@@ -960,8 +960,8 @@ def bivariate_normal(X, Y, sigmax=1.0, sigmay=1.0,
     Ymu = Y-muy
 
     rho = sigmaxy/(sigmax*sigmay)
-    z = Xmu**2/sigmax**2 + Ymu**2/sigmay - 2*rho*Xmu*Ymu/(sigmax*sigmay)
-    return 1.0/(2*pi*sigmax*sigmay*(1-rho**2)) * exp( -z/(2*(1-rho**2)))
+    z = Xmu**2/sigmax**2 + Ymu**2/sigmay**2 - 2*rho*Xmu*Ymu/(sigmax*sigmay)
+    return 1.0/(2*pi*sigmax*sigmay*numerix.mlab.sqrt(1-rho**2)) * exp( -z/(2*(1-rho**2)))
 
 
 
