
s=tf('s');

kp=1;
tp=0.5;

p=0;
for tc=2:0.01:5
theta = tc/tp;
p=p+1;
theta = tc/tp;
b(p,1)=tc;
b(p,2) = 0.31 + 0.6*(1/theta);

b(p,3) = tp*(0.777+0.45*theta);

b(p,4) = tp*(0.44-0.56*((1/theta)^(2.2)));

end


