import numpy as np
import control

m = 10  # [kg] 
k = 4  # [N/m] 
d = 2  # [N/(m/s)] 

# %% System matrices as 2D arrays:
A = np.array([[0, 1], [-k/m, -d/m]])
B = np.array([[0], [1/m]])
C = np.array([[1, 0]])
D = np.array([[0]])

# %% Creating the state space model:
S = control.ss(A, B, C, D)

# %% Displaying S: 
print('S =', S)