t = [0:0.001:10.239]
figure(1)
plot(t,((32-e)*e.**(-t))+e+40,t,(16./(t+0.5))+40)
hchidx = (16./(t+0.5)).+40
hchirp = cos(hchidx.*(2*pi).*t)
echidx = ((32-e)*e.**(-t))+e+40
echirp = cos(echidx.*(2*pi).*t)

figure(2)
plot(t,echirp, t,hchirp)
audiowrite('eh_chirp.wav', transpose([hchirp;echirp]), 19200, 'BitsPerSample', 16)

figure(3)
plot(t,echirp.-hchirp, t,hchirp.-echirp)
audiowrite('ehdelta_chirp.wav', transpose([hchirp-echirp;echirp-hchirp]), 19200, 'BitsPerSample', 16)