P=(2J+1)exp(-(BJ(J+1))/kT) where B, k and