Given a probability density, how can I sample from the induced distribution? Let $f$ be an integratable function such that $\int f(x) dx=1$. If we want to take random samples from this, using whatever programming language one pleases, we should compute $F(t)=\int_{-\infty}^t f(x) dx$, invert this function and feed it numbers drawn from a uniform distribution on $[0,1]$.
However I now want to sample coordinates $(u,v)\in \mathbb{R}^2 \backslash \lbrace 0 \rbrace$ such that the probability of $(u,v)$ lying in a set $E$ is given by
$$
\int_E (u^2+v^2)^{-\frac{3}{2}}du dv
$$
I don't see how I can now try to find
$$F(s,t)=\int_{-\infty}^s\int_{-\infty}^t (u^2+v^2)^{-\frac32} du dv$$
as I get into troubles close to the zero.
What other way is there to obtain a sample following this distribution?
A note for the context: If we consider all lines in the Euclidean plane with Cartesian coordinates, not passing through the origin, they can be represented via $ux+vx=1$. If we impose the condition that the probability densitiy should be invariant under Euclidean transformations, then we arrive at the above distribution. See here
 A: Your probability density on $u,v$ is not a valid probability density.
A probability density should have $\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f(u,v)\mathrm{d}u\mathrm{d}v=1$
If you convert to polar coordinates you will see that
$\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}(u^2+v^2)^{-3/2}\mathrm{d}u\mathrm{d}v=\int_0^{\infty} 2\pi r^{-2} \mathrm{d}r=\infty$
Removing the point $\{0\}$ from $\mathbb{R}^2$ won't help either - you have to remove a neighborhood of zero in order to make your density normalisable.
A: The book you cite describes
the Haar measure on the space of lines in the plane, with a view to making it the intensity measure for a Poisson process.  In the pages before the bit you cite it  mentions  the fact that an invariant measure might have infinite mass, and that if you want a probability distribution come from it you have to condition by restricting yourself to sets of finite measure.  It then goes on to describe the polar coordinate representation that makes this process practical, described in Angela's answer.
An intuition behind the non-finiteness of the integral is: the infinitely overwhelming  majority of lines in the plane are very far from the origin. Such a line is described by an equation $ux+vy=1$ for tiny $u$ and $v$.
