import mir.test: shouldApprox;

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