In [1]:
import numpy as np
import sympy as sy

In [2]:
kv, Tv = sy.symbols("kv,Tv", real=True, positive=True)

In [3]:
s = sy.symbols('s', real=False)

In [6]:
sy.apart(kv/(s*s*(Tv*s+1)), s)


Out[6]:
Tv**2*kv/(Tv*s + 1) - Tv*kv/s + kv/s**2

In [ ]: