import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import lognorm
from scipy.stats import norm
# Compare lognormal and gaussian
mean = 5
sigma = 4
x = np.linspace( max(0, mean-5*sigma), mean+5*sigma, 1000)
ax=plt.subplot(111)
ax.plot(x, lognorm.pdf(x, np.log(1+sigma/mean), scale=mean), label="lognormal")
ax.plot(x, norm.pdf(x, scale=sigma, loc=mean), label="gaussian")
ax.grid()
ax.legend()
ax.set_xlim([x[0], x[-1]])
_ = ax.set_ylim(bottom=0)