Changing the order of integration: surefire method?

[v0id]

Is there a way to transform the limits of integration for a multivariable integral without appealing to geometrical manipulations? For example:
$$\int_a^b \int_{y_1(x)}^{y_2(x)} \int_{z_1(x,y)}^{z_2(x,y)} f(x,y,z) \; dz \; dy\; dx \rightarrow \int_c^d \int_{y_3(z)}^{y_4(z)} \int_{x_3(y,z)}^{x_4(y,z)} f(x,y,z) \; dx \; dy \; dz$$
But how does one determine the limits for the second integral (most importantly, the functions $$x_3(y,z)$$ & $$x_4(y,z)$$, etc.) without actually resorting to an illustration of the function and its integration subset? Clearly illustrations are of no use in higher dimensions. Are there any algorithms that generate these interval transformations?P.S: For simplicity, let us assume that $$f(x,y,z)$$ is continuous, bounded, and generally "nice".

 

[Hurkyl]

The region of integration is simply the set defined by the system of inequalities

$$a \leq x \leq b,$$
$$y_1(x) \leq y \leq y_2(x),$$
$$z_1(x, y) \leq z \leq z_2(x, y).$$

So, all you need to do is to formally manipulate these inequalities to put them into the form you want in the new integral.

I.E. You would first solve this system of inequalities for z, then for y in terms of z, then for x.

 

[gammamcc]

I would not expect a simple general algorithm since you may be forced to split the integral into several regions: That can happen say when one of the z_i(x,y),y_j(x) functions is not one-to-one.

 

[v0id]

Thanks Hurkyl and gammamcc.