%
% clear previous plot
clg
%
% load output from Fortran program for plot of component
% concentrations vs time
load ma3c.out;
%
% store time
np=101
for i=1:np;
   tp(i)=ma3c(i,1);
end;
%
% store component 1 concentration at z = 0
for i=1:np;
   cp(i,1)=ma3c(i,2);
end;
%
% store component 2 concentration at z = 0
for i=1:np;
   cp(i,2)=ma3c(i,3);
end;
%
% store component 3 concentration at z = 0
for i=1:np;
   cp(i,3)=ma3c(i,4);
end;
%
% plot concentrations vs time
figure(1);
plot(tp,cp(:,1),'-',tp,cp(:,2),'-',tp,cp(:,3),'-')
%
% label axes
xlabel(' t');
ylabel(' x1, x2, x3 at z = 0');
title(' Transient Multicomponent Diffusion')
pause;

%
% Encapsulated Postscript output file
  print ma3p2.eps -deps
%
% store component 1 concentration at z = zl
for i=1:np;
   cp(i,4)=ma3c(i,5);
end;
%
% store component 2 concentration at z = zl
for i=1:np;
   cp(i,5)=ma3c(i,6);
end;
%
% store component 3 concentration at z = zl
for i=1:np;
   cp(i,6)=ma3c(i,7);
end;
%
% plot concentrations vs time
figure(2);
plot(tp,cp(:,4),'-',tp,cp(:,5),'-',tp,cp(:,6),'-')
%
% label axes
xlabel(' t');
ylabel(' x1, x2, x3 at z = zl');
title(' Transient Multicomponent Diffusion')
pause;
%
% Encapsulated Postscript output file
  print ma3p3.eps -deps
%
% load output from Fortran program for plot of component
% fluxes vs time
load ma3f.out;
%
% store component 1 flux 
for i=1:np;
   fp(i,1)=ma3f(i,2);
end;
%
% store component 2 flux 
for i=1:np;
   fp(i,2)=ma3f(i,3);
end;
%
% store component 3 flux 
for i=1:np;
   fp(i,3)=ma3f(i,4);
end;
%
% plot fluxes vs time at z = 0
figure(3);
plot(tp,fp(:,1),'-',tp,fp(:,2),'-',tp,fp(:,3),'-')
%
% label axes
xlabel(' t');
ylabel(' n1, n2, n3');
title(' Transient Multicomponent Diffusion')
pause;
%
% Encapsulated Postscript output file
  print ma3p4.eps -deps
%
% send Postscript file to ihd09 Postscript printer
 !op -s 1 -d cs1 -q ihb103 ma3p2.eps
 !op -s 1 -d cs1 -q ihb103 ma3p3.eps
 !op -s 1 -d cs1 -q ihb103 ma3p4.eps