function F = steadystate(lam,Cstarr,thetstar,psii,m,mi,alphaa,a,deltaa,chii,rhoo,Rw)
F=[lam(1)-(thetstar/lam(3))^alphaa*(lam(2))^(1-2*a)*Cstarr;
    lam(2)-(lam(4)/((1-alphaa)*lam(1)))^(1/a)*lam(3)^(-alphaa/a);
    lam(3)-(lam(5)/(alphaa*lam(1)))^(1/(alphaa-1));
    (1/m)-Rw^rhoo*(lam(1)*lam(7)^psii)^(1-rhoo)*lam(3)^alphaa*lam(2)^a/lam(8)^chii;
    (1/mi)-lam(1)*lam(6)^psii*lam(3)^(alphaa-1)/(deltaa*lam(8))^chii;
    lam(6)-lam(5)/(deltaa*lam(8))^chii;
    lam(7)-lam(4)/lam(8)^chii;
lam(8)-0.18*lam(4)];
end