fid1 = fopen('SHEAR.TXT');
shearData = fscanf(fid1,'%g %g %g %g %g',[5, inf])';
fid2 = fopen('TRANSVERSE.TXT');
transverseData = fscanf(fid2,'%g %g %g %g %g',[5, inf])';

hFigure = figure; 
plot(shearData(:,1)*100,shearData(:,2),'k-','linewidth',2); hold on;
plot(transverseData(:,1)*100,transverseData(:,2),'k--','linewidth',2); hold on;
legend_h=legend('inplane shear','transverse','location','best');
% set(legend_h,'FontSize',14); %don't do it, see:
% http://www.mathworks.com/support/solutions/en/data/1-16U2V/?solution=1-16U2V
xlabel('\epsilon_2, \epsilon_{12} [%]')
ylabel('\sigma_2, \sigma_6 [MPa]');
hold off;
h = get(gca,'xlabel'); set(h,'FontSize',14);
h = get(gca,'ylabel'); set(h,'FontSize',14);
print(hFigure,'-dps','Fig8_9.eps')
print(hFigure,'-dpng','Fig8_9.png')

fclose(fid1);
fclose(fid2);
