matlab,用int法求解

用两次int:
syms x y
>> a=int(int(2*(1+x+y^2),y,0,sqrt(2*x-x^2)),x,0,1)
a =
-2/3+9/8*pi