.TITLE DLOG	 COMPUTE NATURAL LOGARITHM FOR DOUBLE 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. 01754
/FOR HARDWARE FLOATING POINT PROCESSOR (FPP), DEFINE 90 FPP,
/EDIT 3   1-27-71
/ EDIT 4	28-FEB-75	R.K. HYATT	INSERTED CHECK FOR ARG = 1
/ EDIT 5	20-AUG-75	M. HEBENSTREIT	DISCLAIMER
       .GLOBL DLOG,.DA,.DE,.CAT1D
	.IFUND %FPP
	.GLOBL .AO,.AS
	.ENDC
	.IFUND %FPP
	.DEFIN DLD%,A
	JMS* A			/DOUBLE LOAD (.AO)
	.ENDM
	.ENDC
	.IFDEF %FPP
DLD=713150
	.DEFIN DLD%
	DLD
	.ENDM
	.ENDC
	.IFUND %FPP
	.DEFIN DMP%,A
	JMS* A			/DOUBLE MULTIPLY (.AS)
	.ENDM
	.ENDC
	.IFDEF %FPP
DMP=711540
	.DEFIN DMP%
	DMP
	.ENDM
	.ENDC
/	    CALLING SEQUENCE
/      JMS*   (DLOG)	SUBR CALL
/      JMP    .+2		SUBR RETURN
/      CAL/XCTADDR		ADDR OF ARGUMENT (XCT IF INDIRECT)
/
DLOG   CAL    0	         /ENTRY-EXIT
       JMS*   .DA	         /GET ARGUMENT (.DA)
       JMP    DLOG02
DLOG01 CAL    0	         /ARGUMENT
DLOG02	LAC	DLOG01		/CHECK FOR ARG =1 **004**
	JMS*	.CAT1D			/**004**
	SKP				/ARG DOESN'T =1 SO COMPUTE LOG
	JMP*	DLOG			/ARG = 1 SO RETURN (FAC SET TO 0.) **004**
	DLD%   .AO	         /LOAD DOUBLE (.AO), **004**
       .DSA   DLOG01+400000
       JMS*   .DE	         /COMPUTE LOG2(ARG) (.DE)
       DMP%   .AS	         /MULTIPLY DOUBLE (.AS)
	.DSA	DLOG07	/ *** DDS JAN69 ***
       JMP*   DLOG	         /EXIT
DLOG07	    0	         / LOGE(2) (1)  (0.6931471806)
	    261344         /         (2)
	    137700         /         (3)
       .END