function[MrtS]=mrtsrnd

mu=27;
sigma=13;
mu1=log(mu^2/sqrt(sigma^2+mu^2));
sigma1=sqrt(log((sigma^2+mu^2)/mu^2));

Mrtperc=101;
while Mrtperc>100
    Mrtperc=lognrnd(mu1,sigma1);
end

MrtS=1-Mrtperc/100;