fig2_article.m 969 Bytes
clear all;
close all;



fig=figure;

fig.Units='inches';
fig.Position=[7.0312 6.6875 4.0312 3.3021];
fig.OuterPosition=[7.0312 6.6875 4.0312 4.3021];

fig.PaperUnits = 'inches';
fig.PaperPosition = [0 0 4.0312 3.3021];
fig.PaperPositionMode = 'manual';
fig.PaperOrientation='landscape';
fig.PaperSize=[4.0312 3.3021];


strain=[0:1e-6:3e-4];
prob=exp(-2*3e-5^2./strain.^2);


 plot(strain(find(prob<0.37)),prob(find(prob<0.37)),'k--','linewidth',1.5)
hold on
plot(strain(find(prob>0.37)),prob(find(prob>0.37)),'k-','linewidth',1.5)

a=xlabel('$E_s$');
b=ylabel('$P_{\varepsilon}=P(E_w>\varepsilon_c)$');
set(a,'interpreter','latex','fontsize',11)
set(b,'interpreter','latex','fontsize',11)

scrit=3e-5;
probcrit=0.37;
plot(xlim,[probcrit probcrit])
plot([scrit scrit],ylim)
c=text(1.5e-4,0.4,'$P_c$');
d=text(3.5e-5,0.65,'$\varepsilon_c$');
set(c,'interpreter','latex','fontsize',11)
set(d,'interpreter','latex','rotation',90,'fontsize',11)

saveas(fig,'Pe_Es','pdf')