The equation of motion is Equation 4.19 in Thornton and Marion, with the frequency omega set to unity. Notice that the limit cycle depends on mu and a but is independent of the initial conditions on position and velocity.
By Robert Johnson, using the Physlet Animator and dataGraph classes from Wolfgang Christian and Mario Belloni.