clc; clear all; close all; 
pkg load control;
% Sistema de segundo orden respuesta a la frecuencia
s = tf('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);

t = 0:0.1:100;
A = 1; % amplitud de la seņal de entrada
w = input('frecuencia =  ')
u = A*sin(w*t);
y = lsim(g,u,t);

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