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

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