import mir.test: shouldApprox;

0.3918266.fInvCDF(1, 1).shouldApprox == 0.50; 
0.522233.fInvCDF(1, 2).shouldApprox == 0.75;
0.5183719.fInvCDF(0.5, 4).shouldApprox == 0.25;
0.08712907.fInvCDF(2, 1).shouldApprox == 0.10;
0.0.fInvCDF(1, 1).shouldApprox == 0;
1.0.fInvCDF(1, 1).shouldApprox == double.infinity;