function [xt] = SDOF_IStepF(t,m,k,Po,tr)
% Resolucion de la ecuacion para vibracion forzada para una carga
% incrementales sin amortiguamiento
%   Los datos requeridos son: [xt] = SDOF_IStepF(t,m,k,Po,tr)
%   [t] → intervalo de tiempo
%   m → masa del sistema
%   k → rigidez del sistema
%   xt → respuesta de desplazamiento en el tiempo sin amortiguamiento
%   tr → tiempo de aplicacion de la carga en la fase incremental

wn = (k/m)^0.5;

xt=zeros(1,length(t));
for i=1:length(t)
    if t(i)<=tr
    xt(i) = (Po/k)*((t(i)/tr)-((sin(wn*t(i)))/(wn*tr)));
    else
    xt(i)=(Po/k)*(1-(1/(wn*tr))*(sin(wn*t(i))-sin(wn*(t(i)-tr))));
    end

end

end