% semi-infinite heat conduction
% response to a   pulse
clear

dt=0.1
for n=1:200
   t(n)=n*dt;tn=1/t(n);
   y(n)=sqrt(tn)*exp(-tn);
end

figure(1); clf;
 plot(t,y,'linewidth', 1.0)
 xlabel('t^* = \tau/t = x^2/4\alphat')
 ylabel('\eta')
 grid on

% fig=figure(1);
% print(fig,fullfile(pwd, 'Figure6_3.tif'),'-dtiff','-r600')