%
% clear previous plot meta file
clg
%
% access file of data to be plotted 
load ma2b.out
%
% read first column of data file into array x
n1i=1
n1f=17
x1=ma2b(n1i:n1f,1);
%
% read second column of data file into array y1
y1=ma2b(n1i:n1f,2);
%
load ma2a.out
%
% read first column of data file into array x
n1i=1
n1f=81
x2=ma2a(n1i:n1f,1);
%
% read second column of data file into array y2
y2=ma2a(n1i:n1f,2);
%
% plot x vs y in Cartesian coordinates
axis([0,2,-0.05,1.05])
axis('square')
plot(x1(:,1),y1(:,1),'o',x2(:,1),y2(:,1),'-')
%
% label x and y axes
ylabel('phi1(1,t)')
xlabel('t')
%
% add title to top of plot
title('analytical - solid, numerical - o')
%
% print output
  print ma2p2.eps -deps
%
% send eps file to ihb103 Postscript printer
 !op  -s 1 -d cs1 -q ihb103 ma2p2.eps