[casinhresult,z]=casinh(z);

function [casinhresult,z]=casinh(z);
casinhresult=[];
persistent ci firstCall ; if isempty(firstCall),firstCall=1;end; 

;
%***BEGIN PROLOGUE  CASINH
%***PURPOSE  Compute the arc hyperbolic sine.
%***LIBRARY   SLATEC (FNLIB)
%***CATEGORY  C4C
%***TYPE      COMPLEX (ASINH-S, DASINH-D, CASINH-C)
%***KEYWORDS  ARC HYPERBOLIC SINE, ASINH, ELEMENTARY FUNCTIONS, FNLIB,
%             INVERSE HYPERBOLIC SINE
%***AUTHOR  Fullerton, W., (LANL)
%***DESCRIPTION
%
% CASINH(Z) calculates the complex arc hyperbolic sine of Z.
%
%***REFERENCES  (NONE)
%***ROUTINES CALLED  CASIN
%***REVISION HISTORY  (YYMMDD)
%   770401  DATE WRITTEN
%   861211  REVISION DATE from Version 3.2
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%***end PROLOGUE  CASINH
if isempty(ci), ci=0; end;
if firstCall,   ci=[complex(0.,1.)];  end;
firstCall=0;
%***FIRST EXECUTABLE STATEMENT  CASINH
casinhresult = -ci.*casin(ci.*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 CASYI