Function rRayleigh

double rRayleigh(RGen) (
  double mode,
  ref RGen gen = rndGen
);