clc; clear all; close all; 
pkg load control;
% Sistema de segundo orden
s = tf('s');

k = 1 ; % ganancia (adimensional en este caso)
tau1 = 3; tau2 = 15; % constante de tiempo (min)
taun = input('tau n =  ');

g = k*(taun*s+1)/(tau1*s+1)/(tau2*s+1);

t = 0:0.1:100; Du = 1 ; % altura del escalón
y=Du*step(g,t);

plot(t,y)
title('sistema de segundo orden con dinámica en numerador- respuesta a un escalón de altura 1')
xlabel('t')
ylabel('y')

pause
y = impulse(g,t);
plot(t,y)
title('sistema de segundo orden con dinámica en numerador - respuesta a un pulso unitario')
xlabel('t')
ylabel('y')