import mir.math.common: approxEqual;

assert(0.5.betaProportionCDF(0.5, 2).approxEqual(0.5));
assert(0.75.betaProportionCDF((1.0 / 3), 3).approxEqual(0.9375));
assert(0.25.betaProportionCDF((1.0 / 9), 4.5).approxEqual(0.8588867));