Sharp PC-1500 - Page 31

All and more abo ut Sharp PC-1500 at http://w ww .PC-1500.info
PROGRAM
FIRST ORDER DIFFERENTI AL
EQUATION
PROGRAM NO.
3
TI T LE P5 -A-7
[ Program List I ( Memory Contents I
10: "A" : IN PU T "X0=
1 + )112
";X ,
11
Y0=-" ; Y, "H
A
= '' ; H, " T:
11
; T B 1 -/ 1 12
20:A= i+/ . S:B=I -~.
c
5:U SING
: Q=0
25:Z =X+T:S=X
D
30 : GOSUB
500
E
40:K=HfF:R=<K-2*Q
)/2:Y = Y+R
F
J (;r, y)
50:Q=Q+3*R-K/2
G
60:X=X+H/2:GOSUB
500
H h
70:K=H* F:R= B* < K-Q
I
>: Y=Y+R
J
80:Q=Q ... 3 *R - B *K
90:GO SUB
500
K
v
J00:K=H*F:P=A *<K - Q
L
):Y=Y+R
IJ 0:Q=0+3* R-A *K
M
1 20:X =X ... H/2:GOSUB
N
500
0
130:K=H*F:R=<K-2*0
)/ 6:Y = Y+R
p
140: Q=Q+ 3 *R - K/2
0 Qn
150: l F
X< Z LET
S= X:
Yl=Y:GOTO
30
R
'"
160: IF
X= ZGOTO 200
s
Xn- 1
170:Y2=<Z-S>*<Y- Yl
Interval of
)/H+Yl
T
1s0 : BEE P
3 :LPR1NT
solutions
u x .:: .. ; 2
u
!90:LPRINT
"Y= "; Y2
:G OTO
210
v
200:BEEP
3 :LPR1NT
w
" X= .,; X
20S:LPR 1N T
"V= "; y
x
Xn
2J0:Z= Z+T :S=X:Yl=Y
y
Yn
2J5:WAIT
:PRINT
x.
z
v
y
2 20: GOTO
30
Y1
Yn-1
500: F=-X*Y
Y2
v.,
5J0:RETURN
STATUS I
474
.
Do not sale this PDF !!!
- 25 -