function [BmRecr]=recrrnd

mu=4;
sigma=8;
mu1=log(mu^2/sqrt(sigma^2+mu^2));
sigma1=sqrt(log((sigma^2+mu^2)/mu^2));
    
BmRecr=lognrnd(mu1,sigma1);