Следующий код сделает это:
import scipy.stats as ss
import numpy as np
import matplotlib.pyplot as plt
center = 100
n = 1001
std = 20
x = np.arange(0, n)
prob = ss.norm.pdf(x,loc=center, scale = std )
prob = prob / prob.sum() #normalize the probabilities so their sum is 1
nums = np.random.choice(x, 20, p=prob)
nums