function out = MGrnd(x)
a = 1-cosh(x)+x.*sinh(x);
b = sqrt(1-exp(-x));

out = a./((x.^2).*(b.^3));
