function out = code05b(t,r,p) % % integrand for code05.m % integrand = exp(i*t.^2); if(p==0) out = real(integrand); else out = imag(integrand); end