Gauss-Jackson Eighth-Order ODE Solver (Fixed Step Size)
Description:
The Gauss Jackson Eighth Order fixed-step ODE solver is finally available!
You have probably noticed that is next to impossible to find source code for this solver.
The only code publicly available is FORTRAN 90 (NETLIB), which is fine if you want to write your own force models in that code language. If you need to call your force models in MATLAB, JAVA or C; FORTRAN code won't work well without an overhaul.
I wrote this function generically enough to allow for additional applications aside from its obvious choice as a numerical orbit propagator.
Please read the comments in the code if you would like to learn how to call this routine. Also look at the demo function I have included for a template on how to write the first and second order ODE calling functions.
This function has been tested successfully on the two body force equation, and the JAT force model (available in the NASA Orbit Determination Toolbox (ODTBX)).
See GJ8_Demo.m for sample orbit propagation comparisons using GJ8 vs ODE113.
引用格式
Darin Koblick (2024). Gauss-Jackson Eighth-Order ODE Solver (Fixed Step Size) (https://www.mathworks.com/matlabcentral/fileexchange/36542-gauss-jackson-eighth-order-ode-solver-fixed-step-size), MATLAB Central File Exchange. 检索来源 .
MATLAB 版本兼容性
平台兼容性
Windows macOS Linux类别
- MATLAB > Mathematics > Numerical Integration and Differential Equations > Ordinary Differential Equations >
标签
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!GJ8/
版本 | 已发布 | 发行说明 | |
---|---|---|---|
1.0.0.0 |