[deiresult,x]=dei(x);

function [deiresult,x]=dei(x);
deiresult=[];

;
%***BEGIN PROLOGUE  DEI
%***PURPOSE  Compute the exponential integral Ei(X).
%***LIBRARY   SLATEC (FNLIB)
%***CATEGORY  C5
%***TYPE      doubleprecision (EI-S, DEI-D)
%***KEYWORDS  EI FUNCTION, EXPONENTIAL INTEGRAL, FNLIB,
%             SPECIAL FUNCTIONS
%***AUTHOR  Fullerton, W., (LANL)
%***DESCRIPTION
%
% DEI calculates the doubleprecision exponential integral, Ei(X), for
% positive doubleprecision argument X and the Cauchy principal value
% for negative X.  If principal values are used everywhere, then, for
% all X,
%
%    Ei(X) = -E1(-X)
% or
%    E1(X) = -Ei(-X).
%
%***REFERENCES  (NONE)
%***ROUTINES CALLED  DE1
%***REVISION HISTORY  (YYMMDD)
%   770701  DATE WRITTEN
%   891115  Modified prologue description.  (WRB)
%   891115  REVISION DATE from Version 3.2
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%***end PROLOGUE  DEI
%***FIRST EXECUTABLE STATEMENT  DEI
deiresult = -de1(-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 DENORM