Hi,
I am a bit confused about the state of solving ODE's in sympy. It seems
that I cannot directly solve initial value problems, right?
Is there a good way to solve the resulting equation system in simple cases?
e.g. my code right now:
t = sym.Symbol('t')
i = sym.Function('i')
prob = Lsym.Derivative( i(t), t, 2 ) + sym.Derivative(i(t), t) + i(t)
sol = sym.dsolve(prob, i(t))
now I'd like to do something like:
equations = [sol.subs(t,0) == var(1), sym.diff(sol,t).subs(t,0) == var(2)]
solve(equations)
sym.diff(sol,t) does not return a fully substituted result however. Instead
I receive:
Subs(Derivative(i(_t) == (_t*C2 + C1)*exp(-_t), _t), (_t,), (0,))
Can I still solve with this somehow or do I need to rethink my approach?
Best Regards,
Matthias
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To view this discussion on the web visit
https://groups.google.com/d/msg/sympy/-/0eDn9XsXbZEJ.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sympy?hl=en.