function z_ss=gwgrowth_steady_state_helper(params)
sigma = 0.5;
alpha = 1/3; 
r = 0.065;
beta = 1/(1+r);
psi = 1;
delta = 0.025;
rho = 0.95;
g = 1/rho - 1;
theta_ss = 1;  %(1)
M = 3;

params = [sigma alpha r beta psi delta rho g M];


options=optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);

z_ss = fsolve(@(z_ss)steadygrow(params, z_ss),  zeros(params(9),1), options); 

