from scipy.stats import logistic
mean = 0.5
scale = 1.2
x = 0.8
pdf = logistic.pdf(x, loc=mean, scale=scale)
print("PDF at", x, ":", pdf)
cdf = logistic.cdf(x, loc=mean, scale=scale)
print("CDF at", x, ":", cdf)
size = 1000
random_numbers = logistic.rvs(loc=mean, scale=scale, size=size)