k1 = 6;
k2 = 13;
t1 = 0.45;
t2 = 0.95;
ksi1=0.45;
ksi2=0.95;

%g1 = tf([k1*t1 k1],[0.01 1]);      % w(s) = 6/s
%g2 = tf([k2*t1 k2],[0.01 1]);
%g3 = tf([k1*t2 k1],[0.01 1]);

%step(g1, 'r', g2, 'g', 0.1);% h(t)
ylabel('h(t)')
title('')
grid on
%bode(g1, 'r', g2, 'g', {0.01, 1000})    % AFC, PFC
title('')
grid on


s=tf('s');
%G = k1*s/(0.01*s+1);
%G2 = k2*s/(0.01*s+1);
%G = k1*exp(-0.3*s);
%G2 = k2*exp(-0.3*s);
%T = feedback(ss(G),.5);
%step(G, 'b', G2, 'g');
%bode(G, 'r', G2, 'b', {1, 10000000})
G = 400/(s*(0.0143*s+1)*(0.005*s+1)*(0.0125*s+1));
%step(G, 'b');
bode(G, 'r')
grid on