def METHOD_NAME(dt, eps, df, Xi, kappa):
    '''
    Simulates the trajectory of a cavity driven by a measurement pulse
    dt = time step
    eps = array with drive pulse amplitude
    df is detuning of drive relative to cavity resonance
    '''
    alpha_m = np.zeros(len(eps)+1)+0.j
    alpha_p = np.zeros(len(eps)+1)+0.j
    #print alpha_p
    k=1
    print(2.*np.pi*df+2.*np.pi*Xi)
    print(2.*np.pi*df-2.*np.pi*Xi)
    for epst in eps:
        #print epst
        alpha_p[k] = (-1.j*epst -
                        (1.j*(2.*np.pi*df+2.*np.pi*Xi) + 2.*np.pi*kappa/2.)*alpha_p[k-1]
                        )*dt+alpha_p[k-1]
        alpha_m[k] = (-1.j*epst -
                        (1.j*(2.*np.pi*df-2.*np.pi*Xi) + 2.*np.pi*kappa/2)*alpha_m[k-1]
                        )*dt+alpha_m[k-1]
        k+=1
    return [alpha_p,alpha_m]