import control

# %% Creating the s-transfer function:

s = control.tf('s')
H_cont = 3/(2*s + 1)

# %% Discretizing:

Ts = 0.1  # Time step [s]
H_disc = control.sample_system(H_cont, Ts, method='zoh')

print('H_disc(z) =', H_disc)