l=6;
tau=100*(10^(-9));
u1=(5.3/8.68)*(10^(-3));
t=(10^(-9)):(10^(-10)):(2*tau);
e=sin((pi.*t)/(2*tau)).^2;
f1=((4*sqrt(pi))/((u1*l)^2)).*((((u1*l)^2)./(4*pi.*t)).^(3/2)).*exp(-(((u1*l)^2)./(4*pi.*t)));
s2=conv(e,f1);
plot(s2);
%plot(f1);
%plot(e);
grid