f=[0.8760 0 0;0.2546 0.6621 -0.5701;0.1508 0.4221 1];
g=[0.2105;0.1033;0.1768];
c=[0 1 3.5];
d=0;
n=3;
CAM=ctrb(f,g);
rcam=rank(CAM);
if rcam==n
    disp('System is controlled')
elseif rcam<n
    disp('System is no controlled')