def METHOD_NAME(freqs, decays, amplitudes, phases, time_step):
    t = np.arange(NTIME)
    complex_amplitudes = amplitudes * np.exp(1j * phases)
    complex_freqs = 2 * np.pi * freqs - 1j * decays
    signal = np.zeros(len(t), dtype=complex)
    for i in range(len(freqs)):
        signal += complex_amplitudes[i] * np.exp(-1j * complex_freqs[i] * t * time_step)
    return signal