Solve non-linear non homogeneous differential equation with python (Duffing oscillator)

Question

I try to solve the Duffing equation using odeint:

def func(z, t):
    q, p = z
    return [p, F*np.cos(OMEGA * t) + Fm*np.cos(3*OMEGA * t) + 2*gamma*