"Verify equation (1)" f = exp(i p / hbar r - r / a) - exp(-i p / hbar r - r / a) F = integral(f,r) I = 0 - eval(F,r,0) check(I == -1 / (i p / hbar - 1 / a) - 1 / (i p / hbar + 1 / a)) f = -m Z e^2 / (4 pi epsilon0 hbar^2) hbar / (i p) I check(f == -m Z e^2 / (2 pi epsilon0 (p^2 + (hbar / a)^2))) "ok" "Verify equation (2)" e = sqrt(4 pi epsilon0 alpha hbar c) p = sqrt(4 m E (1 - cos(theta))) check(f == -2 m Z alpha hbar c / (4 m E (1 - cos(theta)) + (hbar/a)^2)) "ok" Run