¿Qué función cúbica puedo utilizar para obtener una distribución aleatoria sesgada hacia el centro?

Quiero un círculo con más círculos dentro (no necesariamente adentro). La posición de los círculos interiores se determina al azar de tal manera que debería haber más círculos en el centro y menos y menos a medida que se aleja del borde del círculo.

De esta pregunta, deduje que los números pueden estar sesgados usando f(x) en lugar de solo x, x siendo el número aleatorio, por supuesto. Aquí está el código:

def distBiasedMiddle(_min, _max):
r = lambda : random.uniform(0, 1)
r1, r2 = r(), r()
bias = lambda _r: _r**3 * (1 + _max – _min) + _min
return (bias(r1), bias(r2))

testBias = distBiasedMiddle

rps = []
for i in range(100):
rps.append(testBias(-50, 50))

while True:

for p in rps:
     pygame.draw.circle(window, (255, 0, 0), (p[0] + 300, p[1] + 300), 2)
pygame.draw.circle(window, (0,255,0), (300, 300), 50, 1)

Para las personas no familiarizadas con pygame, el círculo verde tiene un radio de 50 y se dibuja en (300, 300) y los círculos rojos se dibujan en (p[0] + 300, p[1] + 300)

Aquí usé x ** 3 pero no funcionó. Necesito ayuda para elegir la función correcta. Aquí está lo que tengo.

Sé que el resto del código es correcto porque utilicé otra función mencionada en esa respuesta (f(X)=2X/(1+X)) y da la distribución correcta, sesgada hacia los valores más altos:

def distBiasedUpper(_min, _max):
r = lambda : random.uniform(0, 1)
r1, r2 = r(), r()
bias = lambda _r: (2 * _r / (1 + _r)) * (1 + _max – _min) + _min
return (bias(r1), bias(r2))