from sympy import *
s, Vpk, C1, C2, C3, R1, R2, R3 = symbols('s Vpk C1 C2 C3 R1 R2 R3', real=True)
init_printing(use_unicode=True)
zC1 = 1 / (s*C1)
zC2 = 1 / (s*C2)
zC3 = 1 / (s*C3)
z22 = R2 + zC2
z33 = R3 + zC3
z2233 = z33*z22/(z33+z22)
z22331 = zC1+z2233
V1 = Vpk / s
VM1 = V1 * z22331/(R1+z22331)
VM1
import numpy as np
import pandas as pd
Vsimul = VM1 \
.subs([(R1, 10**6), (R2, 10**6), (R3, 0)]) \
.subs([(C1, 10*10**-12), (C2, 10*10**-12), (C3, 10*10**-12)]) \
.subs([(Vpk, 3.3)])
# Simplify expression
Vsimul = simplify(Vsimul, rational=True)
t = symbols('t', positive=True)
Vt = inverse_laplace_transform(Vsimul, s, t, noconds=True)
Vt
ts = np.linspace(0, 100*10**-6, 10000)
f = lambdify(t, Vt, "numpy")
pd.Series(f(ts), index=(ts*10**6)).plot(xlabel="us", ylabel="volts")
<AxesSubplot: xlabel='us', ylabel='volts'>