I would like to determine the range (x-axis) spanning 68% of the distribution data for a numpy array:
This prints the values of +/- sigma, but what is the best way to out
