"Atomic transitions 6" Rab = pi e^2 / (3 epsilon0 hbar^2) Iab^2 rho Rba = Rab Exp = 1 / (exp(hbar omega0 / (k T)) - 1) rho = hbar omega0^3 / (pi^2 c^3) Exp "Verify equation (1)" check(Rab == e^2 omega0^3 / (3 pi epsilon0 hbar c^3) Iab^2 Exp) "ok" "Verify equation (2)" Aba = Rba (exp(hbar omega0 / (k T)) - 1) check(Aba == e^2 omega0^3 / (3 pi epsilon0 hbar c^3) Iab^2) "ok" "Verify dimensions" joule = kilogram meter^2 / second^2 e = coulomb omega0 = 1/second epsilon0 = coulomb^2 / (joule meter) hbar = joule second c = meter / second Iab = meter check(e^2 omega0^3 / (epsilon0 hbar c^3) Iab^2 == 1/second) "ok" Run