Q=100;
f_c = 0.843;
alpha = 2*pi*f_c;
betta = pi*f_c/Q;

G = [[-betta -alpha]; [alpha -betta]];
F = -G;



%G=[[-3 1 0]; [0 -2 1]; [0 0 -1]]
%F=[0; 0; 1]

C=[[1 0];[0 1]]


sys=ss(G,F,C,0);
step(sys)

impulse(sys)



null(G)

det(G)

p=charpoly(G)
r=roots(p)


syms t; 
expm(t*G)