[New-bugs-announce] [issue36113] Problem With SciPy Computation of sigma

[N. Srinivasan]

New submission from N. Srinivasan <nisthesecond at verizon.net>:

"""
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

----------
components: Library (Lib)
files: test_normal.py
messages: 336604
nosy: nisthesecond
priority: normal
severity: normal
status: open
title: Problem With SciPy Computation of sigma
versions: Python 3.7
Added file: https://bugs.python.org/file48170/test_normal.py

_______________________________________
Python tracker <report at bugs.python.org>
<https://bugs.python.org/issue36113>
_______________________________________