ism3d.maths.stats.custom_pdf

ism3d.maths.stats.custom_pdf(func, x, sersic_n=1)

note:

the pdf values is not normalized: int(pdf*dx) is not one!

Although we can use scipy.stats, we write a customized code here to improve performance.