Чем scipy.stats.multivariate_normal.pdf отличается от той же функции, написанной с использованием numpy?

Question

Мне нужно использовать многомерное нормальное распределение в скрипте. Я заметил, что моя версия дает ответ, отличный от метода Сципи. Я не могу понять, почему ...

Вот моя функция:

def gauss(x, mu, sigma):
    assert np.linalg.det(sigma)!=0, "determinant of sigma is 0"
    y = np.exp((-1/2)*(x-mu).T.dot(np.linalg.inv(sigma)).dot(x-mu))/np.sqrt(
      np.power(2*np.pi, len(x))*np.linalg.det(sigma)
    )
    return y

Вот сравнение результатов:

from scipy.stats import multivariate_normal
import numpy as np

x = np.array([-0.54849176, 6.39530657])
mu = np.array([15,20])
sigma = np.array([
  [2,3],
  [4,10]
])

print(gauss(x, mu, sigma))
# output is 1.8781656851138248e-37

print(multivariate_normal.pdf(x, mu, sigma))
# output is 2.698549423643947e-61

Есть кто-нибудьзаметил это? Моя функция неверна? Любая помощь будет принята с благодарностью!