import numpy as np
import control

# %% Model parameters:

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 and printing the state space model:

S = control.ss(A, B, C, D)

print('S =', S)