clg
load h2.out
 ni=35
 n0=1
 n1=n0+ni
 t1=h2(n0:n1,1);
y11=h2(n0:n1,2);
y12=h2(n0:n1,3);
y13=h2(n0:n1,4);
y14=h2(n0:n1,5);
y15=h2(n0:n1,6);
y16=h2(n0:n1,7);
n2=n1+1
n3=n2+ni
y21=h2(n2:n3,2);
y22=h2(n2:n3,3);
y23=h2(n2:n3,4);
y24=h2(n2:n3,5);
y25=h2(n2:n3,6);
y26=h2(n2:n3,7);
n4=n3+1
n5=n4+ni
y31=h2(n4:n5,2);
y32=h2(n4:n5,3);
y33=h2(n4:n5,4);
y34=h2(n4:n5,5);
y35=h2(n4:n5,6);
y36=h2(n4:n5,7);
n6=n5+1
n7=n6+ni
y41=h2(n6:n7,2);
y42=h2(n6:n7,3);
y43=h2(n6:n7,4);
y44=h2(n6:n7,5);
y45=h2(n6:n7,6);
y46=h2(n6:n7,7);
n8=n7+1
n9=n8+ni
y51=h2(n8:n9,2);
y52=h2(n8:n9,3);
y53=h2(n8:n9,4);
y54=h2(n8:n9,5);
y55=h2(n8:n9,6);
y56=h2(n8:n9,7);
% Use second time scale
ni=32
n10=n9+1
n11=n10+ni
 t2=h2(n10:n11,1);
y61=h2(n10:n11,2);
y62=h2(n10:n11,3);
y63=h2(n10:n11,4);
y64=h2(n10:n11,5);
y65=h2(n10:n11,6);
% Note that y66 is not used since this would require log10(0)
semilogy(t1,[y11,y12,y13,y14,y15,y16,...
             y21,y22,y23,y24,y25,y26,...
             y31,y32,y33,y34,y35,y36,...
             y41,y42,y43,y44,y45,y46,...
             y51,y52,y53,y54,y55,y56],'-',...
         t2,[y61,y62,y63,y64,y65],'-')
ylabel('(T - Ts)/(To - Ts)')
xlabel('t*alpha/(ro**2)')
title('Gurney Lurie chart for a cylinder')
pause
%
% write .ps or .eps file
% print output.ps
  print h2cylr.esp -deps
%
% send .ps or .eps file to printer
%!op -s 1 -d cs1 -q ihb103 h2cylr.eps
%!op -s 1 -d cs1 -q ihd109 h2cylr.eps
%printfile output.ps