code: def kl_mc(self, other, size=500000):
      fun = lambda x: self.logpdf(x) - other.logpdf(x)
      rvs = self.rvs(size=size)
      return fun(rvs).mean()
  
imports: import numpy
  from scipy import special
  from statsmodels.sandbox.distributions.multivariate import mvstdtprob
  from .extras import mvnormcdf
label_desc: Use image node_id 8 for calling the BivariateNormal obj's underlying member method code with example usage: obj.kl_mc(other, size) and returns: fun