The self-starting Radau IIA integrator is at least comparable with the best of other integrators in speed and accuracy, and it is often superior, particularly at high accuracies. Here, integration of the normalized two-body problem from t0 = 0 [s] to t = 3600 [s] for an eccentricity of e = 0.1 is implemented.
Meysam Mahooti (2021). Radau IIA (https://www.mathworks.com/matlabcentral/fileexchange/56162-radau-iia), MATLAB Central File Exchange. Retrieved .
Just solved an index-1 DAE problem with singular mass matrix, and obtained the same solution as found using with the ode15s integrator. Really useful code.
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!