plot3d(Heaviside(exp(-(y+x^3)^2)-y**2/32-exp(sin(y))),x=-1..1.5,y=-3..0);with(plots);implicitplot(exp(-(y+x^3)^2)-y^2/32-exp(sin(y)),x=-1.5..2,y=-3..0);Int(Heaviside(exp(-(y+x^3)^2)-y^2/32-exp(sin(y))),[x=-1..1.5,y=-3..0],'method'=_MonteCarlo,'epsilon'=0.5e-4);evalf(%);