%
% clear previous plot
clg
%
% load output from Fortran program for plot of solution
load h4r.out;
%
% number of t and r values in file h4r.out
nt=11;
nr=11;
%
% store r
for i=1:nr;
   rp(i)=h4r(i,1);
end;
%
% store solution as a function of r and t for z = zl
for j=1:nt;
for i=1:nr;
    up(i,j)=h4r((j-1)*nr+i,2);
end;
end;
%
% plot solution
figure(1);
plot(rp,up,'-')
hold on;
%
% label axes
xlabel(' r');
ylabel(' u(r,zl,t)');
title(' Dynamic Graetz Problem, zl = 30')
%
% Postscript or Encapsulated Postscript output file
% print output.ps
  print h4p2a.eps -deps
%
% store solution as a function of r and t for z = zl
for j=1:nt;
for i=1:nr;
    up(i,j)=h4r(nr*nt+(j-1)*nr+i,2);
end;
end;
%
% plot solution
figure(2);
plot(rp,up,'-')
hold on;
%
% label axes
xlabel(' r');
ylabel(' u(r,zl,t)');
title(' Dynamic Graetz Problem, zl = 90')
%
% Postscript or Encapsulated Postscript output file
% print output.ps
  print h4p2b.eps -deps
%
% load output from Fortran program for plot of solution
load h4z.out;
%
% number of t and z values in file h4z.out
nt=11;
nz=31;
%
% store z
for i=1:nz;
   zp(i)=h4z(i,1);
end;
%
% store solution as a function of z and t for r = 0
for j=1:nt;
for i=1:nz;
    up(i,j)=h4z((j-1)*nz+i,2);
end;
end;
%
% plot solution
figure(3);
plot(zp,up,'-')
hold on;
%
% label axes
xlabel(' z');
ylabel(' u(0,z,t)');
title(' Dynamic Graetz Problem, zl = 30')
%
% Postscript or Encapsulated Postscript output file
% print output.ps
  print h4p3a.eps -deps
%
% store solution as a function of z and t for r = 0
for j=1:nt;
for i=1:nz;
    up(i,j)=h4z(nz*nt+(j-1)*nz+i,2);
end;
end;
%
% plot solution
figure(4);
plot(zp,up,'-')
hold on;
%
% label axes
xlabel(' z');
ylabel(' u(0,z,t)');
title(' Dynamic Graetz Problem, zl = 90')
%
% Postscript or Encapsulated Postscript output file
% print output.ps
  print h4p3b.eps -deps
%
% send Encapsulated Postscript file to ihb103 Postscript printer
!op -s 1 -d cs1 -q ihb103 h4p2a.eps
!op -s 1 -d cs1 -q ihb103 h4p2b.eps
!op -s 1 -d cs1 -q ihb103 h4p3a.eps
!op -s 1 -d cs1 -q ihb103 h4p3b.eps