function[ys,check] = simple_steadystate(ys,exe);
	global M_
	
		alpha = M_params(1);
		beta = M_.params(2);
		sigma = M_.params(3);
		rho = M_.params(4);
		cs = M_.params(5);
		as = M_.params(6);
		ks = M_.params(7);
		rs = M_.params(8);
		
		check = 0;
		
		rs = (1-beta)/beta;
		ks = (rs/alpha)^(1/(alpha-1));
		cs = rs*ks;
		as = 0;
		c = 0;
		a = 0;
		k = 0;
		r = 0;
		
		ys = [
			cs
			as
			ks
			rs
			c
			a
			k
			r ];