De Moivre's theorem states for integer . For , expanding the cube and taking its real part gives
Rearranging gives the required identity with
No division by is needed, so the proof also covers its zeros.
Dividing Newton's equation by the mass gives the stated model: is the acceleration, represents viscous damping, is the spring's restoring acceleration, and the right side is the applied driving force per unit mass. The complex notation is a convenience; the physical displacement is the real part. Put . The homogeneous characteristic roots are
For and damping , oscillations occur when , the underdamped case of the damped harmonic oscillator.
For a steady response , substitution gives
Taking the modulus and phase gives
The phase is between zero and for positive damping and frequency; an unqualified arctangent loses the correct quadrant above resonance.
With the real drive and real response , the average driving power equals the dissipated viscous power:
For , write this as
It is largest when the squared term vanishes, giving
This power peak is at even when the displacement-amplitude peak is shifted. The frequency bandwidth of viscous oscillator power absorption follows from the half-power condition :
In fact the bandwidth identity is exact in this model. Sharp resonance additionally gives . With zero damping at exact resonance, there is no bounded steady harmonic response to which these finite-power formulas apply.
Figure 1.
Driven-oscillator power peaks at the natural frequency; the exact separation of the half-power frequencies is the damping coefficient
. Details.
One scalar initial condition is needed in each case, because each equation is first order and each regular solution has one integration constant. The usual Picard-Lindelöf theorem applies locally at a regular initial point. For part (c), this requires and ; a condition at an excluded singular endpoint cannot specify a classical solution there. Arbitrary distant boundary data also need not be compatible with one common interval of existence.
Here , so its even powers alternate between and its odd powers between . Collecting even and odd terms of the matrix exponential yields
Thus , , and . It is the rotation matrix for a counterclockwise rotation through angle , preserving lengths and orientation.

There are unlisted articles, also show them or only show them.