% drawTRDiffusion.m
% needs
%	Uo
%	alpha
%	time
%	rmax, zmax, dr, dz
%	c


if time > 0


subplot(2,1,1)
clear u y
hold on
u  =(-rmax:dr:rmax);
y = 2*Uo*exp(-u.^2/(4*alpha*time)) / (4*pi*alpha*time)^(3/2);
plot(u,y,'b-')

subplot(2,1,2)
clear u y
hold on
u=(0:dz:zmax);
y = 2*Uo*exp(-u.^2/(4*alpha*time)) / (4*pi*alpha*time)^(3/2);
hold on
plot(u,y,'b-')

end
