close all;
clc;
a=[0 1;-25 -4];
b=[0;25];
c=[1 0];
d=[0];
td=0.2;
[num den]=ss2tf(a,b,c,d);
system=tf(num,den,'inputdelay',td);
sys=pade(system,2);
[GM PM Wgc Wpc]= margin(sys)
figure(1);
margin(sys);
grid;
figure(2);
rlocus(sys);
grid;
figure(3);
nyquist(sys);
grid;
figure(4);
sys=feedback(sys,1);
step(sys,'-');
grid;
if(PM>0 && GM>0)
fprintf('System is stable');
else
fprintf('System is not stable')
end