%El presente código tiene como finalidad obtener la respuesta al escalón
%unitario de un sistema con realimentación, conociendo la expresión de la
%función transferencia sin realimentación y las funciones de
%realimentación de acuerdo al esquema planteado en el ejercicio resuelto nº
%3 del capítulo 5

z=0;
while (z==0);
     z=1;
     a=input('Introduzca entre corchetes [] los coeficionetes del numerador de la función transferencia:');
     b=input('Introduzca entre corchetes [] los coeficionetes del denominador de la función transferencia:');
     f=input('La función transferencia introducida es (presione ENTER):');
     T=tf(a,b)
     R1=input('Introduzca el valor de realimentación 1:');
     R2=input('Introduzca el valor de realimentación 2:');
     T1=T/(1+T*R1);
     g=input('La función transferencia equivalente es (presione ENTER):');
     Tt=T1/(1+T1*R2)
     resp=input('La función transferencia total es la deseada?[y/n]','s');
        if isempty(resp);
            resp = 'y';
        end
        if resp=='y';
        tiem=input('Ingrese el tiempo de simulación en segundos:');
        g=input('Presione la tecla ENTER para obtener la grafica de la respuesta','s');
        step(Tt,tiem);
        grid;
        resp=input('Desea realizar una nueva simulación?[y/n]:','s');
        if isempty(resp);
            resp = 'y';
        end
        if resp=='y';
            z=0;
            end
        else
        z=0;    
        end 
end