The Qlib Library: function Berli MMA implementation

Berli[n_, a_] :=
  Module[{e = 1}, Do[e = 1 + i*e/a, {i, 1, n}]; 1/e ]