stats.norm.cdf(2,loc=mean, scale=std_dev)

#calculating the probability or the area under curve to the left of this z value
import scipy.stats as stats
stats.norm.pdf(x, loc=mean, scale=std_dev)




# The probability (area) to the right is calculated as (1 - probability to the left)
import scipy.stats as stats
1 - stats.norm.pdf(x, loc=mean, scale=std_dev)