import numpy as np
import control

# %% Creating a state space model:
A = np.array([[0, 1], [2, 3]])
B = np.array([[4], [5]])
C = np.array([[6, 7]])
D = np.array([[8]])

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

# %% Getting the model matrices as lists and then as arrays:
(A_list, B_list, C_list, D_list) = control.ssdata(S)
A_array = np.array(A_list)
B_array = np.array(B_list)
C_array = np.array(C_list)
D_array = np.array(D_list)

# %% Displaying the matrices as arrays:
print('A_array =', A_array)
print('B_array =', B_array)
print('C_array =', C_array)
print('D_array =', D_array)