with(intpakX):
read "Solve_intpakX.m":
`ALIAS/fepsilon`:=1e-30:
GradientSolve_intpakX([x*cos(y)+x*y+1,y*sin(y)+x-2],[x,y],[[
-10,10],[-10,10]]);