# First define the function we want to integrate.var('x,y')f=x^5*y^3*exp(x*y)# Compute the integral of f(x,y) over R = [0,1]x[0,1] as an iterated integral.# Note the usage: integral(function, variable, lower limit, upper limit)vol=integral(integral(f,y,0,1),x,0,1)vol
21*e - 57
# The decimal expansion to 5 significant digits.vol.n(digits=5)
0.083923
# Here's a plot of the function on R.plot3d(f,(x,0,1),(y,0,1),color='orange')