第一问
using TyOptimization
root2d = x-> begin
f1=x[1]^2-x[2]-2
f2=x[2]^2-2*x[1]-4
return [f1,f2]
end

a0=[2,2]
sol,fval=fsolve(root2d, a0)
x0=sol[1]
y0=sol[2]

得到答案：
x=2.2143197433781077
y=2.903211925912013

第二问
用三重积分算出质量即可
m=∭V​ρ(x,y,z)dV=∭V​(2+z)dV

using TyMath
fun(z,y,x)=2+z
xmin=-1
xmax=1
ymin(x)=-sqrt(1-x^2)
ymax(x)=sqrt(1-x^2)
zmin(x,y)=sqrt(x^2+y^2)
zmax(x,y)=4-x^2-y.^2
q,e=integral3(fun,xmin,xmax,ymin,ymax,zmin,zmax)

答案是36.39011486646299