function [xt,xt_p,xt_h] = SDOF_ANA(t,m,k,xo,vo,Po)
%   Resolucion de la ecuacion para vibracion forzada: Armonica sin amortiguamiento
%   La carga esta representada por una funcion sinusoidak=l p(t)=Po.sin(w.t)
%   Los datos requeridos son: [xt] = SDOF_ANA(t,m,k,xo,vo,Po)
%   [t] → intervalo de tiempo
%   m → masa del sistema
%   k → rigidez del sistema
%   wn → frecuencia natural del sistema
%   xo → desplazamiento inicial
%   vo → velocidad inicial
%   xi → ratio de amortiguamiento
%   Po → clase de la carga
%   Po.f → magnitud de la carga
%   Po.w → frecuencia de la carga

wn=(k/m)^0.5;
beta=Po.w/wn;

A=xo;
B=(vo/wn)-(Po.f/k)*(beta/(1-beta^2));
C=(Po.f/k)*(1/(1-beta^2));

% Solucion homogenea

xt_p = A*cos(wn*t)+B*sin(wn*t);

% Solucion particular

xt_h = C*sin(Po.w*t);

% Solucion total

xt=xt_p+xt_h;

end