function F = findsteady_RL(x,ypsilon,gamma_l,beta_p,RDss,rho_s,eta)
cell1 = (ypsilon*(1-gamma_l))*((x(1)-RDss)/(1-beta_p*gamma_l*(1-rho_s)))-x(1);
 F = cell1;