Numerical solution of second-order ode, Numerical solution of second-order ode ,16-61 – HP 50g Graphing Calculator User Manual

Page 538

Advertising
background image

Page 16-61

LL@)PICT To recover menu and return to PICT environment.
@(X,Y)@

To determine coordinates of any point on the graph.

Use the

š™ keys to move the cursor around the plot area. At the bottom of

the screen you will see the coordinates of the cursor as (X,Y), i.e., the calculator
uses X and Y as the default names for the horizontal and vertical axes,
respectively. Press

L@CANCL to recover the menu and return to the PLOT

WINDOW environment. Finally, press

$ to return to normal display.

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, x" = - 18.75 x - 1.962 x',

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' ]

T

, and A is the 2 x 2 matrix

shown below.

The initial conditions are now written as w = [0 6]

T

, for t = 0. (Note: The

symbol [ ]

T

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

dt

dx

x

dt

x

d

=

962

.

1

75

.

18

2

2

=

'

962

.

1

75

.

18

1

0

'

'

x

x

x

x

Advertising