Multivariable calculus

Double Integral of e^(x + y) Over the Unit Square

The double integral of e to the x plus y over the unit square equals (e minus 1) squared, about 2.95249. Because e to the x plus y equals e to the x times e to the y, the integrand factors, and the answer is the square of the integral of e to the x from 0 to 1, which is e minus 1.

f(x,y)=ex+yf(x,y) = e^{x+y}
∫01∫01ex+y dy dx=(e−1)2≈2.95249\int_{0}^{1}\int_{0}^{1} e^{x+y}\,dy\,dx = (e-1)^2 \approx 2.95249

Numerically 2.952492, confirmed by quadrature on every build.

The exponent law does the factoring for you

A sum in the exponent is a product of exponentials, which is exactly the shape the separable shortcut wants.

ex+y=exeye^{x+y} = e^x e^y
∫01∫01exey dy dx=(∫01ex dx)(∫01ey dy)=(e−1)2\int_{0}^{1}\int_{0}^{1} e^x e^y\,dy\,dx = \left(\int_{0}^{1} e^x\,dx\right)\left(\int_{0}^{1} e^y\,dy\right) = (e-1)^2

Numerically e−1≈1.71828e - 1 \approx 1.71828, so the answer is about 1.718282≈2.952491.71828^2 \approx 2.95249. Both factors are the same integral because the square is symmetric in xx and yy.

The direct route, for when the law is not obvious

You do not need the exponent law. Integrate in yy with xx frozen and the constant exe^x appears on its own.

∫01ex+y dy=[ex+y]y=0y=1=ex+1−ex=ex(e−1)\int_{0}^{1} e^{x+y}\,dy = \left[e^{x+y}\right]_{y=0}^{y=1} = e^{x+1} - e^{x} = e^{x}(e-1)
∫01ex(e−1) dx=(e−1)[ex]01=(e−1)2\int_{0}^{1} e^{x}(e-1)\,dx = (e-1)\left[e^x\right]_{0}^{1} = (e-1)^2

The constant e−1e - 1 pulls out of the outer integral because it holds no xx. Recognising which pieces are constants at each stage is the whole skill in iterated integration.

The mistake: dividing by the exponent

Writing ∫ex+y dy=ex+yx+y\int e^{x+y}\,dy = \frac{e^{x+y}}{x+y} is wrong. You divide by the derivative of the exponent with respect to the integration variable, not by the exponent itself, and here that derivative is 1.

  • ∫ex+y dy=ex+y\int e^{x+y}\,dy = e^{x+y}, since ∂∂y(x+y)=1\frac{\partial}{\partial y}(x+y) = 1.
  • ∫ex+3y dy=ex+3y3\int e^{x+3y}\,dy = \frac{e^{x+3y}}{3}, since the derivative is 3.
  • ∫exy dy=exyx\int e^{xy}\,dy = \frac{e^{xy}}{x} for x≠0x \neq 0, since the derivative is xx.

Compare this with sin⁡(x+y)\sin(x+y), which does not factor at all. Exponentials turn sums into products and sines do not, so the same looking argument gets very different treatment.

Frequently asked questions

Why does e^(x + y) factor when sin(x + y) does not?

The exponent law ea+b=eaebe^{a+b} = e^a e^b splits the exponential into one factor per variable. Sine has no such law: sin⁡(x+y)\sin(x+y) expands to sin⁡xcos⁡y+cos⁡xsin⁡y\sin x\cos y + \cos x \sin y, a sum of two products rather than one product.

What is the value over the rectangle 0 to 2 by 0 to 3?

The same factoring gives (e2−1)(e3−1)≈121.94\left(e^2 - 1\right)\left(e^3 - 1\right) \approx 121.94. The integrand grows fast, so most of that total comes from the corner near (2,3)(2,3).