%Sistema dinámico lineal de segundo orden
%Para su correcta ejecución requiere del
%programa cap2_simuejemplo27
%Entradas:  t  tiempo
%           x  vector de estados
% Salida  xp=Ax+Bu
% la matriz A está formada por  rho factor
% de amortiguamiento y omega_n la frecuencia  natural
% de resonancia mecánica.
function xp =cap2_ejemplo27(t,x)
% frecuencia natural de resonancia
wn=1; 

% factor de amortiguamiento
rho=0.4;
A=[ 0, 1; 
   -wn^2, -2*rho*wn];
%matriz B
B=[0;
   wn];

%señal de entrada
u=1;
%sistema dinámico lineal
xp=A*x+B*u;
end