fcA=30;
gammamA=1.5;
bienA=2.44*sqrt(fcA/gammamA)/10000;
khacA=0.67*fcA/gammamA;
x1=0:0.0001:0.0035;
paramsA(1)=bienA;
paramsA(2)=fcA;
paramsA(3)=gammamA;
paramsA(4)=khacA;
fcP=30;
gammamP=1.5;
bienP=0.002;
khacP=0.85*fcP/gammamP;
paramsP(1)=bienP;
paramsP(2)=fcP;
paramsP(3)=gammamP;
paramsP(4)=khacP;
plot(x1,Anh(x1,paramsA),'o-r','LineWidth',1,'MarkerSize',4);
text(0.0004,Anh(0.0004,paramsA),'Anh\rightarrow',...
	'HorizontalAlignment','right');
hold on
plot(x1,Phap(x1,paramsP),'g*-','LineWidth',1,'MarkerSize',4);
text(0.0007,Phap(0.0007,paramsP),'\leftarrowPhap');
%text(0.0003,Phap(0.0003,paramsP),'\bullet\leftarrow\fontname...
%{times}fcu/gama{\ite}^{0.67{\itt}}at{\itt}=0.0003);
axis('auto');
hold off;
xlabel('Strain');
ylabel('Stress');
title('Short term design stress - strain curve for norman weight concrete');
grid off;
zoom on;