1499/CH7/EX7.33/q33.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

A=[-2 0;1 -1];
B=[0;1];

[r c]=size(A)
I=eye(r,c)
P=cont_mat(A,B);
disp(P,"Controllability Matrix=")
d=det(P)
if d==0
  printf("matrix is singular, so system is uncontrollable");
else
  printf("system is controllable");
end;