int 1/r^(2*m+2)*(2*(r/a)^2-(r/a)^4) dr