Figure 5.
These graphs show the evolution of the orbital parameter of model M3. The stronger mutual interaction between the two stars results in strong mode of VHS mechanism, resulting in a constant ΔΩ = 0°, which is a complete overlap for integration of both secular and equations of motion. There are tiny inclination oscillations for both integration of equations of motion (darker) and secular equations (lighter), which have slightly different amplitude and period.