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}
0101ex+ydydx=(e1)22.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
0101exeydydx=(01exdx)(01eydy)=(e1)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 e11.71828e - 1 \approx 1.71828, so the answer is about 1.7182822.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+ydy=[ex+y]y=0y=1=ex+1ex=ex(e1)\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(e1)dx=(e1)[ex]01=(e1)2\int_{0}^{1} e^{x}(e-1)\,dx = (e-1)\left[e^x\right]_{0}^{1} = (e-1)^2

The constant e1e - 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+ydy=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+ydy=ex+y\int e^{x+y}\,dy = e^{x+y}, since y(x+y)=1\frac{\partial}{\partial y}(x+y) = 1.
  • ex+3ydy=ex+3y3\int e^{x+3y}\,dy = \frac{e^{x+3y}}{3}, since the derivative is 3.
  • exydy=exyx\int e^{xy}\,dy = \frac{e^{xy}}{x} for x0x \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 sinxcosy+cosxsiny\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 (e21)(e31)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).