.TITLE ALOG   COMPUTE NATURAL LOGARITHM FOR REAL ARGUMENT
/
/COPYRIGHT (C) 1975
/DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASS.
/
/THIS SOFTWARE IS FURNISHED UNDER A LICENSE FOR USE ONLY
/ON A SINGLE COMPUTER SYSTEM AND MAY BE COPIED ONLY WITH
/THE INCLUSION OF THE ABOVE COPYRIGHT NOTICE.  THIS
/SOFTWARE, OR ANY OTHER COPIES THEREOF, MAY NOT BE PRO-
/VIDED OR OTHERWISE MADE AVAILABLE TO ANY OTHER PERSON
/EXCEPT FOR USE ON SUCH SYSTEM AND TO ONE WHO AGREES TO
/THESE LICENSE TERMS.  TITLE TO AND OWNERSHIP OF THE
/SOFTWARE SHALL AT ALL TIMES REMAIN IN DEC.
/
/THE INFORMATION IN THIS DOCUMENT IS SUBJECT TO CHANGE
/WITHOUT NOTICE AND SHOULD NOT BE CONSTRUED AS A COM-
/MITMENT BY DIGITAL EQUIPMENT CORPORATION.
/
/DEC ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY
/OF ITS SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DEC.
	.EJECT
/COPYRIGHT 1970,1971, DIGITAL EQUIPMENT CORP., MAYNARD, MASS. 10754
/FOR HARDWARE FLOATING POINT PROCESSOR (FPP), DEFINE 90 FPP,
/EDIT 2   1-26-71
/ EDIT 3	28-FEB-75	R.K. HYATT	INSERTED CHECK FOR ARG OF 1
/ EDIT 4	20-AUG-75	M. HEBENSTREIT	DISCLAIMER
       .GLOBL ALOG,.DA,.EE,.CAT1R
	.IFUND %FPP
	.GLOBL .AG,.AK
	.ENDC
	.IFUND %FPP
	.DEFIN FLD%,A		/FLOATING LOAD (.AG)
	JMS* A
	.ENDM
	.ENDC
	.IFDEF %FPP
FLD=713050
	.DEFIN FLD%
	FLD
	.ENDM
	.ENDC
	.IFUND %FPP
	.DEFIN FMP%,A
	JMS* A			/FLOATING MULTIPLY (.AK)
	.ENDM
	.ENDC
	.IFDEF %FPP
FMP=711440
	.DEFIN FMP%
	FMP
	.ENDM
	.ENDC
/             CALLING SEQUENCE
/      JMS*   (ALOG)          SUBR CALL
/      JMP    .+2             SUBR RETURN
/      CAL/XCTADDR            ADDR OF ARGUMENT (XCT IF INDIRECT)
/
ALOG   CAL    0              /ENTRY EXIT
       JMS*   .DA            /GET ARGUMENT (.DA)
       JMP    ALOG02
ALOG01 CAL    0              /ARGUMENT
ALOG02	LAC	ALOG01	/SEE IF ARG = 1 **003**
	JMS*	.CAT1R	/**003**
	SKP			/ARG DOESN'T =1 SO COMPUTE LOG
	JMP*	ALOG		/RETURN IF ARG =1 (FAC SET TO 0.) **003**
	FLD%   .AG            /LOAD REAL (.AG), **003**
       .DSA   ALOG01+400000  / (ARG)
       JMS*   .EE            /COMPUTE LOG2(ARG)
       FMP%   .AK            /MULTIPLY REAL (.AK)
       .DSA   ALOG07         / (LOGE2)
       JMP*   ALOG           /EXIT
ALOG07        140000         /LOGE(2) (1)  (0.6931471806)
              261344         /        (2)
       .END