"Scattering part 4" check((1 - cos(theta))^2 == 4 sin(theta/2)^4) "Verify cross section formula" f = -2 m beta / (p^2 + epsilon^2 hbar^2) epsilon = 0 beta = -Z alpha hbar c p = sqrt(4 m E (1 - cos(theta))) dsigma = f conj(f) check(dsigma == (Z alpha hbar c / (2 E))^2 / (4 sin(theta/2)^4)) "ok" Run