function y=envelope(t)
y=(log(t+0.001)-log(0.001)).*exp(-8*t);
