Not sure if that's exactly what you want, but it should give you some distribution with more points per volume in the center of the sphere:
You can use this simple formula to find a random evenly distributed point on a spheres surfaces: Generate a random sample of points distributed on the surface of a unit sphere
When you take this point and move it inwards by multiplying it with a random value and dividing by r:
(x,y,z) *= randomValue
where (x,y,z) is the vector pointing at the random point on the surface and randomValue is some random value between 0 and 1.
This gives an even distribution with respect to r, but since the surface are gets smaller with smaller radius, your points will concentrate in the middle.
If you want to further tweak the distribution, you can add some power to randomValue:
(x,y,z) *= pow(randomValue, exponent)
where the higher the exponent, the more points you will get in the middle of the sphere.