"Fine structure 2" "Verify equation (1)" En = -mu c^2 alpha^2 / (2 n^2) Er = -En^2 / (2 mu c^2) (4 n / (l + 1/2) - 3) Eso = n En^2 / (mu c^2) (j (j + 1) - l (l + 1) - 3/4) / (l (l + 1/2) (l + 1)) T = En + Er + Eso - En (1 + alpha^2 / (n (j + 1/2)) - 3 alpha^2 / (4 n^2)) check(eval(T, j, l + 1/2) == 0) check(eval(T, j, l - 1/2) == 0) check(eval(T, l, j + 1/2) == 0) check(eval(T, l, j - 1/2) == 0) "ok" Run