function lns = lnorms(means,vars,rawelems) % lnorms (mean,vars,rawelems) % Converts standard normal random values (rawelems) to % lognormals with the defined means and variances (means % and vars). For 1 random value you would call: % lnorms(mean, var, randn) % Inputting vectors makes multiple lognormals simultaneously % Note that vars are Variances, not st. deviations nmeans = log(means)-0.5.*log(vars./means.^2 + 1); nvars = log(vars./means.^2 + 1); normals = rawelems.*sqrt(nvars) + nmeans; lns = exp(normals);