def G(x, n):
    return np.sin(n * np.pi * x) * g(x)

def do_integral(N):
    res = np.zeros(N)
err = np.zeros(N)
for n in range(1, N):
    res[n], err[n] = integ.quad(G, 0, 1, args=(n,))
return res, err
c, err = do_integral(500)
S = np.cumsum(c[1:] ** 2)
Написание бесконечной суммы (функции, имеющей интеграл) в Python