## Gibberish

### 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')