The only safe and sure way to convert the iterated integral ${\int }_a^b{\int }_{y\=y_L\left(x\right)}^{y\=y_T\left(x\right)}f\left(x\,y\right) \mathrm{dy} \mathrm{dx}$ to an equivalent iterated integral in the opposite order, that is, to ${\int }_c^d{\int }_{x\=x_L\left(y\right)}^{x\=x_R\left(y\right)}f\left(x\,y\right) \mathrm{dx} \mathrm{dy}$, is to infer the region of integration from the first form of the integral, draw that region, and from that sketch, build the second form of the integral with the order of integration reversed. There is no automatic way to accomplish this reversal.

There are two Maple tools that are useful in deducing the region of integration from a given iterated double integral. The first is a task template for visualizing regions of integration; the second, the inequal command in the plots package. Both tools are demonstrated in the examples below.

Sometimes reversing the order of integration makes it possible to evaluate an iterated integral that otherwise could not be evaluated in closed form. Sometimes reversing the order of integration will mean that the sum of two or more iterated integrals is the the equivalent of one iteration taken in the opposite order.