This is an adapted Runge Kutta method (theoretical background e.g. in Hairer, Lubich, Roche "The numerical solution of differential-algebraic systems by runge-kutta methods") to solve DAEs. Already implemented is the Radau II A method of stage 1, 2 and 3 but basically every Butcher tableau can be implemented (see code for details). A simple example (the mathematical pendulum) is provided to show the usage. The nonlinear system solver is the newton method, but it can also be exchanged.
Stefan Schiessl (2021). Runge Kutta DAE solver (https://www.mathworks.com/matlabcentral/fileexchange/39857-runge-kutta-dae-solver), MATLAB Central File Exchange. Retrieved .
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!