% Determinar la respuesta en el tiempo del ejercicio marcado en la figura

clear all
clc

% Datos
g=981; %[cm/s^2]
L=60; %[cm]

% Determinamos condiciones de frontera
theta_o=[5 20 40 60];
vo=0;

% Calculo de la frecuencia natural
wn=(g/L)^0.5 %[1/s]\\
Tn=2*pi()/wn %[s]

% Determinamos el tiempo de la respuesta
t=[0:0.01:4];

% Calculo de la respuesta en el tiempo
    % Resolucion de la ecuacion para vibracion libre sin amortiguamiento
    %   Los datos requeridos son: [xt,vt,at] = SDOF_VL(t,wn,xo,vo)
    %   [t] → intervalo de tiempo
    %   wn → frecuencia natural del sistema
    %   xo → desplazamiento inicial
    %   vo → velocidad inicial
for i=1:length(theta_o)
    [theta{i}] = SDOF_VL(t,wn,theta_o(1,i),vo);
end

figure
for i=1:length(theta_o)
plot(t,theta{i})
hold on
end
xlabel('Tiempo [s]')
ylabel('Angulo')
title('Respuesta en el tiempo')