First i looked at the case of ## \epsilon = 0## and came to the conclusion, that this oscillator has a circular limit cycle in a phase space trajectory, when plotted with the axes x and ##\dot{x}##.
I have found that ##x_0^p (t) = a_1 \cos(t)## which implies that all other Fourier- coefficients...
How did you find PF: I am studying physics and looked for a place where i can ask questions and talk about interesting new developments, some quick googling got me here.