%
% store solution
load h1.out
%
% store solution over time and space
ni=1;
nf=51;
%
% space grid, temperatures
   xp=h1(ni:nf,1);
   u1=h1(ni:nf,2);
   ni=nf+1;
   nf=ni+50;
   u2=h1(ni:nf,2);
   ni=nf+1;
   nf=ni+50;
   u3=h1(ni:nf,2);
   ni=nf+1;
   nf=ni+50;
   u4=h1(ni:nf,2);
   ni=nf+1;
   nf=ni+50;
   u5=h1(ni:nf,2);
%
% plot solution
plot(xp,u1,'-',xp,u2,'-',xp,u3,'-',xp,u4,'-',xp,u5,'-');
hold on;
%
% label axes
xlabel(' x');
ylabel(' T(x,t)');
title(' erf, T0 = 25, Ts = 250, alpha = 0.001, t = 100, 200, ..., 500');
%
% write .ps or .eps file
% print output.ps
  print h1p1.eps -deps
%
% send .ps or .eps file to ihb103 Postscript printer
!op -s 1 -d cs1 -q ihb103 h1p1.eps
%printfile output.ps