%display latex
var('t, s, alpha, delta, A, L, theta, omega')
k = function('k')(t)
de = diff(k, t) == s * k ^ alpha - k * (delta + A + L + 2 * theta * sin(2 * pi * t / omega))
show(de)
solution = desolve(de, k, ivar=t, show_method=True)
show(solution)