Bottom Page

• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 Hydrogen ground state integral Eduard Unladen Swallow Posts: 4 Threads: 2 Joined: Aug 2018 Reputation: 0 Likes received: 0 #1 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. ```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 Larz60+ aetate et sapientia Posts: 9,180 Threads: 331 Joined: Sep 2016 Reputation: 398 Likes received: 1202 #2 Sep-26-2018, 08:07 PM (This post was last modified: Sep-26-2018, 08:08 PM by Larz60+. Edited 2 times in total.) is it really zero, or .99999999999999999999999991 ? Eduard Unladen Swallow Posts: 4 Threads: 2 Joined: Aug 2018 Reputation: 0 Likes received: 0 #3 Sep-26-2018, 10:19 PM The output should be 1 or close to one like 0.99999999... for all values of 'a' in the code. But as i make a smaller in this equal to the bohr radius a_0 it becomes zero. « Next Oldest | Next Newest »

Top Page

 Possibly Related Threads... Thread Author Replies Views Last Post finding the integral of probability density function Staph 3 561 Aug-11-2019, 09:19 AM Last Post: buran

Forum Jump:

Users browsing this thread: 1 Guest(s)