using TyOptimization
F = z -> [z[1]^2 - z[2] - 2, z[2]^2 - 2*z[1] - 4]
x0 = [2.0, 2.0]
x, _ = fsolve(F, x0)
println(x)


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
q, e = integral3(fun, xmin, xmax, ymin, ymax, zmin, zmax)
println("零件总质量: ", q, " g")
println("积分误差: ", e)