function out = code05a(t,r,p) % % integrand for code05.m % integrand = (-2*i/pi)*exp(i*r*cosh(t)); if(p==0) out = real(integrand); else out = imag(integrand); end