function [dE] = splitting(n,l,j,coef)

% coef = m*c^2*(Z*alpha)^4

dE = -coef./2 ./(n.^3) .* (1./(j+1/2) - 3./(4.*n));
