HP 15c - Manual (Page 77)

Section 3: Calculating in Complex Mode 77
The program listed below is set up to compute the values of y
k
from evenly spaced values of
x
k
. You must provide a subroutine labeled "B" that places Im(P(z)) in the real X-register. The
program uses inputs that specify the step size h, the number of points n along the real axis,
and z
0
= x
0
+ iy
0
, the initial point on the streamline. You must enter n, h, and z
0
into the Z-,
Y-, and X-registers before running the program.
The program computes the values of z
k
and stores them in matrix A in the form a
k1
= x
k-1
and
a
k2
= y
k-1
for k = 1, 2, ... , n.
Keystrokes Display
|¥ Program mode.
´CLEARM 000-
´bA 001-42,21,11
) 002- 33
O4 003- 44 4 Stores h in R
4
.
) 004- 33
2 005- 2
´mA 006-42,23,11 Dimensions matrix A to be n×2.
|` 007- 43 35
O>A 008-44,16,11 Makes all elements of A be
zero.
OV 009- 44 25 Stores zero in Index register.
´>1 010-42,16, 1 Sets R
0
= R
1
= 1.
|( 011- 43 33 Recalls z
0
to X-registers.
O2 012- 44 2 Stores x
0
in R
2
.
´UOA
´U 013u 44 11 Sets a
11
= x
0
.
´} 014- 42 30
O3 015- 44 3 Stores y
0
in R
3
.
´UOA 016u 44 11 Sets a
12
= y
0
.
´Ut1 017- 22 1 Brances if matrix A not
full (n > 1).
´b0 018-42,21, 0
l>A 019-45,16,11 Recalls descriptor of matrix A.
|n 020- 43 32
´b1 021-42,21, 1
´} 022- 42 30 Restores z
0
.
GB 023- 32 12 Calculates Im(P(z
0
)) (or
Re(P(z
0
)) for equipotential line.)
O5 024- 44 5 Stores c in R
5
.
´b2 025-42,21, 2 Loop for finding y
k
.
1 026- 1