View source for Fubini's theorem

{{short description|Conditions for switching order of integration in calculus}}
{{for|the Fubini theorem for [[Meagre set#Terminology|category]]|Kuratowski–Ulam theorem}}
'''Fubini's theorem''' is a [[theorem]] in [[measure theory]] that gives conditions under which a [[double integral]] can be computed as an [[iterated integral]], i.e. by integrating in one variable at a time. Intuitively, just as the volume of a loaf of bread is the same whether one sums over standard slices or over long thin slices, the value of a double integral does not depend on the order of integration when the hypotheses of the theorem are satisfied. The theorem is named after [[Guido Fubini]], who proved a general result in 1907; special cases were known earlier through results such as [[Cavalieri's principle]], which was used by [[Leonhard Euler]].

More formally, the theorem states that if a function is [[Lebesgue integral|Lebesgue integrable]] on a rectangle <math>X\times Y</math>, then one can evaluate the double integral as an iterated integral:<math display="block">\, \iint\limits_{X\times Y} f(x,y)\,\mathrm{d}(x,y) = \int_X\left(\int_Y f(x,y)\,\mathrm{d}y\right)\mathrm{d}x=\int_Y\left(\int_X f(x,y) \, \mathrm{d}x \right)  \mathrm{d}y.</math> 
This formula is generally not true for the [[Riemann integral]] (however, it is true if the function is continuous on the rectangle; in [[multivariable calculus]], this weaker result is sometimes also called Fubini's theorem).

'''Tonelli's theorem''', introduced by [[Leonida Tonelli]] in 1909, is similar but is applied to a non-negative [[measurable function]] rather than to an integrable function over its domain. The Fubini and Tonelli theorems are usually combined and form the '''Fubini–Tonelli theorem''', which gives the conditions under which it is possible to switch the [[Order of integration (calculus)|order of integration]] in an iterated integral.

A related theorem is often called '''Fubini's theorem for infinite series''',<ref>{{citation |last= Tao |first= Terence |date= 2016 |title= Analysis I |page= 188 |publisher= Springer |isbn= 9789811017896 }}</ref> although it is due to [[Alfred Pringsheim]].<ref>{{citation |last1= Whittaker |first1= E. T. |last2= Watson |first2= G. N. |date= 1902 |title= A Course of Modern Analysis |publisher= Cambridge University Press }}</ref> It states that if <math display="inline">\{a_{m,n}\}_{m=1,n=1}^{\infty}</math> is a double-indexed sequence of real numbers, and if <math display="inline">\displaystyle \sum_{(m,n)\in \N\times \N} a_{m,n} </math> is [[absolutely convergent]], then
<math display=block> \sum_{(m,n)\in\N\times \N}a_{m,n} = \sum_{m=1}^\infty\sum_{n=1}^{\infty} a_{m,n} = \sum_{n=1}^\infty \sum_{m=1}^\infty a_{m,n}. </math>

==History==
A special case of Fubini's theorem for continuous functions on the product of closed, bounded subsets of real vector spaces was known to [[Leonhard Euler]] in the 18th century. In 1904, [[Henri Lebesgue]] extended this result to bounded measurable functions on a product of intervals.<ref>{{citation |last= Lebesgue |first= Henri |authorlink= Henri Lebesgue |date= 1904 |title= Leçons sur l'intégration et la recherche des fonctions primitives |trans-title= Lessons on integration and the search for primitive functions |language= fr |publisher= Gauthier-Villars |place= Paris |url= https://archive.org/details/leconegrarecher00leberich}}</ref> Levi conjectured that the theorem could be extended to functions that are integrable rather than bounded{{Citation needed|date=April 2023}} and this was proven by Fubini in 1907.<ref>{{citation |last= Fubini |first= Guido |date= 1907 |authorlink= Guido Fubini |title= Sugli integrali multipli |trans-title= On multiple integrals |language= it |journal= Rom. Acc. L. Rend. (5) |volume= 16 |issue= 1 |pages= 608–614 |jfm= 38.0343.02 }}; reprinted in {{citation |last= Fubini |first= G. |date= 1958 |title= Opere scelte |trans-title= Selected works |volume= 2 |language= it |publisher= Cremonese |pages= 243–249}}</ref> In 1909, Leonida Tonelli gave a variation of the Fubini theorem that applies to non-negative functions rather than integrable functions.<ref>{{citation |last= Tonelli |first= Leonida |date= 1909 |authorlink= Leonida Tonelli |title= Sull'integrazione per parti |trans-title= On integration by parts |journal= Atti della Accademia Nazionale dei Lincei |trans-journal= Proceedings of the National Lincean Academy |language= it |series= (5) |volume= 18 |issue= 2 |pages= 246–253}}</ref>

==Product measures==
If {{mvar|X}} and {{mvar|Y}} are [[measure space]]s, there are several natural ways to define a [[product measure]] on the product {{math|''X'' × ''Y''}}.

In [[Product (category theory)|the sense of category theory]], measurable sets in the product {{math|''X'' × ''Y''}} of measure spaces are the elements of the [[σ-algebra]] generated by the products {{math|''A'' × ''B''}}, where {{mvar|A}} is measurable in {{mvar|X}}, and {{mvar|B}} is measurable in {{mvar|Y}}.

A measure {{mvar|μ}} on {{math|''X'' × ''Y''}} is called a '''{{dfn|product measure}}''' if {{math|1= ''μ''(''A'' × ''B'') = ''μ''{{sub|1}}(''A'')''μ''{{sub|2}}(''B'')}} for measurable subsets {{math|''A'' ⊂ ''X''}} and {{math|''B'' ⊂ ''Y''}} and measures {{math|''μ''{{sub|1}}}} on {{mvar|X}} and {{math|''μ''{{sub|2}}}} on {{mvar|Y}}. In general, there may be many different product measures on {{math|''X'' × ''Y''}}. Fubini's theorem and Tonelli's theorem both require technical conditions to avoid this complication; the most common approach is to assume that all measure spaces are [[σ-finite measure|σ-finite]], in which case there is a unique product measure on {{math|''X'' × ''Y''}}. There is always a unique maximal product measure on {{math|''X'' × ''Y''}}, where the measure of a measurable set is the [[Infimum and supremum|inf]] of the measures of sets containing it that are countable unions of products of measurable sets. The maximal product measure can be constructed by applying [[Carathéodory's extension theorem]] to the additive function {{mvar|μ}} such that {{math|1= ''μ''(''A'' × ''B'') = ''μ''{{sub|1}}(''A'')''μ''{{sub|2}}(''B'')}} on the ring of sets generated by products of measurable sets. (Carathéodory's extension theorem gives a measure on a measure space that in general contains more measurable sets than the measure space {{math|''X'' × ''Y''}}, so strictly speaking, the measure should be restricted to the [[σ-algebra]] generated by the products {{math|''A'' × ''B''}} of measurable subsets of {{mvar|X}} and {{mvar|Y}}.)

The product of two [[complete measure space]]s is not usually complete. For example, the product of the [[Lebesgue measure]] on the unit interval {{mvar|I}} with itself is not the Lebesgue measure on the square {{math|''I'' × ''I''}}. There is a variation of Fubini's theorem for complete measures, which uses the completion of the product of measures rather than the uncompleted product.

==For integrable functions==
Suppose {{mvar|X}} and {{mvar|Y}} are [[sigma-finite measure|σ-finite]] measure spaces and suppose that {{math|''X'' × ''Y''}} is given the product measure (which is unique as {{mvar|X}} and {{mvar|Y}} are σ-finite). Fubini's theorem states that if {{mvar|f}} is {{math|''X'' × ''Y''}} integrable, meaning that {{mvar|f}} is a [[measurable function]] and
<math display="block">\int_{X\times Y} |f(x,y)|\,\mathrm{d}(x,y) < \infty,</math>
then
<math display="block">\int_X\left(\int_Y f(x,y)\,\mathrm{d}y\right)\,\mathrm{d}x = \int_Y\left(\int_X f(x,y)\,\mathrm{d}x\right)\,\mathrm{d}y = \int_{X\times Y} f(x,y)\,\mathrm{d}(x,y).</math>

The first two integrals are iterated integrals with respect to two measures, respectively, and the third is an integral with respect to the product measure. The partial integrals <math display="inline">\int_Y f(x,y)\,\mathrm{d}y</math> and <math display="inline">\int_X f(x,y)\,\mathrm{d}x</math> need not be defined everywhere, but this does not matter as the points where they are not defined form a set of measure 0.

If the above integral of the absolute value is not finite, then the two iterated integrals may have different values {{See below|{{section link||Counterexamples}} below for an illustration of this possibility}}.

The condition that {{mvar|X}} and {{mvar|Y}} are σ-finite is usually harmless because almost all measure spaces for which one wishes to use Fubini's theorem are σ-finite. 
Fubini's theorem has some rather technical extensions to the case when {{mvar|X}} and {{mvar|Y}} are not assumed to be σ-finite.{{sfnp|Fremlin|2003}} The main extra complication in this case is that there may be more than one product measure on {{math|''X'' × ''Y''}}. Fubini's theorem continues to hold for the maximal product measure but can fail for other product measures. For example, there is a product measure and a non-negative measurable function {{mvar|f}} for which the double integral of {{math|{{abs|''f''}}}} is zero but the two iterated integrals have different values {{see below|{{slink||Counterexamples}} below for an example of this}}. Tonelli's theorem and the Fubini–Tonelli theorem ({{xref-print|text= stated below}}) can fail on non σ-finite spaces, even for the maximal product measure.

==Tonelli's theorem for non-negative measurable functions==
'''{{visible anchor|Tonelli's theorem}}''', named after [[Leonida Tonelli]], is a successor of Fubini's theorem. The conclusion of Tonelli's theorem is identical to that of Fubini's theorem, but the assumption that {{math|{{abs|''f''}}}} has a finite integral is replaced by the assumption that {{mvar|f}} is a non-negative measurable function.

Tonelli's theorem states that if {{math|(''X'', ''A'', ''μ'')}} and {{math|(''Y'', ''B'', ''ν'')}} are [[Sigma-finite measure|σ-finite measure spaces]], while <math display="inline">f:X\times Y \to [0,\infty]</math> is a non-negative measurable function, then
<math display="block">\int_X\left(\int_Y f(x,y)\,\mathrm{d}y\right)\,\mathrm{d}x = \int_Y\left(\int_X f(x,y)\,\mathrm{d}x\right)\,\mathrm{d}y = \int_{X\times Y} f(x,y)\,\mathrm{d}(x,y).</math>

A special case of Tonelli's theorem is in the interchange of the summations, as in {{nowrap|<math display="inline">\sum_x \sum_y a_{xy} = \sum_y \sum_x a_{xy}</math>,}} where {{mvar|a{{sub|xy}}}} are non-negative for all {{mvar|x}} and {{mvar|y}}. The crux of the theorem is that the interchange of order of summation holds even if the series diverges. In effect, the only way a change in order of summation can change the sum is when there exist some subsequences that diverge to {{math|+∞}} and others diverging to {{math|−∞}}. With all elements non-negative, this does not happen in the stated example.

Without the condition that the measure spaces are σ-finite, all three of these integrals can have different values. Some authors give generalizations of Tonelli's theorem to some measure spaces that are not σ-finite, but these generalizations often add conditions that immediately reduce the problem to the σ-finite case.  For example, one could take the σ-algebra on {{math|''A'' × ''B''}} to be that generated by the product of subsets of finite measure, rather than that generated by all products of measurable subsets, though this has the undesirable consequence that the projections from the product to its factors {{mvar|A}} and {{mvar|B}} are not measurable. Another way is to add the condition that the support of {{mvar|f}} is contained in a countable union of products of sets of finite measures. D.&nbsp;H. Fremlin gives some rather technical extensions of Tonelli's theorem to some non σ-finite spaces. None of these generalizations have found any significant applications outside of abstract measure theory, largely because almost all measure spaces of practical interest are σ-finite.{{sfnp|Fremlin|2003}}

==Fubini–Tonelli theorem==
Combining Fubini's theorem with Tonelli's theorem gives the '''{{dfn|Fubini–Tonelli theorem}}'''. Often just called Fubini's theorem, it states that if {{mvar|X}} and {{mvar|Y}} are [[σ-finite measure]] spaces, and if {{mvar|f}} is a measurable function, then
<math display="block">\int_X\left(\int_Y |f(x,y)|\,\mathrm{d}y\right)\,\mathrm{d}x=\int_Y\left(\int_X |f(x,y)|\,\mathrm{d}x\right)\,\mathrm{d}y=\int_{X\times Y} |f(x,y)|\,\mathrm{d}(x,y)</math>
Furthermore, if any one of these integrals is finite, then
<math display="block">\int_X\left(\int_Y f(x,y)\,\mathrm{d}y\right)\,\mathrm{d}x=\int_Y\left(\int_X f(x,y)\,\mathrm{d}x\right)\,\mathrm{d}y=\int_{X\times Y} f(x,y)\,\mathrm{d}(x,y).</math>

The absolute value of {{mvar|f}} in the conditions above can be replaced by either the positive or the negative part of {{mvar|f}}; these forms include Tonelli's theorem as a special case as the negative part of a non-negative function is zero and so has finite integral. Informally, all these conditions say that the double integral of {{mvar|f}} is well defined, though possibly infinite.

The advantage of the Fubini–Tonelli over Fubini's theorem is that the repeated integrals of {{math|{{abs|''f''}}}} may be easier to study than the double integral. As in Fubini's theorem, the single integrals may fail to be defined on a measure 0 set.

==For complete measures==
The versions of Fubini's and Tonelli's theorems above do not apply to integration on the product of the real line {{math|ℝ}} with itself with Lebesgue measure.  The problem is that Lebesgue measure on {{math|ℝ × ℝ}} is not the product of Lebesgue measure on {{math|ℝ}} with itself, but rather the completion of this: a product of two complete measure spaces {{mvar|X}} and {{mvar|Y}} is not in general complete. For this reason, one sometimes uses versions of Fubini's theorem for complete measures: roughly speaking, one replaces all measures with their completions. The various versions of Fubini's theorem are similar to the versions above, with the following minor differences:
* Instead of taking a product {{math|''X'' × ''Y''}} of two measure spaces, one takes the completion of the product.
* If {{mvar|f}} is measurable on the completion of {{math|''X'' × ''Y''}} then its restrictions to vertical or horizontal lines may be non-measurable for a measure zero subset of lines, so one has to allow for the possibility that the vertical or horizontal integrals are undefined on a set of measure 0 because they involve integrating non-measurable functions. This makes little difference, because they can already be undefined due to the functions not being integrable. 
* One generally also assumes that the measures on {{mvar|X}} and {{mvar|Y}} are complete, otherwise the two partial integrals along vertical or horizontal lines may be well-defined but not measurable. For example, if {{mvar|f}} is the characteristic function of a product of a measurable set and a non-measurable set contained in a measure 0 set then its single integral is well defined everywhere but non-measurable.

==Proofs==

Proofs of the Fubini and Tonelli theorems are necessarily somewhat technical, as they have to use a hypothesis related to '''σ-finiteness'''. Most proofs involve building up to the full theorems by proving them for increasingly complicated functions, with the steps as follows.
# Use the fact that the measure on the product is multiplicative for rectangles to prove the theorems for the characteristic functions of rectangles.
# Use the condition that the spaces are σ-finite (or some related condition) to prove the theorem for the characteristic functions of measurable sets. This also covers the case of simple measurable functions (measurable functions taking only a finite number of values).
# Use the condition that the functions are measurable to prove the theorems for positive measurable functions by approximating them by simple measurable functions. This proves Tonelli's theorem.
# Use the condition that the functions are integrable to write them as the difference of two positive integrable functions and apply Tonelli's theorem to each of these. This proves Fubini's theorem.

===Riemann integrals===
For [[Riemann integral]]s, Fubini's theorem is proven by refining the partitions along the {{mvar|x}}-axis and {{mvar|y}}-axis as to create a joint partition of the form {{math|[''x''{{sub|''i''}}, ''x''{{sub|''i''+1}}] × [''y''{{sub|''j''}}, ''y''{{sub|''j''+1}}]}}, which is a partition over {{math|''X'' × ''Y''}}. This is used to show that the double integrals of either order are equal to the integral over {{math|''X'' × ''Y''}}.

==Counterexamples==

The following examples show how Fubini's theorem and Tonelli's theorem can fail if any of their hypotheses are omitted.

===Failure of Tonelli's theorem for non σ-finite spaces===

Suppose that {{mvar|X}} is the unit interval with the Lebesgue measurable sets and Lebesgue measure, and {{mvar|Y}} is the unit interval with all the subsets measurable and the [[counting measure]], so that {{mvar|Y}} is not σ-finite. If {{mvar|f}} is the characteristic function of the diagonal of {{math|''X'' × ''Y''}}, then integrating {{mvar|f}} along {{mvar|X}} gives the {{math|0}} function on {{mvar|Y}}, but integrating {{mvar|f}} along {{mvar|Y}} gives the function {{math|1}} on {{mvar|X}}. So, the two iterated integrals are different. This shows that Tonelli's theorem can fail for spaces that are not σ-finite no matter which product measure is chosen. The measures are both [[decomposable measure|decomposable]], showing that Tonelli's theorem fails for decomposable measures (which are slightly more general than σ-finite measures).

===Failure of Fubini's theorem for non-maximal product measures===

Fubini's theorem holds for spaces even if they are not assumed to be σ-finite provided one uses the maximal product measure. In the example above, for the maximal product measure, the diagonal has infinite measure so the double integral of {{math|{{abs|''f''}}}} is infinite, and Fubini's theorem holds vacuously.

However, if we give {{math|''X'' × ''Y''}} the product measure such that the measure of a set is the sum of the Lebesgue measures of its horizontal sections, then the double integral of {{math|{{abs|''f''}}}} is zero, but the two iterated integrals still have different values. This gives an example of a product measure where Fubini's theorem fails.

This gives an example of two different product measures on the same product of two measure spaces. For products of two σ-finite measure spaces, there is only one product measure.

===Failure of Tonelli's theorem for non-measurable functions===

Suppose that {{mvar|X}} is the first uncountable ordinal, with the finite measure where the measurable sets are either countable (with measure 0) or the sets of countable complement (with measure 1). The (non-measurable) subset {{mvar|E}} of {{math|''X'' × ''X''}} given by pairs {{math|(''x'', ''y'')}} with {{math|''x'' < ''y''}} is countable on every horizontal line and has countable complement on every vertical line. If {{mvar|f}} is the characteristic function of {{mvar|E}} then the two iterated integrals of {{mvar|f}} are defined and have different values 1 and 0. The function {{mvar|f}} is not measurable. This shows that Tonelli's theorem can fail for non-measurable functions.

===Failure of Fubini's theorem for non-measurable functions===

A variation of the example above shows that Fubini's theorem can fail for non-measurable functions even if {{math|{{abs|''f''}}}} is integrable and both repeated integrals are well defined: if we take {{mvar|f}} to be 1 on {{mvar|E}} and –1 on the complement of {{mvar|E}}, then {{math|{{abs|''f''}}}} is integrable on the product with integral 1, and both repeated integrals are well defined, but have different values 1 and –1.

[[Wacław Sierpiński]] demonstrated that assuming the continuum hypothesis, one can identify {{mvar|X}} with the unit interval {{mvar|I}}, so there is a bounded non-negative function on {{math|''I'' × ''I''}} whose two iterated integrals (using Lebesgue measure) are both defined but unequal.<ref>{{citation |last= Sierpiński |first= Wacław |authorlink= Wacław Sierpiński |date= 1920 |title= Sur un problème concernant les ensembles mesurables superficiellement |trans-title= On a problem concerning superficially measurable sets |journal= [[Fundamenta Mathematicae]] |lang= fr |volume= 1 |issue= 1 |pages= 112–115 |doi= 10.4064/fm-1-1-112-115 |url= https://eudml.org/doc/212592 |doi-access= free }}</ref> The stronger versions of Fubini's theorem on a product of two unit intervals with Lebesgue measure, where the function is no longer assumed to be measurable but merely that the two iterated integrals are well defined and exist, are independent of the standard [[Zermelo–Fraenkel axioms]] of [[set theory]] with the axiom of choice (ZFC).  The continuum hypothesis and [[Martin's axiom]] both imply that there exists a function on the unit square whose iterated integrals are not equal, while [[Harvey Friedman (mathematician)|Harvey Friedman]] showed that it is consistent with {{abbr|ZFC|Zermelo–Fraenkel set theory with the axiom of choice}} that a strong Fubini-type theorem for {{closed-closed|0, 1}} does hold, and whenever the two iterated integrals exist they are equal.<ref>{{citation |last= Friedman |first= Harvey |authorlink= Harvey Friedman (mathematician) |date= 1980 |title= A Consistent Fubini-Tonelli Theorem for Nonmeasurable Functions |journal= [[Illinois Journal of Mathematics]] |volume= 24 |issue= 3 |pages= 390–395 |doi= 10.1215/ijm/1256047607 |mr= 573474 |url= http://projecteuclid.org/euclid.ijm/1256047607 |doi-access= free }}</ref> {{hatnote inline|See [[List of statements undecidable in ZFC]].}}

===Failure of Fubini's theorem for non-integrable functions===
Fubini's theorem tells us that (for measurable functions on a product of σ-finite measure spaces) if the integral of the absolute value is finite, then the order of integration does not matter; if we integrate first with respect to {{mvar|x}} and then with respect to {{mvar|y}}, we get the same result as if we integrate first with respect to {{mvar|y}} and then with respect to {{mvar|x}}. The assumption that the integral of the absolute value is finite is "[[Lebesgue integral|Lebesgue integrability]]", and without it the two repeated integrals can have different values.

A simple example to show that the repeated integrals can be different in general is to take the two measure spaces to be the positive integers, and to take the function {{math|''f''(''x'', ''y'')}} to be {{math|1}} if {{math|1=''x'' = ''y''}}, {{math|−1}} if {{math|1=''x'' = ''y'' + 1}}, and {{math|0}} otherwise. Then the two repeated integrals have different values {{math|0}} and {{math|1}}.

Another example is as follows for the function
<math display="block">\frac{x^2-y^2}{(x^2+y^2)^2} = -\frac{\partial^2}{\partial x\,\partial y} \arctan(y/x).</math>

The [[double integral|iterated integral]]s
<math display="block">\begin{align}
\int_{x=0}^1\left(\int_{y=0}^1\frac{x^2-y^2}{(x^2+y^2)^2}\,\mathrm{d}y\right)\,\mathrm{d}x &= \frac{\pi}{4} \\
\int_{y=0}^1\left(\int_{x=0}^1\frac{x^2-y^2}{(x^2+y^2)^2}\,\mathrm{d}x\right)\,\mathrm{d}y &=-\frac{\pi}{4}
\end{align}</math>
have different values. The corresponding double integral does not [[absolute convergence|converge absolutely]] (in other words the integral of the [[absolute value]] is not finite):
<math display="block">\int_0^1\int_0^1 \left|\frac{x^2-y^2}{\left(x^2 + y^2\right)^2}\right|\,\mathrm{d}y\,\mathrm{d}x=\infty.</math>

== Fubini's theorem in multiplications of integrals ==

=== Product of two integrals ===

For the product of two integrals with lower limit zero and a common upper limit we have the following formula:
{{Equation box 1 |indent=: |border= 1 |border colour= #a2a9b1 |equation=
<math>\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right] = \int_{0}^{1} \int_{0}^{u} x\,v(xy) \,w(x) + x\,v(x) \,w(xy) \,\mathrm{d}x \,\mathrm{d}y </math>
}}

=== Proof ===

Let {{math|''V''(''x'')}} and {{math|''W''(''x'')}} are primitive functions of the functions {{math|''v''(''x'')}} and {{math|''w''(''x'')}} respectively, which pass through the origin:
<math display="block">\int_{0}^{u} v(x) \,\mathrm{d}x = V(u), 
\quad \quad 
\int_{0}^{u} w(x) \,\mathrm{d}x = W(u).</math>

Therefore, we have
<math display="block">\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right] = V(u) W(u) </math>

By the [[product rule]], the derivative of the right-hand side is
<math display="block">\frac{\mathrm{d}}{\mathrm{d}x} \left[V(x) W(x)\right] = V(x)w(x) + v(x)W(x)</math>

and by integrating we have:
<math display="block">\int_{0}^{u} V(x)w(x) + v(x)W(x) \,\mathrm{d}x=V(u) W(u) </math>

Thus, the equation from the beginning we get:
<math display="block">\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right]= \int_{0}^{u} V(x)w(x) + v(x)W(x) \,\mathrm{d}x </math>

Now, we introduce a second integration parameter {{mvar|y}} for the description of the antiderivatives {{math|''V''(''x'')}} and {{math|''W''(''x'')}}:
<math display="block">\begin{alignat}{3}
\int_{0}^{1} x\,v(xy) \,\mathrm{d}y &= \left[V(xy)\right]_{y = 0}^{y = 1} &&= V(x) \\
\int_{0}^{1} x\,w(xy) \,\mathrm{d}y &= \left[W(xy)\right]_{y = 0}^{y = 1} &&= W(x)
\end{alignat}</math>

By insertion, a double integral appears:
<math display="block">\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right]= \int_{0}^{u} \left[\int_{0}^{1} x\,v(xy) \,\mathrm{d}y\right]w(x) + v(x)\left[\int_{0}^{1} x\,w(xy) \,\mathrm{d}y\right] \,\mathrm{d}x </math>

Functions that are foreign to the concerned integration parameter can be imported into the inner function as a factor:
<math display="block">\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right]= \int_{0}^{u} \left[\int_{0}^{1} x\,v(xy) \,w(x) \,\mathrm{d}y\right] + \left[\int_{0}^{1} x\,v(x)\,w(xy) \,\mathrm{d}y\right] \,\mathrm{d}x </math>

In the next step, the [[Integral#Properties|sum rule]] is applied to the integrals:
<math display="block">\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right]= \int_{0}^{u} \int_{0}^{1} x\,v(xy) \,w(x) + x\,v(x)\,w(xy) \,\mathrm{d}y \,\mathrm{d}x </math>

And finally, we use Fubini's theorem
<math display="block">\left[\int_{0}^{u} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{u} w(x) \,\mathrm{d}x\right]= \int_{0}^{1} \int_{0}^{u} x\,v(xy) \,w(x) + x\,v(x)\,w(xy) \,\mathrm{d}x \,\mathrm{d}y </math>

== Calculation examples ==
{{tone|section|date=April 2024}}

=== Arcsine integral ===

The arcsine integral, also called the inverse sine integral, is a function that cannot be represented by [[elementary function]]s. However, some of the values of the arcsine integral can be expressed with elementary functions. These values can be determined by integrating the derivative of the arcsine integral, which is the quotient of the arcsine divided by the [[identity function]]{{snd}} the cardinalized arcsine.{{Clarify|date=February 2025|reason= Isn't this just the definition of the "arcsine integral" (by analogy with the sine integral)?|text=The arcsine integral is exactly the original antiderivative of the cardinalized arcsine.}} To integrate this function, Fubini's theorem serves as a key, which unlocks the integral by exchanging the order of the integration parameters. When applied correctly, Fubini's theorem leads directly to an antiderivative function that can be integrated in an elementary way (which is shown in cyan in the following equation chain):

<math display="block">\begin{align}
\operatorname{Si}_{2}(1) &= \int_{0}^{1} \frac{1}{x}\arcsin(x) \,\mathrm{d}x = {\color{blue}\int_{0}^{1}} {\color{green}\int_{0}^{1}} \frac{\sqrt{1-x^2}\,y}{(1- x^2 y^2)\sqrt{1-y^2}} \,{\color{green}\mathrm{d}y} \,{\color{blue}\mathrm{d}x} \\
&= {\color{green}\int_{0}^{1}} {\color{blue}\int_{0}^{1}} \frac{\sqrt{1-x^2}\,y}{(1-x^2 y^2)\sqrt{1-y^2}} {\color{blue}\,\mathrm{d}x} {\color{green}\,\mathrm{ d}y} =\int_{0}^{1} \frac{\pi\,y}{2\sqrt{1-y^2}(1+\sqrt{1-y^2}\,)} \,\mathrm{d}y \\
&= {\color{RoyalBlue}\left\{ \frac{\pi}{2} \ln\left[2 \left(1 + \sqrt{1 - y^2}\,\right)^ {-1}\right] \right\}_{y = 0}^{y = 1}} = \frac{\pi}{2}\ln(2)
\end{align}</math>

=== Dirichlet eta Function ===
The [[Dirichlet series]] defines the [[Dirichlet eta function]] as follows:
<math display="block"> \eta(s) = \sum_{n=1}^\infty \frac{(-1)^{n-1}}{n^s} = 1 - \frac{1}{2^ s} + \frac{1}{3^s} - \frac{1}{4^s} + \frac{1}{5^s} - \frac{1}{6^s} \pm \cdots </math>

The value {{math|''η''(2)}} is equal to {{math|{{sfrac|''π''{{sup|2}}|12}}}} and this can be proven with Fubini's theorem{{dubious|| the condition for using Fubini is not fulfilled since |f(x,n)|=\frac{1}{n} {x}^{n-1} is NOT (Lebesgue) integrable because the harmonic series diverges!|date=January 2024}} in this way:
<math display="block"> \eta(2) = \sum_{n = 1}^\infty (-1)^{n-1}\frac{1}{n^2} = \sum_{n = 1}^\infty \int_{0}^{1} (-1)^{n-1}\frac{1}{n} {x}^{n-1} \,\mathrm{d}x = \int_{0 }^{1} \sum_{n = 1}^\infty (-1)^{n-1}\frac{1}{n} {x}^{n-1} \,\mathrm{d}x = \int_{0}^{1} \frac{1}{x}\ln(x+1) \,\mathrm{d}x </math>

The integral of the product of the [[Multiplicative inverse|reciprocal function]] and the [[Logarithm|natural logarithm]] of the ''successor function'' is a [[Polylogarithm|polylogarithmic integral]] and it cannot be represented by elementary function expressions. Fubini's theorem again unlocks this integral in a combinatorial way. This works by carrying out double integration on the basis of Fubini's theorem used on an additive combination of fractionally rational functions with fractions of linear and square denominators:
<math display="block">\begin{align}
&\int_{0}^{1} \frac{1}{x}\ln(x+1) \,\mathrm{d}x = {\color{blue}\int_{0}^{ 1}} {\color{green}\int_{0}^{1}} \frac{4}{3(x^2+2xy+1)} + \frac{2x}{3(x^2y+1 )} - \frac{1}{3(xy+1)} \,{\color{green}\mathrm{d}y} \,{\color{blue}\mathrm{d}x} \\
&\hphantom{\quad\quad\quad} = {\color{green}\int_{0}^{1}} {\color{blue}\int_{0}^{1}} \frac{4}{3(x^2+2xy +1)} + \frac{2x}{3(x^2y+1)} - \frac{1}{3(xy+1)} \,{\color{blue}\mathrm{d}x} \,{\color{green}\mathrm{d}y} \\
&\hphantom{\quad\quad\quad} = \int_{0}^{1} \frac{2\arccos(y)}{3\sqrt{1-y^2}} \,\mathrm {d}y = {\color{RoyalBlue}\left[\frac{\pi^2}{12}-\frac{1}{3}\arccos(y)^2 \right]_{y = 0} ^{y = 1}} \\
&\hphantom{\quad\quad\quad} = \frac{\pi^2}{12}
\end{align}</math>

This way of working out the integral of the ''cardinalized'' natural logarithm{{definition needed|date=August 2025}} of the successor function was discovered by James Harper.<ref>{{cite journal |last= Harper |first= James D. |date= 2003 |title= Another Simple Proof of <math display="inline">1 + \frac{1}{2^{2}} + \frac{1}{3^{2}} + \cdot \cdot \cdot = \frac{\pi^{2}}{6}</math> |journal= The American Mathematical Monthly |volume= 110 |number= 6 |pages= 540–541 |doi= 10.2307/3647912 |jstor= 3647912 }}</ref>

The original antiderivative, shown here in cyan, leads directly to the value of {{math|''η''(2)}}:
<math display="block"> \eta(2) = \frac{\pi^2}{12}</math>

=== Integrals of Complete Elliptic Integrals ===

The [[improper integral]] of the [[Elliptic integral#Complete elliptic integral of the first kind|complete elliptic integral of first kind]], {{math|''K''(''k'')}}, takes the value of twice the [[Catalan constant]] accurately. The antiderivative of that K-integral belongs to the so-called ''elliptic polylogarithms''.{{definition needed|date=August 2025}} The Catalan constant can only be obtained via the [[Inverse tangent integral|arctangent integral]], which results from the application of Fubini's theorem:
<math display="block">\begin{align}
\int_{0}^{1} K(x) \,\mathrm{d}x &= {\color{blue}\int_{0}^{1}} {\color{green}\int_ {0}^{1}} \frac{1}{\sqrt{(1 - x^2 y^2)(1 - y^2)}} \,{\color{green}\mathrm{d}y }\,{\color{blue}\mathrm{d}x} \\
&= {\color{green}\int_{0}^{1}}{\color{blue}\int_{0}^{1}} \frac{1}{\sqrt{(1 - x^2 y^2)(1 - y^2)}} \,{\color{blue}\mathrm{d}x} \,{\color{green} \mathrm{d}y} \\
&= \int_{0}^{1} \frac{\arcsin(y)}{y\sqrt{1 - y^2}} \,\mathrm{d}y \\
&= {\color{RoyalBlue} \left\{ 2\,\mathrm{Ti}_{2} \left[ y\left(1 + \sqrt{1 - y^2}\,\right)^{-1} \right] \right\}_{y = 0}^{y = 1}} \\
&= 2\,\mathrm{Ti}_{2}(1) =2\beta(2) =2\,C
\end{align}</math>

This time, the expression now in royal cyan color tone is not elementary, but it leads directly to the equally non-elementary value of the "Catalan constant" using the arctangent integral, also called ''inverse tangent integral''.

The same procedure also works for the [[Elliptic integral#Complete elliptic integral of the first kind|complete elliptic integral of second kind]], {{math|''E''(''k'')}}, in the following way:
<math display="block">\begin{align}
\int_{0}^{1} E(x) \,\mathrm{d}x &={\color{blue}\int_{0}^{1}} {\color{green}\int_ {0}^{1}} \frac{\sqrt{1 - x^2 y^2}}{\sqrt{1 - y^2}} \,{\color{green}\mathrm{d}y} \,{\color{blue}\mathrm{d}x} \\
&= {\color{green}\int_{0}^{1}}{\color{blue}\int_{0}^{1}} \frac {\sqrt{1 - x^2 y^2}}{\sqrt{1 - y^2}} \,{\color{blue}\mathrm{d}x} \,{\color{green}\mathrm {d}y} \\
&= \int_{0}^{1} \left[\frac{\arcsin(y)}{2y\sqrt{1 - y^2}} + \frac{1}{2}\right] \,\mathrm{d}y \\
&= {\color{RoyalBlue}\left\{ \mathrm{Ti}_{2} \left[ y\left(1 + \sqrt{1 - y^2}\,\right )^{-1} \right] + \frac{1}{2} y \right\}_{y = 0}^{y = 1}} \\
&= \mathrm{Ti}_{2}(1) + \frac{1}{2} =\beta(2) + \frac{1}{2} =C + \frac{1}{2}
\end{align}</math>

=== Double execution for the Exponential Integral Function ===
The [[Euler-Mascheroni constant]] emerges as the [[improper integral]] from zero to infinity at the integration on the product of negative [[natural logarithm]] and the [[Exponential function|exponential reciprocal]]. But it is also the improper integral within the same limits on the ''cardinalized difference'' of the reciprocal of the successor function and the exponential reciprocal:
<math display="block">
\definecolor{cerulean}{rgb}{0.0, 0.48, 0.65}
\gamma = {\color{WildStrawberry}\int_0^\infty \frac{-\ln(x)}{\exp(x)}\,\mathrm{d}x} = {\color{cerulean}\int_{0}^{\infty} \frac{1}{x}\left[\frac{1}{x + 1}-\exp(-x)\right] \,\mathrm{d}x } </math>

The concord of these two integrals can be shown by successively executing Fubini's theorem twice and by leading this double execution of that theorem over the identity to an integral of the complementary [[Exponential integral|exponential integral function]].

The complementary integral exponential function is defined as:
<math display="block">\mathrm{E}_{1}(x) = \exp(-x)\int_{0}^{\infty} \frac{\exp(-xy)}{y+1} \, \mathrm{d}y </math>

The derivative of the complementary integral exponential function is:
<math display="block">\frac{\mathrm{d}}{\mathrm{d}x} \,\mathrm{E}_{1}(x) = -\frac{1}{x}\exp(-x) </math>

First implementation of Fubini's theorem:

This integral from a construction of the integral exponential function leads to the integral from the negative Natural Logarithm and the Exponential Reciprocal:
<math display="block">\begin{align}
\gamma &= {\color{WildStrawberry}\int_{0}^{\infty} -\exp(-y)\ln(y) \,\mathrm{d}y} \\
&= {\color{green}\int_{0}^{\infty}} {\color{blue}\int_{0}^{\infty}} \exp(-y)\left(\frac{1}{x + y} - \frac{1}{x + 1}\right) \,{\color{blue}\mathrm{d}x} \,{\color{green}\mathrm{d}y} \\
&= {\color{blue}\int_{0}^{\infty}}{\color{green}\int_{0}^{\infty}} \exp(-y)\left(\frac{1}{x + y} - \frac{1}{x + 1}\right) \,{\color{green}\mathrm{d}y} \,{\color{blue}\mathrm{d}x} \\
&= {\color{brown}\int_{0}^{\infty} \left[\exp(x)\,\mathrm{E}_{1}(x) - \frac{1}{x + 1}\right] \,\mathrm{d}x}
\end{align}</math>

Second implementation of Fubini's theorem:

The previously described integral from the described cardinalized difference leads to the previously mentioned integral from the exponential integral function:
<math display="block">\definecolor{cerulean}{rgb}{0.0, 0.48, 0.65}
\begin{align}
\gamma &= {\color{brown}\int_{0}^{\infty} \left[\exp(x)\,\mathrm{E}_{1}(x) - \frac{1}{x + 1}\right] \,\mathrm{d}x} \\
&= {\color{blue}\int_{0}^{\infty}}{\color{blueviolet}\int_{0}^{\infty}} \exp(-xz)\left[\frac{1}{z + 1}-\exp(-z)\right] \,{\color{blueviolet}\mathrm{d}z} \,{\color{blue}\mathrm{d}x} \\
&= {\color{blueviolet}\int_{0}^{\infty}}{\color{blue}\int_{0}^{\infty}} \exp(-xz)\left[\frac{1}{z + 1}-\exp(-z)\right] \,{\color{blue}\mathrm{d}x} \,{\color{blueviolet}\mathrm{d}z} \\
&= {\color{cerulean}\int_{0}^{\infty} \frac{1}{z}\left[\frac{1}{z + 1}-\exp(-z)\right] \,\mathrm{d}z}
\end{align}</math>

In principle, products from exponential functions and fractionally rational functions can be integrated like this:

<math display="block">\frac{\exp(-ax)}{bx + c} = \frac{\mathrm{d}}{\mathrm{d}x} \left\{-\frac{1}{b} \exp\left(\frac{ac}{b}\right)\,\mathrm{E}_{1}\left[\frac{a}{b}\left(bx + c\right)\right] \right\} </math>
<math display="block">\int_0^{\infty} \frac{\exp(-ax)}{bx + c} \,\mathrm{d}x = \left\{-\frac{1}{b}\exp \left(\frac{ac}{b}\right)\,\mathrm{E}_{1}\left[\frac{a}{b}\left(bx + c\right)\right]\right \}_{x = 0}^{x = \infty} = \frac{1}{b}\exp\left(\frac{ac}{b}\right) \,\mathrm{E}_{1 }\left(\frac{ac}{b}\right) </math>

In this way it is shown accurately by using Fubini's theorem twice that these integrals are indeed identical to each other.

=== Normal distribution integral ===

Now this formula for the squaring of an integral is set up:
{{Equation box 1 |indent=: |border= 1 |border colour= #a2a9b1 |equation=
<math>\left[\int_{0}^{\infty} f(x) \,\mathrm{d}x\right]^2 = \int_{0}^{1} \int_{0}^{\infty} 2x\,f(x) \,f(xy)\,\mathrm{d}x \,\mathrm{d}y </math>
}}

This chain of equations can then be generated accordingly:
<math display="block">\begin{align}
\left[\int_{0}^{\infty} \exp(-x^2) \,\mathrm{d}x\right]^2 &= \int_{0}^{1} \int_{0}^{\infty} 2x\exp(-x^2)\exp(-x^2 y^2) \,\mathrm{d}x\,\mathrm{d}y \\
&= \int_{0}^{1} \int_{0}^{\infty} 2x\exp\left[-x^2 (y^2 + 1)\right] \,\mathrm{d}x\,\mathrm{d}y \\
&= \int_{0}^{1} \left\{- \frac{1}{y^2 + 1}\exp\left[-x^2(y^2 + 1)\right]\right\}_{x = 0}^{x = \infty} \,\mathrm{d}y \\
&= \int_{0}^{1} \frac{1}{y^2 + 1} \,\mathrm{d}y \\
&= \arctan(1) \\
&= \frac{\pi}{4}
\end{align}</math>

For the integral of the [[normal distribution]] this value can be generated:
<math display="block">\int_{0}^{\infty} \exp(-x^2) \,\mathrm{d}x = \frac{1}{2}\sqrt{\pi} </math>

=== Dilogarithm of one ===

Now another formula for the squaring of an integral is set up again:

{{Equation box 1 |indent=: |border= 1 |border colour= #a2a9b1 |equation=
<math>\left[\int_{0}^{1} g(x) \,\mathrm{d}x\right]^2 = \int_{0}^{1} \int_{0}^{1} 2x\,g(x) \,g(xy)\,\mathrm{d}x \,\mathrm{d}y </math>
}}

So this chain of equations applies as a new example:
<math display="block">\begin{align}
\frac{\pi^2}{4} &= \left(\arcsin(1)\right)^2 = \left[\int_{0}^{1} \frac{1}{\sqrt{1 - x^2}} \,\mathrm{d}x\right]^2 \\
&= \int_{0}^{1} \int_{0}^{1} \frac{2x}{\sqrt{(1 - x^2)(1 - x^2 y^2)}} \,\mathrm{d}x \,\mathrm{d}y \\
&= \int_{0}^{1} \left[\frac{2}{y}\operatorname{artanh}(y) - \frac{2}{y}\operatorname{artanh}\left(\frac{\sqrt{1 - x^2}\,y}{\sqrt{1 - x^2 y^2}}\right)\right]_{x = 0}^{x = 1}\,\mathrm{d}y \\
&= \int_{0}^{1} \frac{2}{y}\operatorname{artanh}(y) \,\mathrm{d}y \\
&= \left[2\,\mathrm{Li}_{2}(y) - \frac{1}{2}\,\mathrm{Li}_{2}(y^2)\right]_{y = 0}^{y = 1} \\
&= \frac{3}{2}\,\mathrm{Li}_{2}(1)
\end{align}</math>

For the [[Dilogarithm]] of one this value appears:
<math display="block">\mathrm{Li}_{2}(1) = \frac{\pi^2}{6} </math>

In this way the [[Basel problem]] can be solved.

=== Legendre's relation ===

In this next example, the more generalized form of the equation is used again as a mold:
{{Equation box 1 |indent=: |border= 1 |border colour= #a2a9b1 |equation=
<math>\left[\int_{0}^{1} v(x) \,\mathrm{d}x\right]\left[\int_{0}^{1} w(x) \,\mathrm{d}x\right]= \int_{0}^{1} \int_{0}^{1} \left(x\,v(xy) \,w(x) + x\,v(x)\,w(xy)\right) \,\mathrm{d}x \,\mathrm{d}y </math>
}}

The following integrals can be computed by using the incomplete [[elliptic integral]]s of the first and second kind as antiderivatives and these integrals have values that can be represented with ''complete elliptic integrals'':
<math display="block">\begin{align}
\int_{0}^{1} \frac{1}{\sqrt{1 - x^4}} \,\mathrm{d}x &= \left\{ - \frac{1}{2}\sqrt{2}\,F\left[\arccos(x);\frac{1}{2}\sqrt{2}\right] \right\}_{x = 0}^{x = 1} \\
&= \frac{1}{2}\sqrt{2}\,K\left(\frac{1}{2}\sqrt{2}\right) \\
\int_{0}^{1} \frac{x^2}{\sqrt{1 - x^4}} \,\mathrm{d}x &= \left\{\frac{1}{2}\sqrt{2}\,F\left[\arccos(x);\frac{1}{2}\sqrt{2}\right] - \sqrt{2}\,E\left[\arccos(x);\frac{1}{2}\sqrt{2}\right] \right\}_{x = 0}^{x = 1} \\
&= \frac{1}{2}\sqrt{2}\left[2\,E\left(\frac{1}{2}\sqrt{2}\right) - K\left(\frac{1}{2}\sqrt{2}\right)\right]
\end{align}</math>

Inserting these two integrals into the above form gives:
<math display="block">\begin{align}
&\left[\int_{0}^{1} \frac{1}{\sqrt{1 - x^4}} \,\mathrm{d}x\right]\left[\int_{0}^{1} \frac{x^2}{\sqrt{1 - x^4}} \,\mathrm{d}x\right] = \int_{0}^{1} \int_{0}^{1} \frac{x^3 (y^2+1)}{\sqrt{(1 - x^4)(1 - x^4 y^4)}} \,\mathrm{d}x\,\mathrm{d}y \\
&\hphantom{\quad\quad\quad}= \int_{0}^{1} \left\{\frac{y^2 + 1}{2\,y^2}\left[\text{artanh}\left(y^2\right) - \text{artanh}\left(\frac{\sqrt{1 - x^4}\,y^2}{\sqrt{1 - x^4 y^4}}\right)\right]\right\}_{x = 0}^{x = 1}\,\mathrm{d}y \\
&\hphantom{\quad\quad\quad}= \int_{0}^{1} \frac{y^2 + 1}{2\,y^2}\,\mathrm{artanh}\left(y^2\right) \,\mathrm{d}y \\
&\hphantom{\quad\quad\quad}= \left[\arctan(y) - \frac{1 - y^2}{2\,y}\,\mathrm{artanh}\left(y^2\right)\right]_{y = 0}^{y = 1} \\
&\hphantom{\quad\quad\quad}= \arctan(1) \\
&\hphantom{\quad\quad\quad}= \frac{\pi}{4}
\end{align}</math>

For the [[Lemniscate elliptic functions|lemniscatic]] special case of [[Legendre's relation]], this result emerges:
<math display="block">K\left(\frac{1}{2}\sqrt{2}\right)\left[2\,E\left(\frac{1}{2}\sqrt{2}\right) - K\left(\frac{1}{2}\sqrt{2}\right)\right] = \frac{\pi}{2} </math>

==See also==

* {{annotated link|Cavalieri's principle}} − an early particular case
* {{annotated link|Coarea formula}} − generalization to geometric measure theory
* [[Disintegration theorem]] – theorem in measure theory, a restricted converse to Fubini's theorem
* [[Fubini's theorem for distributions]] - version of Fubini's theorem for distributions, that is, generalized functions
* [[Kuratowski–Ulam theorem]] – analog of Fubini's theorem for arbitrary second countable Baire spaces
* [[Symmetry of second derivatives]] − analogue for differentiation
* {{annotated link|Fubini's nightmare}}

==References==
{{reflist}}

==Further reading==
*{{citation |last= Billingsley |first= Patrick |year= 1995 |authorlink= Patrick Billingsley |title= Probability and Measure |location= New York |publisher= Wiley |chapter= Product Measure and Fubini's Theorem |pages= 231–240 |isbn= 0-471-00710-2 }}
*{{citation |last= DiBenedetto |first= Emmanuele |date= 2002 |title=Real Analysis |series= Birkhäuser Advanced Texts: Basler Lehrbücher |publisher= Birkhäuser |place= Boston |isbn= 0-8176-4231-5 |doi= 10.1007/978-1-4612-0117-5 |mr= 1897317 }}
*{{citation |last= Fremlin |first= D. H. |date= 2003 |title= Measure theory |volume= 2 |publisher= Torres Fremlin |place= Colchester |isbn= 0-9538129-2-8 |mr= 2462280 }}
*{{citation |last= Weir |first= Alan J. |date= 1973 |title= Lebesgue Integration and Measure |location= Cambridge |publisher= Cambridge University Press |pages= 83–92 |chapter= Fubini's Theorem |isbn= 0-521-08728-7 }}


==External links==
* {{SpringerEOM|title=Fubini theorem |id=F/f041870 |first=L. D. |last=Kudryavtsev}}

{{Measure theory}}

{{DEFAULTSORT:Fubini's Theorem}}
[[Category:Theorems in measure theory]]
[[Category:Theorems in calculus]]
[[Category:Articles containing proofs]]