function F = somevarsolve(x,DD,omega,N_nrss,C_nrss,w_rss,epsilon_wr,sigma,phi)
    F=[DD*omega*x(2)+DD*(1-omega)*N_nrss-omega*x(1)-(1-omega)*C_nrss;
        w_rss*((epsilon_wr-1)/epsilon_wr)*(x(1)^(-sigma))-(x(2)^phi)];
end