Reason: Proof of the Cauchy-Schwarz and triangle inequalities via the quadratic-in-t argument, with the degenerate case handled through the null-support argument. Approved by Aaron.
Proof
Throughout we use the linearity and monotonicity of expectation from Linearity and Monotonicity of the Lebesgue Integral and the closure properties of square-integrable random variables from that definition (sums and scalar multiples are square-integrable; products are integrable; in particular all expectations below are defined and finite).
Part 1. Write a=β₯Xβ₯22β=E[X2], b=E[XY], and c=β₯Yβ₯22β=E[Y2].
Case c>0. For every real t, the random variable (X+tY)2 is nonnegative, and expanding pointwise, (X+tY)2=X2+2tXY+t2Y2, so by linearity and monotonicity of expectation
0β€E[(X+tY)2]=a+2tb+t2c.
Taking t=βb/c gives 0β€aβ2b2/c+b2/c=aβb2/c, hence b2β€ac. Since the nonnegative square root is order preserving on nonnegative reals (if 0β€uβ€v and uβ>vβ then squaring the inequality of nonnegative numbers would give u>v), we get β£bβ£=b2ββ€aβcβ=β₯Xβ₯2ββ₯Yβ₯2β, using acβ=aβcβ (both sides are nonnegative with square ac, so they agree by the uniqueness in the cited theorem) and b2β=β£bβ£ (both sides are nonnegative with square b2).
where the middle step is Part 1 together with E[XY]β€β£E[XY]β£. Both β₯X+Yβ₯2β and β₯Xβ₯2β+β₯Yβ₯2β are nonnegative, so the order-preservation of the nonnegative square root recalled in Part 1 yields β₯X+Yβ₯2ββ€β₯Xβ₯2β+β₯Yβ₯2β. β