View Source Integrator.SampleEqns (Integrator v0.1.2)
Functions to be used in testing
Summary
Functions
The Euler equations of a rigid body without external forces. This is a standard test problem proposed by Krogh for solvers intended for nonstiff problems [see below]. Based on "rigidode.m" from Matlab/Octave. The analytical solutions are Jacobian elliptic functions.
Simulates a point mass or particle falling through pass affected by gravity. Used for comparisons
with the Matlab/Octave ballode.m
routine
From octave b
Functions
The Euler equations of a rigid body without external forces. This is a standard test problem proposed by Krogh for solvers intended for nonstiff problems [see below]. Based on "rigidode.m" from Matlab/Octave. The analytical solutions are Jacobian elliptic functions.
See rigidode.m
in Matlab.
Shampine, L. F., and M. K. Gordon, Computer Solution of Ordinary Differential Equations, W.H. Freeman & Co., 1975
Simulates a point mass or particle falling through pass affected by gravity. Used for comparisons
with the Matlab/Octave ballode.m
routine
From octave b
fvdp = @(t,x) [x(2); (1 - x(1)^2) * x(2) - x(1)];