Numerical Solution Of Second-Order Ode - HP 49g+ User Manual

Numerical solution of second-order ODE

Integration of second-order ODEs can be accomplished by defining the
solution as a vector. As an example, suppose that a spring-mass system is
subject to a damping force proportional to its speed, so that the resulting
differential equation is:
or,
subject to the initial conditions, v = x' = 6, x = 0, at t = 0. We want to find x,
x' at t = 2.
Re-write the ODE as: w' = Aw, where w = [ x x' ]
matrix shown below.
x
The initial conditions are now written as w = [0 6]
T
symbol [ ] means the transpose of the vector or matrix).
To solve this problem, first, create and store the matrix A, e.g., in ALG mode:
Then, activate the numerical differential equation solver by using: ‚ Ï
˜ @@@OK@@@ . To solve the differential equation with starting time t = 0 and
final time t = 2, the input form for the differential equation solver should look
as follows (notice that the Init: value for the Soln: is a vector [0, 6]):
2
d x
18 . 75 . 1
x
2
dt
x" = - 18.75 x - 1.962 x',
'
x 0 1
' 18 . 75 . 1 962
dx
962
dt
T
, and A is the 2 x 2
x
x '
T
, for t = 0. (Note: The
Page 16-64