Evolution of the orbital eccentricities e_b and e_c and apsidal alignment angle Δϖ = ϖ_b - ϖ_c. Red curves represent the evolution obtained from the Wisdom-Holman N-body integration evolution of the best fit (see Fig. 5), while red curves are the Wisdom-Holman N-body integration with GR precession included. Gray dashed and dotted curves represent the evolution of the best fit but derived from the octupole secular perturbation theory and the octupole evolution with a GR precession correction, respectively.