Sep-26-2018, 05:29 PM

The integral is supposed to return one but it returns zero.

https://quantummechanics.ucsd.edu/ph130a...de233.html

The integral I'm doing is shown in the link above in the section that says "The radial wavefunctions should be normalised as below." I'm using R_10 which is the ground state.

https://quantummechanics.ucsd.edu/ph130a...de233.html

The integral I'm doing is shown in the link above in the section that says "The radial wavefunctions should be normalised as below." I'm using R_10 which is the ground state.

import scipy from scipy import constants import math import numpy as np from numpy.polynomial import Laguerre from scipy import special import sympy from sympy import * from sympy import var from sympy import integrals from scipy import integrate #constants a_0 = scipy.constants.value('Bohr radius') e = scipy.constants.value('elementary charge') hbar = scipy.constants.value('Planck constant over 2 pi') a = a_0 def int3(x): c1 = (2/a)**(1.5) y1 = c1*np.exp(-x/a) m = np.conjugate(y1)*y1*x**2 return m I3= integrate.quad(int3, 0, np.inf) print(I3)Thank you in advance