rigidbody¶

probnum.problems.zoo.diffeq.rigidbody(t0=0.0, tmax=20.0, y0=None, params=(- 2.0, 1.25, - 0.5))[source]

Initial value problem (IVP) for rigid body dynamics without external forces

The rigid body dynamics without external forces is defined through

$\begin{split}f(t, y) = \begin{pmatrix} y_2 y_3 \\ -y_1 y_3 \\ -0.51 \cdot y_1 y_2 \end{pmatrix}\end{split}$

The ODE system has no parameters. This implementation includes the Jacobian $$J_f$$ of $$f$$.

Parameters
• t0 – Initial time. Default is 0.0

• tmax – Final time. Default is 20.0

• y0(shape=(3, )) – Initial value. Default is [1., 0., 0.9].

• params – Parameter of the rigid body problem. Default is (-2.0, 1.25, -0.5).

Returns

InitialValueProblem object describing the rigid body dynamics IVP with the prescribed configuration.

Return type

InitialValueProblem