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