[cgammaresult,z]=cgamma(z);

function [cgammaresult,z]=cgamma(z);
cgammaresult=[];

;
%***BEGIN PROLOGUE  CGAMMA
%***PURPOSE  Compute the complete Gamma function.
%***LIBRARY   SLATEC (FNLIB)
%***CATEGORY  C7A
%***TYPE      COMPLEX (GAMMA-S, DGAMMA-D, CGAMMA-C)
%***KEYWORDS  COMPLETE GAMMA FUNCTION, FNLIB, SPECIAL FUNCTIONS
%***AUTHOR  Fullerton, W., (LANL)
%***DESCRIPTION
%
% CGAMMA(Z) calculates the complete gamma function for COMPLEX
% argument Z.  This is a preliminary version that is portable
% but not accurate.
%
%***REFERENCES  (NONE)
%***ROUTINES CALLED  CLNGAM
%***REVISION HISTORY  (YYMMDD)
%   770701  DATE WRITTEN
%   861211  REVISION DATE from Version 3.2
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%***end PROLOGUE  CGAMMA
%***FIRST EXECUTABLE STATEMENT  CGAMMA
cgammaresult = exp(clngam(z));
%
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',z); evalin('caller',[inputname(1),'=FUntemp;']); end
end
%DECK CGAMR