#!/usr/bin/env python
# coding: utf-8

# In[1]:


from sympy import *


# In[2]:


s, Vpk, C1, C2, C3, R1, R2 = symbols('s Vpk C1 C2 C3 R1 R2')
init_printing(use_unicode=True)

zC1:1/(s·C1);
zC2:1/(s·C2);
zC3:1/(s·C3);
z22:R2+zC2;
z223:zC3·z22/(zC3+z22);
z2231:zC1+z223;
VM1:V1·z2231/(R1+z2231);
# In[3]:


zC1 = 1 / (s*C1)
zC2 = 1 / (s*C2)
zC3 = 1 / (s*C3)
z22 = R2 + zC2
z223 = zC3*z22/(zC3+z22)
z2231 = zC1+z223
V1 = Vpk / s
# V1 = symbols('V1')
VM1 = V1 * z2231/(R1+z2231)
for e in [zC1, zC2, zC3, z22, z223, z2231, VM1]:
    pprint(e)
    print('***')


# In[4]:


C = Symbol('C', real=True)
R = Symbol('R', real=True)
VM1 = VM1.subs([(c, C) for c in [C1, C2, C3]])
VM1 = VM1.subs([(r, R) for r in [R1, R2]])
VM1 = VM1.simplify()
VM1


# In[7]:


t = symbols('t')
VM1t = inverse_laplace_transform(VM1, s, t, noconds=True)
simplify(VM1t)