def analytic_solution(t, A0, alpha, delta, g, K0, n, N0, theta, **params):
k0 = K0 / (A0 * N0)<EOL>lmbda = (g + n + delta) * (<NUM_LIT:1> - alpha)<EOL>ks = ((k0**(<NUM_LIT:1> - alpha) * np.exp(-lmbda * t) +<EOL>(<NUM_LIT:1> / (theta * (g + n + delta))) * (<NUM_LIT:1> - np.exp(-lmbda * t)))**(<NUM_LIT:1> / (<NUM_LIT:1> - alpha)))<EOL>ys = cobb_douglas_output(ks, alpha, **params)<EOL>cs = ((theta - <NUM_LIT:1>) / theta) * ys<EOL>return [ks, cs]<EOL>
Analytic solution for model with Cobb-Douglas production.