def du_dx(x, y):
    wa = 1
    wp = 0.6
    g = 0.6
    n = 1
    A = ((n * wp) + (wa / 2))
    B = (((1 + n) * wp) - (wa / 2))
    X = (n + 1)
    C = np.sqrt(X)
    dydx_1 = ((A * y[1]) + ((g * C) * y[3]))
    dydx_2 = (((- A) * y[0]) - ((g * C) * y[2]))
    dydx_3 = ((B * y[3]) + ((g * C) * y[1]))
    dydx_4 = (((- B) * y[2]) - ((g * C) * y[0]))
    return [dydx_1, dydx_2, dydx_3, dydx_4]
def du_dx(x, y):
    wa = 1
    wp = 0.6
    g = 0.6
    n = 1
    A = ((n * wp) + (wa / 2))
    B = (((1 + n) * wp) - (wa / 2))
    X = (n + 1)
    C = np.sqrt(X)
    dydx_1 = ((A * y[1]) + ((g * C) * y[3]))
    dydx_2 = (((- A) * y[0]) - ((g * C) * y[2]))
    dydx_3 = ((B * y[3]) + ((g * C) * y[1]))
    dydx_4 = (((- B) * y[2]) - ((g * C) * y[0]))
    return [dydx_1, dydx_2, dydx_3, dydx_4]