integrate tan(theta)*legendreP(1,rcos(theta))r^2 sin(theta) dr dphi dtheta, r = 0 to R, phi = 0 to 2*Pi, theta = 0 to Pi