function [dliresult,x]=dli(x);
dliresult=[];
;
%***BEGIN PROLOGUE DLI
%***PURPOSE Compute the logarithmic integral.
%***LIBRARY SLATEC (FNLIB)
%***CATEGORY C5
%***TYPE doubleprecision (ALI-S, DLI-D)
%***KEYWORDS FNLIB, LOGARITHMIC INTEGRAL, SPECIAL FUNCTIONS
%***AUTHOR Fullerton, W., (LANL)
%***DESCRIPTION
%
% DLI(X) calculates the doubleprecision logarithmic integral
% for doubleprecision argument X.
%
%***REFERENCES (NONE)
%***ROUTINES CALLED DEI, XERMSG
%***REVISION HISTORY (YYMMDD)
% 770701 DATE WRITTEN
% 890531 Changed all specific intrinsics to generic. (WRB)
% 890531 REVISION DATE from Version 3.2
% 891214 Prologue converted to Version 4.0 format. (BAB)
% 900315 CALLs to XERROR changed to CALLs to XERMSG. (THJ)
%***end PROLOGUE DLI
%***FIRST EXECUTABLE STATEMENT DLI
if( x<=0.0D0 )
xermsg('SLATEC','DLI','LOG INTEGRAL UNDEFINED FOR X LE 0',1,2);
end;
if( x==1.0D0 )
xermsg('SLATEC','DLI','LOG INTEGRAL UNDEFINED FOR X = 0',2,2);
end;
%
[dliresult ]=dei(log(x));
%
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',x); evalin('caller',[inputname(1),'=FUntemp;']); end
end
%DECK DLLSIA
