// Sistema de segundo orden respuesta a la frecuencia
s = poly(0,'s')

k = 1 // ganancia (adimensional en este caso)
tau = 5 // constante de tiempo (min)
z = input('factor de dumping =  ')

g = k/(tau^2*s^2+2*z*tau*s+1)

sl = syslin('c',g)

t = 0:0.1:100
A = 1 // amplitud de la señal de entrada
w = input('frecuencia =  ')
u = A*sin(w*t)
y = csim(u,t,sl)

clf
plot2d(t,u)
title('sistema de segundo orden - respuesta a la frecuencia')
xlabel('t')
ylabel('y')

plot2d(t,y,style =[color("blue")])
