function z=rbc_steady(x); alp=0.35; bet=0.99; delt=0.025; A=1; gam=1;1.007; mu=0.3; g=0.2; sig=1.5; y=x(1); c=x(2); k=x(3); h=x(4); i=x(5); r=x(6); w=x(7); lam=x(8); z=zeros(8,1); z(1)=((1-mu)*c)/(mu*(1-h))-w; z(2)=r+1-delt-(1/bet); z(3)=alp*A*(k/h)^(alp-1)-r; z(4)=(1-alp)*A*(k/h)^(alp)-w; z(5)=(1-g)*A*h*(k/h)^(alp)-c-i; z(6)=gam-1+delt-(i/k); z(7)=mu*(c)^(mu*(1-sig)-1)*(1-h)^((1-mu)*(1-sig))-lam; z(8)=A*h*(k/h)^(alp)-y;