Integral of Power/Fermat's Proof

Theorem
$\displaystyle \forall n \in \Q_{>0}: \int_0^b x^n \mathrm d x = \frac {b^{n+1}} {n+1}$

Proof
First let $n$ be a positive integer.

Take a real number $r \in \R$ such that $0 < r < 1$ but reasonably close to $1$.

Consider a subdivision $S$ of the closed interval $\left[{0 \,.\,.\, b}\right]$ defined as:
 * $S = \left\{{0, \ldots, r^2 b, r b, b}\right\}$

... that is, by taking as the points of subdivision successive powers of $r$.

Now we take the upper sum $U \left({S}\right)$ over $S$ (starting from the right):

Now we let $r \to 1$ and see that each of the terms on the bottom also approach $1$.

Thus:
 * $\displaystyle \lim_{r \to 1} S = \frac {b^{n+1}}{n+1}$

That is:
 * $\displaystyle \int_0^b x^n \mathrm d x = \frac {b^{n+1}} {n+1}$

for every positive integer $n$.

Now assume $n = \dfrac p q$ be a strictly positive rational number.

We set $s = r^{1/q}$ and proceed:

As $r \to 1$ we have $s \to 1$ and so that last expression shows:

So the expression for the main result still holds for rational $n$.

Historical Note
This method was used by Fermat, and predated the work done by Newton and Leibniz by a considerable period.