syms x y z
%3a
expr_A = x*exp(-y);
expr_A_y = int(expr_A,y,0,log(x));
expr_A_yx = int(expr_A_y,x,0,2*z);
expr_A_yxz = int(expr_A_yx,z,1,2);
integral_a = vpa(expr_A_yxz);
% display results
disp('"Part A"')
"Part A"
fprintf('The result of integral is: %f\n\n',integral_a);
The result of integral is: 1.666667
%3b
expr_B = z/(y+1);
expr_B_x = int(expr_B,x,0,sqrt(1-z^2));
expr_B_xz = int(expr_B_x,z,0,1);
expr_B_xzy = int(expr_B_xz,y,0,1);
integral_b = vpa(expr_B_xzy);
% display results
disp('"Part B"')
"Part B"
fprintf('The result of integral is: %f\n\n',integral_b);
The result of integral is: 0.231049