When mathematicians encounter the expression $e^{xy}$, the question of how to integrate it arises frequently in calculus courses and applied sciences. That said, the integral of $e$ to the $xy$ is not a single fixed formula; rather, the solution depends entirely on which variable serves as the integration variable and whether the other letters represent constants or functions. Understanding this distinction unlocks the ability to handle everything from basic antiderivatives to complex double integrals in physics and engineering And that's really what it comes down to..
Integration with Respect to $x$
Treating $y$ as a constant, the integral of $e^{xy}$ with respect to $x$ follows a straightforward substitution pattern. Set $u = xy$, which implies $du = y,dx$, or equivalently $dx = \frac{du}{y}$. Substituting these into the integral yields:
$\int e^{xy},dx = \int e^u \cdot \frac{du}{y} = \frac{1}{y}e^u + C = \frac{e^{xy}}{y} + C$
This result holds whenever $y \neq 0$. Because of that, to verify, differentiate $\frac{e^{xy}}{y}$ with respect to $x$ using the chain rule. The derivative of the outer exponential brings down the factor $y$, which cancels the $\frac{1}{y}$ coefficient, leaving exactly $e^{xy}$. Always include the constant of integration $C$ when writing an indefinite integral, because the derivative of any constant vanishes.
Integration with Respect to $y$
By symmetry, integrating with respect to $y$ while holding $x$ constant produces an analogous result. Using the substitution $v = xy$ gives $dv = x,dy$, so:
$\int e^{xy},dy = \frac{1}{x}e^{xy} + C$
provided $x \neq 0$. Students sometimes mistakenly apply the power rule here, treating $e^{xy}$ as though it were $x^n$ or $y^n$. Now, remember that the exponential function is its own derivative up to a chain-rule factor, so its integral mirrors that behavior. If the exponent were merely $x$ or $y$ alone, the denominator would simply be $1$; the presence of the product $xy$ introduces the reciprocal of the constant coefficient.
Double Integrals Over Rectangular Regions
When both variables appear in the limits or the integrand itself, a double integral becomes necessary. For a rectangular region $[a,b] \times [c,d]$, the iterated integral takes the form:
$\int_c^d \int_a^b e^{xy},dx,dy$
The inner integral with respect to $x$ produces $\frac{e^{xy}}{y}\Big|_{x=a}^{x=b} = \frac{e^{by} - e^{ay}}{y}$. The outer integral then becomes:
$\int_c^d \frac{e^{by} - e^{ay}}{y},dy$
This outer integral often lacks an elementary closed form, requiring numerical methods or special functions such as the exponential integral $\text{Ei}(z)$. In such cases, switching the order of integration or employing series expansions may offer a more tractable path. Fubini’s theorem guarantees that for continuous functions on rectangular domains, the order of integration can be swapped without changing the result, though one order might simplify the arithmetic significantly Less friction, more output..
Advanced Techniques and Special Cases
Beyond simple substitution, several techniques address more complicated scenarios involving $e^{xy}$.
Integration by Parts: When the integrand includes a polynomial factor such as $x e^{xy}$, integration by parts becomes essential. Let $u = x$ and $dv = e^{xy},dx$, so $du = dx$ and $v = \frac{e^{xy}}{y}$. Applying $\int u,dv = uv - \int v,du$ yields:
$\int x e^{xy},dx = \frac{x e^{xy}}{y} - \int \frac{e^{xy}}{y},dx = \frac{x e^{xy}}{y} - \frac{e^{xy}}{y^2} + C$
Series Expansion: If $xy$ remains small or if closed-form antiderivatives prove elusive, expand the exponential as a Taylor series:
$e^{xy} = \sum_{n=0}^{\infty} \frac{(xy)^n}{n!}$
Integrating term by term
Integrating term‑by‑term with respect to (x) gives
[ \int e^{xy},dx = \sum_{n=0}^{\infty}\frac{(xy)^{n}}{n!So }\int x^{n},dx = \sum_{n=0}^{\infty}\frac{(xy)^{n}}{n! Which means },\frac{x^{,n+1}}{n+1}+C = \sum_{n=0}^{\infty}\frac{x^{,n+1}y^{,n}}{(n+1)! }+C Small thing, real impact..
If the outer integral is taken over a finite interval ([c,d]) in the variable (y), the same reasoning applied to the inner integral with respect to (y) yields
[ \int_{c}^{d} e^{xy},dy = \sum_{n=0}^{\infty}\frac{x^{,n}}{n!},\frac{d^{,n+1}-c^{,n+1}}{n+1}+C . ]
So naturally, the double integral over the rectangle ([a,b]\times[c,d]) can be written as a double series
[ \int_{c}^{d}!\int_{a}^{b} e^{xy},dx,dy = \sum_{n=0}^{\infty}\frac{1}{n!}, \left(\frac{b^{,n+1}-a^{,n+1}}{n+1}\right) \left(\frac{d^{,n+1}-c^{,n+1}}{n+1}\right) . ]
Although this series converges for all real (a,b,c,d), it is rarely used in practice because the number of terms required for high accuracy can be large. In most applications one resorts to numerical quadrature or to the exponential integral (\operatorname{Ei}), which emerges when the inner integral (\int \frac{e^{by}-e^{ay}}{y},dy) is evaluated directly.
Summary
- The indefinite integral of (e^{xy}) with respect to one variable is elementary: (\int e^{xy},dx = \frac{e^{xy}}{y}+C) (for (y\neq0)) and (\int e^{xy},dy = \frac{e^{xy}}{x}+C) (for (x\neq0)).
- When both variables appear in the limits, a double integral over a rectangular region is required. The inner integration yields a quotient (\frac{e^{by}-e^{ay}}{y}), and the resulting outer integral typically does not admit a simple closed form.
- Fubini’s theorem permits swapping the order of integration, which may simplify the computation; in some configurations the swapped order leads to an integrand that is more readily handled by elementary antiderivatives.
- For cases where an elementary antiderivative is unavailable, a Taylor‑series expansion of the exponential provides a systematic way to express the result as an infinite sum, though computational effort must be balanced against desired precision.
- Special functions such as the exponential integral or numerical integration routines are the standard tools for evaluating the outer integral in practical scenarios.
To wrap this up, while the single‑variable antiderivatives of (e^{xy}) are straightforward, the presence of two variables in the integration limits introduces significant complexity. Mastery of substitution, integration by parts, series expansion, and the flexibility offered by Fubini’s theorem equips the reader to tackle both indefinite and definite double integrals involving the exponential of a product.
