l=1.5;
tau=100*(10^(-9));b=5;
u1=(5.3/8.68)*(10^(-3));p=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)));
%plot(s2);
%plot(f1);
%plot(e);
s=zeros(1,length(t));
for i=1:length(t)
	c=e(1,1:i);
	d=rot90(c,-2);h=zeros(1,length(t)-i);w=[d h];
	r=f1.*w;
	s(i)=sum(r);
end
s22=p.*s;
plot(s22)