function [xt_a,xt_sa] = SDOF_StepF(t,m,k,xi,Po)
% Resolucion de la ecuacion para vibracion forzada para una carga
% escalonada (con y sin amortiguamiento) - Step function
%   Los datos requeridos son: [xt_a,xt_na] = SDOF_StepF(t,m,k,xi,Po)
%   [t] → intervalo de tiempo
%   m → masa del sistema
%   k → rigidez del sistema
%   xt_a → respuesta de desplazamiento en el tiempo con amortiguamiento
%   xt_sa → respuesta de sesplazamiento en el tiempo sin amortiguamiento

wn = (k/m)^0.5;
wd=wn*(1-xi^2)^0.5;

xt_sa = (Po/k)*(1-cos(wn*t));
xt_a = (Po/k)*(1-(exp(-xi*wn*t)).*((cos(wd*t))+(xi/(1-xi)^0.5)*(sin(wd*t))));

end