Wow... somebody tell me where I go wrong.
Let a be randomly chosen in the interval [0, 1] and represent r^2. Then let b be chosen in the interval [0, 2PI] (or [0, PI] for that matter, equivalently) be the angle. Then (sqrt(a),b) is a uniform distribution for points on the unit disk.
You have to such pairs (a, b). Call them (a, b) and (a', b'). These are two points selected at random. Rotate both vectors until b = 0 (no problem, right? frame of reference...) So now you have (a, 0) and (a', b'').
The distance between these points can be calculated using the law of cosines.
D^2 = a + a' - 2sqrt(a*a')cos(b'').
Note that a = r^2, a' = r'^2, and b is theta, so this should be true.
Now, try to find the expectation value of D^2.
<D^2> = <a + a' - 2sqrt(a*a')cos(b'')>
<D^2> = <a> + <a'> - 2<sqrt(a*a')cos(b'')>
<D^2> = 1/2 + 1/2 - 2<sqrt(a*a')cos(b'')>
<D^2> = 1 - 2<sqrt(a*a')cos(b'')>
Since for any a and a', b'' can be anything between -PI and PI, cos(b'') has an equal chance of being positive or negative... so
<D^2> = 1 - 2<sqrt(a*a')cos(b'')> = 1 - 0 = 1.
So <D^2> = 1.
Is the average squared distance 1? I know that doesn't really relate to the problem directly, but I thought it was interesting.