In [1]:
%display latex
In [2]:
var('t, s, alpha, delta, A, L, theta, omega')
Out[2]:
In [3]:
k = function('k')(t)
In [4]:
de = diff(k, t) == s * k ^ alpha - k * (delta + A + L + 2 * theta * sin(2 * pi * t / omega))
show(de)
In [5]:
solution = desolve(de, k, ivar=t, show_method=True)
show(solution)