自定义PDF的规范化

Question

我需要实现一些概率密度函数li：

import numpy as np

def bekkers(x, a, m, d):
    P = a*np.exp((-1*(x**(1/3) - m)**2)/(2*d**2))*np.power(x, -2/3)
    return(P)

问题是，它没有内在的标准化，所以必须手工标准化。我尝试使用scipy.integrate.quad来计算从0到无穷大的函数的积分，并将其用作P的分母，但它非常慢并且当我尝试将函数用于非线性时它会干扰优化配件。最优雅，最快速的方法是什么？