Предположим, что вы сохранили исходный cij в массиве numpy, вы можете сначала суммировать такие термины, как c11 + c21 + c31 + c41.Это может быть сделано путем суммирования каждого столбца, попробуйте c.sum (axis = 0)
>>> import numpy as np
>>> c = np.arange(24).reshape(4,6)
>>> c
array([[ 0, 1, 2, 3, 4, 5],
[ 6, 7, 8, 9, 10, 11],
[12, 13, 14, 15, 16, 17],
[18, 19, 20, 21, 22, 23]])
>>> c = c.sum(axis=0)
>>> c
array([36, 40, 44, 48, 52, 56])