% function[sqz] = csqrt(z)
% Complex square root using positive real axis branch cut
% (EJP branch cut)
function[sqz] = csqrt(z)
% Pekeris branch cut
if real(z) < 0
% EJP Branch cut
%if imag(z) < 0
 sqz=i*sqrt(-z);
else
 sqz=sqrt(z);
end
