import mir.math.common: approxEqual;

assert(0.5.betaCDF(1, 1).approxEqual(0.5));
assert(0.75.betaCDF(1, 2).approxEqual(0.9375));
assert(0.25.betaCDF(0.5, 4).approxEqual(0.8588867));