Sin pérdida de generalidad, vamos a tratar sólo con las integrales con un límite inferior de $a = 0$ (teniendo en cuenta que, en general, tenemos que $\int_a^b = \int_0^b - \int_0^a$). Por lo tanto, debemos mostrar que $\int_0^b x^n = \frac{b^{n+1}}{n+1}$
Deje $\mathcal{P}_m = \{t_0,t_1,\dots,t_m\} = \{0, b\cdot\frac{1}{m}, b\cdot\frac{2}{m}, \dots, b\cdot\frac{m-1}{m},b\}$ (al calcular las integrales de Riemann desde cero, este debe ser su partición). Observar que $t_{i} - t_{i-1} = b\cdot\frac{1}{m}$ y $M_{i} = \big(b \cdot \frac{i}{m}\big)^n = (\frac{b}{m})^n \cdot i^n$ (vamos a tratar con la parte superior de sumas). En consecuencia,
$U(\mathcal{P}_m, x^n) = \displaystyle\sum_{i=1}^m M_i(t_i - t_{i-1}) = \displaystyle\sum_{i=1}^m (\frac{b}{m})^n \cdot i^n \frac{b}{m} = (\frac{b}{m})^{n+1}\displaystyle\sum_{i=1}^m i^n$
Hagamos una pausa. Necesitamos lidiar con la suma de $\displaystyle\sum_{i=1}^m i^n$ de alguna manera. Me voy a referir a usted aquí. Ahora,
$
U(\mathcal{P}_m, x^n) = (\frac{b}{m})^{n+1}\displaystyle\sum_{i=1}^m i^n = (\frac{b}{m})^{n+1}\bigg(\frac{1}{n+1}\displaystyle\sum_{k=0}^n (-1)^k \binom{n+1}{k}B_k m^{n+1-k} \bigg)
$
que se simplifica a
$ \frac{b^{n+1}}{n+1} \bigg(B_0 + \text{terms involving negative powers of m} \bigg)$
Tomando nota de que $B_0 = 1$, tenemos el límite de $\displaystyle\lim_{m \to \infty} U(\mathcal{P}_m, x^n) = \frac{b^{n+1}}{n+1}$. Espero que usted estaría de acuerdo en que estamos haciendo ahora.