import mir.math.common: log;
import mir.test: shouldApprox;

0.50.fLPDF(1, 1).shouldApprox == log(0.3001054);
0.75.fLPDF(1, 2).shouldApprox == log(0.2532039);
0.25.fLPDF(0.5, 4).shouldApprox == log(0.4904035);
0.10.fLPDF(2, 1).shouldApprox == log(0.7607258);
0.00.fLPDF(1, 3).shouldApprox == double.infinity;