l=1.2;
tau=80*(10^(-9));b=18;
u1=(5.3/8.68)*(10^(-3));p=2*10^(-10);m=b*tau/p;
t=p:p:(b*tau);t1=t(1,1:2*tau/p);
e1=sin((pi.*t1)/(2*tau)).^2;x=(b-2)*tau/p;v=zeros(1,x);e=[e1 v];
f1=((4*sqrt(pi))/((u1*l)^2)).*((((u1*l)^2)./(4*pi.*t)).^(3/2)).*exp(-(((u1*l)^2)./(4*pi.*t)));


e2= ones(size(t1));e3=[e2 v];
ps4=conv(e3,f1);s5=p.*s4;
%plot(s5)
%hold on
s2=conv(e,f1);s3=p.*s2;
%plot(s3)
%grid on;
% zoom on
s=length(t1);v=s2(1:1,1:s);
plot(t1,v)