синтаксис для решения системы дифференциальных уравнений в sympy

Question

Я новичок в симпы и в процессе обучения. Я просматривал документацию и вопросы в обмене стеками относительно символического решения системы дифференциальных уравнений с начальными условиями с использованием sympy.

У меня есть простая система ОДУ-s

(dV/dt) = -( 1/RC) * (V(t)) + I(t)/C 

(dI/dt) = -(R1/L) * (I(t)) - (1/L) * V(t) + Vs/L 

с начальными условиями V(0) = V0 и I(0) = I0

Я просмотрел много вопросов в обмен стека, а не успешным в поиске подходящего ответа. Было бы очень полезно, если бы кто-нибудь мог показать мне синтаксис, чтобы войти в систему связанных дифференциальных уравнений с начальными условиями.

Answer 1

Поддержка системы ODE поддерживается только в версии SymPy. Он будет добавлен в 0.7.6. Синтаксис будет

V, I = symbols("V I", cls=Function) 
RC, t, C, Vs, L, R1, V0, I0 = symbols("RC t C Vs L R1 V0 I0") 
system = [Eq(V(t).diff(t), -1/RC*V(t) + I(t)/C), Eq(I(t).diff(t), -R1/L*I(t) - 1/L*V(t) + Vs/L)] 
ics = {V(0): V0, I(0): I0} 
dsolve(system, [V(t), I(t)], ics=ics) 

Кажется, что есть ошибка, которая предотвращает это от работы в текущем мастере SymPy, если я что-то неправильно набранном (https://github.com/sympy/sympy/issues/8193).