Upon integrating Eqns. (39-40) numerically, one discovers that the system oscillates very quickly, on the order of one or more oscillations per unit time scale. This is bad for computations, because very small time steps will be needed to accurately resolve the behavior of the system. Secondly, the important dynamics of the system will be happening ``on top of'' these oscillations.
One solution to this problem is to essentially cancel the change in phase of the system. A change of variables can be made in order to do this. The effects of b in our analysis will be neglected with the assumption that the added injection will do little to change the frequency at which the system oscillates.
The original e in Eqn. (1) will be transformed into a new variable E by
where exp represents the exponential function. Substituting this into Eqn. (1), one obtains an expression for the rate of change of the new variable E of the form
Essentially, what has been done is to add a term to Eqn. (1) which will allow one to ``follow'' the system around through its changes in phase. The new variable E will now be called e.
Now a value must be found for . Essentially, it is the frequency at which the system would like to naturally oscillate. The imaginary terms in Eqn. (42) are
So, in order to suppress the imaginary part of Eqn (42), it is neccesary to set
Though n varies with time, the value can be substituted for n in Eqn. (43) where is the value of n at the fixed point of the system. The assumption made here is that the expression for the frequency thus obtained will still hold for values of n away from the fixed point. This argument is reinforced by the fact that oscillations of n, if any, tend to be around this fixed point and have small amplitude (see Figure [1]).
Setting Eqn. (39) equal to zero and ignoring , can be found to be
This allows Eqn. (39) to be written