phi_of_t = []
for ind, val in enumerate(f_of_t):
    phi_of_t += [2 * np.pi * np.cumsum(f_of_t[:ind])]