Issue 36113: Problem With SciPy Computation of sigma
"""
Construct normal probability plot
"""
import math
import scipy.stats as ss
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
import prettytable
from random import seed
seed(100)
################
mu = 0.0
sigma = 1.0
x = np.linspace(mu - 5*sigma, mu + 5*sigma, 10)
norm_pdf_x= ss.norm.pdf(x,mu,sigma)
plt.plot(x,norm_pdf_x)
print("The Std. Dev. of normally Dist. x =", np.std(norm_pdf_x))
print()
###########################################
Produces a wrong Result for sigma:
The Std. Dev. of normally Dist. x = 0.13142071847657633
Should be 1.0