Как получить кумулятивную функцию распределения с помощью NumPy?

Я хочу создать CDF с NumPy , мой код следующий:

histo = np.zeros(4096, dtype = np.int32)
for x in range(0, width):
   for y in range(0, height):
      histo[data[x][y]] += 1
      q = 0 
   cdf = list()
   for i in histo:
      q = q + i
      cdf.append(q)

Я иду по массиву, но выполнение программы занимает много времени. Есть встроенная функция с этой функцией, не так ли?