%
% clear previous plot meta file
clg
%
% access file of data to be plotted
load m2.out
%
% read data into a series of vectors
n1i=1
n1f=31
nin=30
x11=m2(n1i:n1f,1);
y11=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y12=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y13=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y14=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y15=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y16=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y17=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y18=m2(n1i:n1f,2);
n1i=1+n1f;
n1f=n1i+nin;
y19=m2(n1i:n1f,2);
%
% plot x vs y in Cartesian coordinates
axis([0,1,0,1])
axis('square')
plot(x11,y11,'-',..
     x11,y12,'-',..
     x11,y13,'-',..
     x11,y14,'-',..
     x11,y15,'-',..
     x11,y16,'-',..
     x11,y17,'-',..
     x11,y18,'-',..
     x11,y19,'-')
%
% label x and y axes
ylabel('phi(r,tau)')
xlabel('r')
%
% add title to top of plot
title('phi(r,tau) vs r, tau = 0, 0.05, 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 1')
pause
hold
%
% save plot image in meta file figur.met
meta figur.met
%
% check for meta file; ! can be used to execute a Unix command within
% Matlab
!ls *.met
%
% the meta file figur.met created by mat.m is converted to a Postscript
% file, figur.ps
!gpp figur.met -dps
%
% check for .ps file
!ls *.ps
%
% send Postscript file to ad109 Postscript printer
!op -d cs1 -q ihd109 figur.ps
%
% delete .met files
!rm *.met
%
% delete .ps files
!rm *.ps
%
% check if Postscript plotting is finished
!oq -d cs1
