"Yukawa potential" "Verify equation (1)" f = -m V0 / (hbar^2 mu) / (i k) * (exp(i k r - mu r) - exp(-i k r - mu r)) I = integral(f,r) f = 0 - eval(I,r,0) check(f == -2 m V0 / (hbar^2 mu) / (k^2 + mu^2)) "ok" "Verify equation (2)" k = sqrt(4 m E (1 - cos(theta)) / hbar^2) check(f == -2 m V0 / mu / (4 m E (1 - cos(theta)) + mu^2 hbar^2)) "ok" "Verify equation (3)" f = 1 / (4 m E y + mu^2 hbar^2)^2 I = defint(f, y, 0, 2) check(I == 2 / (8 m E mu^2 hbar^2 + mu^4 hbar^4)) "ok" Run