A = [-1, 10, 3; 1, 5, 3; 4, 13, 7] eig(A) B = [1; 1; 1] C = eye(3) P = ss(A, B, C, 0) x0 = B initial(P, x0) eig(A) help initlal help initial [T, D] = eig(A) T\A*T A = [5, -4; 6, -5] B = [1; 1] p = [-3; -4] F = place(A, B, p) A = [1, 2; -3, -4] B = [2; 3] p = [-2; -3] F = place(A, B, p) Uc = ctrb(A, B) det(Uc)