Useful MatLab code for plotting
Increase all text size:
set(gca,'FontSize',14)
figureHandle = gcf;
set(findall(figureHandle,'type','text'),'fontSize',14)
Labelling the plot:
xlabel('$\alpha_\varepsilon$ - direction of strain probe $(\circ)$','interpreter','latex');
ylabel('Second order work $W_2$ and $W_2^{\mu}$ (J)','interpreter','latex');
Hatch region:
harea = area([140 140 150 150 140],[-0.05 0.3 0.3 -0.05 -0.05],3 );
child=get(harea,'Children')
set( harea, 'FaceColor', 'm')
alpha('texture')
Increase text size for plotyy function
[ax h1 h2]=plotyy(e22,ev,e22,q)
set(ax,'FontSize',14)
Subplot
%subplot of (w2 normalized), (compare), (new,lostcontact), (total contact)
subplot(2,2,1);
polar(degtorad(alpha),w2normal,'r--*');
title('Normalized $W_2$ at $\eta=0.92$','interpreter','latex')
subplot(2,2,2);
plot(alpha,ntotal,'b--*');
xlim([0 350]);
xlabel('Angle of stress probe - $\alpha_\sigma$','interpreter','latex');
ylabel('Total contacts - $N_c$','interpreter','latex');
grid on;
title('Total contacts numbers $N_c$ at $\alpha_\sigma$','interpreter','latex')
subplot(2,2,3);
plot(alpha,nlost,'r--o',alpha,nnew,'b--*');
xlim([0 350]);
xlabel('Angle of stress probe - $\alpha_\sigma$','interpreter','latex');
ylabel('New ($N_{new}$) and lost ($N_{lost}$) contacts','interpreter','latex');
legend('N_{lost}','N_{New}');
grid on;
title('Lost contacts and new contacts during the directional research')
subplot(2,2,4);
plot(alpha,w2macro,'b--o',alpha,w2micro,'r--*');
legend('W_2 macro','W_2 micro');
xlabel('Angle of stress probe - $\alpha_\sigma$','interpreter','latex');
ylabel('Second order work (J)','interpreter','latex');
xlim([0 350]);
grid on;
title('Comparison $W_{2}macro$ and $W_{2}micro$','interpreter','latex')
Crocodile & Saviour 