Function rBeta

double rBeta(RGen) (
  double a,
  double b,
  ref RGen gen = rndGen
);