Closed jblake95 closed 5 years ago
Have now fixed this by using averages (between the new and previous values) of the Lagrange coefficients, as in Curtis, 'Orbital Mechanics for Engineering Students'. The example case converges upon the fourth iteration, as expected.
Orbit improvement algorithm for the Gauss method is not working. The iterations lead to a divergence in rho_i instead of a convergence. Need to debug this.