[Python] Sympy running without providing solution given initial conditions

I have 2 differential equations I am trying to solve. Running the code below takes about an hour on my computer before it outputs a very long equation if i do not include the initial conditions. I have yet to see any solutions when running the code with initial conditions after several hours. Is there something I'm doing wrong?

import sympy as sp

# solve = False

r1 = 3570
r2 = 3570
r3 = 6e6
c1 = 200e-12
c2 = 200e-12
c3 = 120e-12

t = sp.symbols("t")

V0 = sp.Function('V0')(t)
V1 = sp.Function('V1')(t)
V2 = sp.Function('V2')(t)

dV1 = sp.Derivative(V1, t)
dV2 = sp.Derivative(V2, t)

R1, R2, R3, C1, C2, C3 = sp.symbols("R1 R2 R3 C1 C2 C3")

eq1 = sp.Eq((C1 + C3)*dV1 + ((1/R1) + (1/R3))*V1 + (-C3)*dV2 + (-1/R3)*V2, (1/R1)*V0)
eq2 = sp.Eq((C3)*dV1 + (1/R3)*V1 + (-C3 - C2)*dV2 + (-1/R2)*V2, 0)

sol = sp.dsolve([eq1, eq2], [V1, V2], ics={V1.subs(t,0): 0, V2.subs(t, 0): 0})

# print(sp.latex(sol[1].subs({"R1": r1, "R2": r2, "R3": r3, "C1": c1, "C2": c2, "C3": c3})))

print(sp.latex(sol[1]))


I've tried simplifying and removing portions of problem to reduce the complexity but I am at the limit of what I can remove.

I found this question that had a comment suggesting wolfram alpha but i don't believe that wolfram can symbolically solve systems of differential equations of this form. Wolfram Alpha Screenshot

Continue reading...