docstr: Kullback-Leibler divergence between this and another distribution
  
  int f(x) (log f(x) - log g(x)) dx
  
  where f is the pdf of self, and g is the pdf of other
  
  uses double integration with scipy.integrate.dblquad
  
  limits currently hardcoded
code: def kl(self, other):
      """Kullback-Leibler divergence between this and another distribution
  
      int f(x) (log f(x) - log g(x)) dx
  
      where f is the pdf of self, and g is the pdf of other
  
      uses double integration with scipy.integrate.dblquad
  
      limits currently hardcoded
  
      """
      fun = lambda x: self.logpdf(x) - other.logpdf(x)
      return self.expect(fun)
  
imports: import numpy
  from scipy import special
  from statsmodels.sandbox.distributions.multivariate import mvstdtprob
  from .extras import mvnormcdf
label_desc: Use image node_id 7 for calling the BivariateNormal obj's underlying member method code with example usage: obj.kl(other) and returns: self