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