%
% clear previous plot meta file
clg
%
% access file of data to be plotted 
load h6.out
%
% read first column of data file into array x
n1i=1
n1f=101
x=h6(n1i:n1f,1);
%
% read second column of data file into array y1
y1=h6(n1i:n1f,2);
%
% read third column of data file into array y2
y2=h6(n1i:n1f,3);
%
% plot x vs y in Cartesian coordinates
axis([0,2,-0.2,1.1])
axis('square')
plot(x(:,1),y1(:,1),'--',x(:,1),y2(:,1),'-')
%
% label x and y axes
ylabel('y(zl,t)')
xlabel('t')
%
% add title to top of plot
 title('DSS020 (--), Exact (solid), N = 101')
%title('DSS024 (--), Exact (solid), N = 101')
pause
%
% print output
  print h6p2.eps -deps
% print h6p3.eps -deps
% print h6p4.eps -deps
% print h6p5.eps -deps
% print h6p6.eps -deps
% print h6p7.eps -deps
% print h6p8.eps -deps
% print h6p9.eps -deps

%
% send eps file to ihd09 Postscript printer
 !op  -s 1 -d cs1 -q ihd109 h6p2.eps
%!op  -s 1 -d cs1 -q ihd109 h6p3.eps
%!op  -s 1 -d cs1 -q ihd109 h6p4.eps
%!op  -s 1 -d cs1 -q ihd109 h6p5.eps
%!op  -s 1 -d cs1 -q ihd109 h6p6.eps
%!op  -s 1 -d cs1 -q ihd109 h6p7.eps
%!op  -s 1 -d cs1 -q ihd109 h6p8.eps
%!op  -s 1 -d cs1 -q ihd109 h6p9.eps