int [arctan(1/sqrt 2) to pi/4] int [cot(theta) csc(theta) to 2 sec(theta)] r^3 dr dtheta