[ar,ai,br,bi]=zexp(ar,ai,br,bi);

function [ar,ai,br,bi]=zexp(ar,ai,br,bi);
%***BEGIN PROLOGUE  ZEXP
%***SUBSIDIARY
%***PURPOSE  Subsidiary to ZBESH, ZBESI, ZBESJ, ZBESK, ZBESY, ZAIRY and
%            ZBIRY
%***LIBRARY   SLATEC
%***TYPE      ALL (ZEXP-A)
%***AUTHOR  Amos, D. E., (SNL)
%***DESCRIPTION
%
%     doubleprecision COMPLEX EXPONENTIAL FUNCTION B=EXP(A)
%
%***SEE ALSO  ZAIRY, ZBESH, ZBESI, ZBESJ, ZBESK, ZBESY, ZBIRY
%***ROUTINES CALLED  (NONE)
%***REVISION HISTORY  (YYMMDD)
%   830501  DATE WRITTEN
%   910415  Prologue converted to Version 4.0 format.  (BAB)
%***end PROLOGUE  ZEXP
persistent ca cb zm ; 

if isempty(zm), zm=0; end;
if isempty(ca), ca=0; end;
if isempty(cb), cb=0; end;
%***FIRST EXECUTABLE STATEMENT  ZEXP
zm = exp(ar);
ca = zm.*cos(ai);
cb = zm.*sin(ai);
br = ca;
bi = cb;
end
%DECK ZKSCL