Solution: The following answers we derived using Maple. The codes are provided following the iterated integral.
mass =
evalf(int(int(int(x^2+y^2,z=0..sqrt(9-y^2)),y=3*x..3),x=0..1),3);
To find the center of mass we must first find the actual mass. The following is the code used by maple to find the exact mass.
Center of mass
int(int(int(x^2+y^2,z=0..sqrt(9-y^2)),y=3*x..3),x=0..1);
evalf((5/56)*int(int(int(x*(x^2+y^2),z=0..sqrt(9-y^2)),y=3*x..3),x=0..1),3);
evalf((5/56)*int(int(int(y*(x^2+y^2),z=0..sqrt(9-y^2)),y=3*x..3),x=0..1),3);
evalf((5/56)*int(int(int(z*(x^2+y^2),z=0..sqrt(9-y^2)),y=3*x..3),x=0..1),3);
Moment of inertia
evalf(int(int(int((x^2+y^2)*(x^2+y^2),z=0..sqrt(9-y^2)),y=3*x..3),x=0..1),3);
by: gm