def equilibrium_capital(alpha, delta, g, n, rho, theta, **params):
    """Steady state value for capital stock (per unit effective labor)."""
    return (alpha / (delta + rho + theta * g))**(1 / (1 - alpha))