function f = ss(x,beta,theta,delta,rho);
% For simplicity, rename the variables from the x vector
c = x(1);
k = x(2);
z = x(3);
% f = [resource constraint at steady-state, euler at steady-state, law of motion]
% Equations organized such that they're equal to zero.
f = [c+k-(z*k^theta + (1-delta)*k);
1/c-beta*((theta*z*k^(theta-1)+(1-delta))/c);
z-(1-rho)-rho*z];