A = [0 1 0 0; 0 0 0 0.0024; 0 0 0 1; 0 -0.0024 0 0]; B = [0 0; 1 0; 0 0; 0 1]; C = [1 0 0 0; 0 0 1 0]; co = ctrb(A, B) num = rank(co); if num == 4 'The system is controllable.' else 'The system is not controllable.' end ob = obsv(A, C) num = rank(ob); if num == 4 'The system is observable.' else 'The system is not observable.' end