4 / AND B IS DOUBLE PRECISION
6 / VERSION 5A 4-26-77 MH
11 FSTA SIGN /SAVE SIGN OF EXPONENT
12 JGE POSINT /ITS POSITIVE
19 FLDA ONE /START WITH ONE
23 LOOP, JEQ YES /DONE IF ITS ZERO
24 FDIV TWO /DIVIDE BY TWO
28 FMUL TWO /IS EXPONENT ODD ?
35 FLDA TEMP /EXPONENT OVER 2
37 JXN LOOP,1+ /LOOP IF MORE BITS
38 YES, FLDA XR1 /DONE, RESTORE XR 1
40 FLDA SIGN /CHECK SIGN OF EXPONENT
41 JLT INVERT /IT WS NEGATIVE, INVERT RESULT
43 FLDA PROD /RETURN ANSWER
46 FLDA ONE /RETURN WITH 1/PROD
51 JA SQUARE /GO SQUARE THE BASE