%Crintea Cezar 444D/1999
l=1;
u1=(5.3/8.68)*(10^(-3));
t=(10^(-9)):(10^(-11)):(100*(10^(-9)));
f1=((4*sqrt(pi))/((u1*l)^2)).*((((u1*l)^2)./(4*pi.*t)).^(3/2)).*exp(-(((u1*l)^2)./(4*pi.*t)));
c1=f1(1,1:10);
[g1,k]=max(f1);
u2=(2.3/8.68)*(10^(-3));
f2=((4*sqrt(pi))/((u2*l)^2)).*((((u2*l)^2)./(4*pi.*t)).^(3/2)).*exp(-(((u2*l)^2)./(4*pi.*t)));
c2=f2(1,1:10);
[g2,k]=max(f2);
dif=1:length(f2);
for i=1:length(f1)
	dif(i)=abs(f1(i)-(g1/2));
end
for i=1:length(dif)
	if (dif(i) <= 1000)
		disp(i);
		disp(f1(i));
	end
end
linie=ones(size(f1)).*(g1/2);
hold on
plot(f1,'r');
plot(f2,'g');
plot(linie,'b');
hold off