6 / VERSION 5A 4-27-77 PT
20 /IF X<=0 - IT IS AN ERROR
95 LA0, 0 /.6118015411059928976
101 LA1, 0 /.5045560107523952859
107 LA2, 7777 /.4054651081081643810
113 LA3, 7777 /.3184537311185346147
119 LA4, 7776 /.2231435513142097553
125 LA5, 7776 /.1718502569266592214
131 LA6, 7775 /.08961215868968712374
137 LA7, 7773 /.030771658666753687
167 FLDA% BPDAL,1 /ADDRESS
171 JLE DALERR /X <= 0 IS ERROR
177 DALA, FADD DAL1 /ADD BACK
178 FSTA XRDAL /STORE AT X
179 /EXPONENT STORED IN XR0
180 /MANTISSA STORED IN XR1-5
181 /PICK UP EXP + MULTIPLY BY LOGE(2)
186 /XRDAL IS NOW FRACTION IN RANGE .5<=F<1.0
188 /LOG(F)=LOG(A(K1)*A(K2)...(F))-(LOG(A(K1))+
190 /FIT F IN A 1/16 RANGE
191 /I.E. 1/2-9/16,9/16-10/16,ETC.
192 /MULTIPLY F BY APPROPRIATE A(K) MULTIPLIER
193 /KEEP RUNNING SUM OF LOG(A(K))
203 FMUL D16 /16 REAL PARTS
206 FLDA A0,1 /GET MULTIPLIER
208 FLDA LA0,1 /ADD LOG(A(K)) TO SUM
213 /NOW F>1. USE TAYLOR SERIES
214 /LOG(T)=Z-(Z^2)/2+(Z^3)/3+... WHERE Z=T-1