def analytic_solution(t, k0, alpha, delta, g, n, s, **params):
lmbda = (g + n + delta) * (<NUM_LIT:1> - alpha)<EOL>ks = (((s / (g + n + delta)) * (<NUM_LIT:1> - np.exp(-lmbda * t)) +<EOL>k0**(<NUM_LIT:1> - alpha) * np.exp(-lmbda * t))**(<NUM_LIT:1> / (<NUM_LIT:1> - alpha)))<EOL>return ks<EOL>
Analytic solution for model with Cobb-Douglas production.