ln.obj <- LogNormalSimulator$new(mu = 5, sigma = 0.5, n_sample = 100000)
ln.sample <- ln.obj$generate_lognormal_sample()
stats <- ln.obj$compute_stats()
print(stats)