module Distribution
  module BivariateNormal
    module GSL_
      class <<self
        def pdf(x,y,rho,s1=1.0,s2=1.0)
          GSL::Ran::bivariate_gaussian_pdf(x, y, s1,s2,rho)
        end
      end
    end
  end
end