clc
clear all

b_w=0.15;
d=0.173
f_c=40;
[M_ED]=M_value;


for i=1:10
    
disp([' step: ', num2str(i)]);  

iter = 0;   
tolerance =0 ;
Epsilon = 0;
M_Eds = M_ED(i)
value=1

     for e_s1 = 0.0: 0.1: 20
                for e_c2 = -0.1: -0.1: -3.5
                    disp([' ec2: ', num2str(e_c2)]); 
                    disp([' es1: ', num2str(e_s1)]);
                    
                    zeta = abs(e_c2)/(abs(e_c2)+abs(e_s1));
                    x= zeta*d;

                    if(e_c2 <=0 && e_c2 >=-2.0)
                    alpha_r=-0.5*(e_c2+(e_c2^2/6));
                    k_a=(e_c2+8)/(4*(e_c2+6));
                    end

                    if(e_c2 <=-2.0 && e_c2 >=-3.5)
                    alpha_r=(2/(3*e_c2))+1;
                    k_a=((3*e_c2^2)+(4*e_c2)+2)/((6*e_c2^2)+(4*e_c2));
                    end


                    F_cd=b_w*x*alpha_r*f_c;
                    a=k_a*x;
                    z=d-a;

                    M_rd=F_cd*z;


                     % CONVERGENCE CHECK

                     tolerance =M_Eds-M_rd 


                         if (abs(tolerance - Epsilon) <= 0.1)
                           disp('value is ok');
                               Es1(i)=e_s1;
                               Ec2(i)=e_c2;
                               [Es1', Ec2']
                               value=0;
                               break;
                                                         
                         else
                          disp('valve is checking');
                         end
                end
        if(value==0)
          break;
        end
     end
      
      
  


end