File RTD

Directory of image this file is from
This file as a plain text file

[RTD
[PROGRAM TO FIND SQUARE ROOT OF X AND ITS RECIPROCAL
[ENTRY RECIPROCAL OF SQUARE ROOT
[JMP 1K
[X
[RESULT
#1K LDA
0
JMP 1X-1
ADD 1R+1
STC p+111
LDA 2
[TAKE ABSOLUTE VALUE
[CHECK FOR ZERO DIVIDE
STC p+110
LDA i2
AZE i
JMP 1H [ZERO DIVIDE
APO
COM
STC p+102+1
LDA i2
APO
COM
STC p+76+2
[CALCULATE RESULT EXPONENT
ADD p+75
SCR i1
APO
JMP p+3
LZE
ADD 1Q+1 [CONSTANT 1
COM
STA 3
ROL i1
APO i
COM
SCR i1
STC p+61
[MODIFY FRACTION IF NECESSARY
LZE
ADD 1N-6 [=-1
STC p+56
ADD 3
STC p+41
[APPROXIMATE RESULT
JMP 1M
p+52
p+70
p+53
JMP 1A
p+70
p+50
p+47
[A=X**2
[LOOP
 JMP 1M
p+45
p+44
p+51
[A=A*N
JMP 1M
p+47
p+35
p+45
[A=3-A
JMP 1S
p+40
p+42
p+41
[X=A*X
JMP 1M
p+31
p+36
p+27
[HALVE RESULT
LDA
p+25
ADD 1N-6 [=-1
STC p+23
SRO i
3567
JMP p-26
[DONE WITH LOOP
SET i3
 0
LDA 3
ADD p+14
STA 3
LDA
p+11+1
STA i3
LDA
p+6+2
STA i3
 JMP RETURN
  0
0
0
880
0
0
 2
3000
0000
  0
0
0
  0
0
0
[CONSTANTS FOR APPROXIMATION OF
[  1|ROOT X
 0001
5440
4060
 0002
2103
3727
[SQUARE ROOT OF N,
[FOUND BY MULTIPLYING
[1|SQUARE ROOT OF N  BY N
[ENTRY:
[JMP 1T
[X
[RESULT
#1T LDA
0
JMP 1X-1
ADD 1R+1
STC p+32
ADD 3
STA
p+22
STA
p+23
STC p+23
LDA 2
STC p-25
LDA i2
AZE
JMP p+5
[STORE ZERO IF ZERO
STA 3
STA i3
STA i3
JMP p+13
STC p-35+1
LDA i2
STC p-37+2
[X=1|SQUARE ROOT OF N
JMP 1K
p-41
 ADDR RESULT
[X=RESULT*N
JMP 1M
p-44
 ADDR RESULT
 ADDR RESULT
 JMP RETURN
[END RTD
[MJM 4-1-70



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search