[debateshang]

% question1
clear all
clc
syms x y z
f=x^2+y^2+z^2;
F=int(int(int(f,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),1,2);
answer1=vpa(F,a);
% answer
%F=1610027357/6563700+14912/4641*2^(1/4)-6072064/348075*2^(1/2)+64/225*2^(3/4)