scipy.integrate error:only size-1 arrays can be converted to Python scalars

Question

I am trying to calculate probabilities from Gaussian Distribution by varying the standard deviation(std). I am expecting by using the Gaussian Quadrature with 21 points, int