%subplot of (w2 normalized), (compare), (new,lostcontact), (total contact)
%need these variable:
%(alpha)(nperm)(nlost)(nnew)(w2macro)(w2micro)
figure;
subplot(2,2,1);
% theta  = [0:10:350];
% 
% h_fake=polar(degtorad(theta),r);

% set(P, 'Visible', 'off');
% hold on;
redLine(1:37)=0.5;
theta=[0:10:360];
polar2(degtorad(alpha),w2normal,[0 1.5],'r--o');
hold on;
polar2(degtorad(theta),redLine,'r-');
hold off;
title('Normalized $W_2^n$ at $\eta=0.92$','interpreter','latex');


subplot(2,2,2);
plot(alpha,nperm,'b--*');
xlim([0 350]);
xlabel('Angle of stress probe - $\alpha_\sigma$','interpreter','latex');
ylabel('Permanent contacts - $N_p$','interpreter','latex');
grid on;
title('Permanent contacts $N_p$ 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','interpreter','latex')

subplot(2,2,4);
plot(alpha,w2macro,'b--o',alpha,w2micro,'r--*');
legend('W_2','W_2^m');
xlabel('Angle of stress probe - $\alpha_\sigma$','interpreter','latex');
ylabel('Second order work (J)','interpreter','latex');
xlim([0 350]);
grid on;
title('Comparison between $W_{2}$ and $W_{2}^m$ (calculated by microscopic variables)','interpreter','latex')


set(gca,'FontSize',14)
figureHandle = gcf;
set(findall(figureHandle,'type','text'),'fontSize',14)