Consider the set of all partitions {P} of [a,b], defined as
{P:P=\{a=x_0<x_1<x_2<...<x_{n-1}<x_n}
The x_ks divide your interval into a countable number of sub-intervals \left[ x_{k-1}, x_k \right]. The definition of the lower/upper integral is rather abstract and useless - it defines them as the supremum/infimum of L(f,P) and U(f,P) respectively on the set of
all possible partitions P. But there exists a nice lemma that if you take a sequence of partitions (P_k), and the limit of their partition sizes \delta = x_{k}-x_{k-1} goes to zero (I think that's the requirement, it might be more general), then the supremum (infimum) of the L(f,P) or U(f,P) is the same as the limit of the sequence U(f,p_k):
sup{L(f,P):P is a partition of [a,b]}=lim U(f,p_k)
Intuitively, as your partition gets "finer and finer", the error of the lower (upper) sum gets less. I forgot how the rigorous version goes here. The point is, instead of working with supremums, you work with a limits of sequences (which are 'nicer').
but i don't understand how how you could make up a sequence of partitions, to get this
U(f)-L(f) \le U(f,P_k) - L(f,P_k) = \sum_{j=1}^n (2 x_{j-1}(x_j-x_{j-1}) + (x_j-x_{j-1})^2)(x_j-x_{j-1})
I'm not entirely sure what you're asking about - let's see, the \sum_{j=1}^n (2 x_{j-1}(x_j-x_{j-1}) + (x_j-x_{j-1})^2)(x_j-x_{j-1})
is unique to f(x)=x^2, under any partition (I derived at the top of my first post, assuming x>0 and that x^2 is strictly increasing; you get something very similar when x<0).
To actually work with this, it would be nice to actually choose a sequence of partitions (that gets succesively finer): an easy one is
P_k=P \{ a=x_0<x_1=a+\frac{b-a}{k}<...<x_j=a+j\frac{b-a}{k}<...<x_k=a+k\frac{b-a}{k}=b \}
Which reduces down to what
mathwonk was using: (I'll quote his post)
for x^2 on [0,a] you have a subdivision of length a/n, and the upper value on [(k-1)a/n, k/n] is (ka)^2/n^2, while the lower value is [(k-1)a]^2/n^2.
so the lower sum is (a/n)[0^2/n^2 + a^2/n^2 +...+[(n-1)a]^2/n^2],
and the upper sum is (a/n)[a^2/n^2 + ...+[(n-1)a]^2/n^2 + [(n)a]^2/n^2
i.e.,
\begin{align*}<br />
lim \left(U(f)-L(f) \right)&=lim\left( U(f,P_k)-L(f,P_k) \right) \\<br />
&=lim \sum_{j=1}^n (2 x_{j-1}(x_j-x_{j-1}) + (x_j-x_{j-1})^2)(x_j-x_{j-1}) \\<br />
&=lim \sum_{j=1}^{k}\left( 2x_j \left( \frac{b-a}{k} \right) + \left( \frac{b-a}{k} \right) ^2 \right) \left( \frac{b-a}{k} \right) \\<br />
&=lim_{k \rightarrow \infty } \frac{(b-a)^3}{k} \\<br />
&=0 \end{align}