function [k] = steadyk(beta, h, ki, delta, theta, phi)


k = ((delta + (1/beta)- 1)/(theta*(ki^phi)*(h^(1-theta-phi))))^(1/(theta-1));


end

