ConceptPythonBeginner Computer algebra: formulas as trees, exact until you ask for digits sol = sp.dsolve(sp.Eq(A(t).diff(t), -k * A(t) ** 2), A(t), ics={A(0): A0})
ToolPythonAdvanced Perturbation theory with SymPy: the Duffing oscillator's frequency shift x1_naive = sp.dsolve(naive, x1(t), ics={x1(0): 0, x1(t).diff(t).subs(t, 0): 0}).rhsx1_sol = sp.dsolve(eq1, x1(tau), ics={x1(0): 0, x1(tau).diff(tau).subs(tau, 0): 0}).rhs