clear all;
close all;

k = 3;
t = 2;
G = tf ( k , [t,1]);
C = k
s=tf('s');
hold on

%opt = stepDataOptions('StepAmplitude',1);
%step(G,opt);

%opt = stepDataOptions('StepAmplitude',2);
%step(G,opt);

%opt = stepDataOptions('StepAmplitude',1);
%step(G,opt);

H=(k/(k+1))/(1+s*t/(k+1)); %Y/Yc
%H=(C*G)/(1+C*G);
%H=feedback(C*G+1);
k=1
H1=(k/(k+1))/(1+s*t/(k+1));
step(H1)
k=2
H2=(k/(k+1))/(1+s*t/(k+1));
step(H2)
k=10
H3=(k/(k+1))/(1+s*t/(k+1));
step(H3);

%------------------

clear all;
t=2;
k=10;
G = tf ( k , [t,1]);
C = k;
U=C/(1+C*G);
Huyc=U/1;
figure
hold on;
step(Huyc);

epsilon=1/(1+C*G);
Hepsilonyc=epsilon/1;
step(Hepsilonyc);