Positions sampling implemented

For all the three components - halo, disk and bulge. The script
works, but the initial conditions it genererates has all the
velocities set to zero.