%provides the functions to solve the plastic situation!


function F= pfun(x)

global bf fc Es;
global d As epcp M;

F(1)=0.57*fc*bf*(0.5*x(1)+x(2))-As*Es*x(3);
F(2)=(x(3)+epcp)*x(1)+epcp*x(2)-epcp*d;
F(3)=0.57*fc*bf*(0.5*x(1)*(d-x(2)-x(1)/3)+x(2)*(d-0.5*x(2)))-M*1000000;
end