// Tanque agitado V = 5m3, v = 1 m3/min 

s = poly(0,'s')

k = 1 // ganancia (adimensional en este caso)
tau = 5 // constante de tiempo (min)
g = k/(tau*s+1)

sl = syslin('c',g)

t = 0:0.1:20; Du = 0.5 // altura del escalón
y = Du*csim('step',t,sl)

clf
plot2d(t,y)
title('sistema de primer orden - respuesta a un escalón de altura 0.5')
xlabel('t')
ylabel('y')

pause
sl1 = tf2ss(sl); disp(sl1)

pause
y = Du*csim('step',t,sl1)
clf
plot2d(t,y)
title('idem pero con state-space')
xlabel('t')
ylabel('y')

pause
disp(ss2tf(sl1))
