How to generate random points inside the cylinder?

I have a cylinder with the known radius r and the known height h and with the known coordinates of the centers – (X1, Y1, Z1) and (X2,Y2,Z2). How can I get random points inside the cylinder? I'm writing my code in Fortran.