% solução numérica de sistema com 1 gdl
clc
clearvars
close all

% limites de tempo
ts=[0 20];

% condições iniciais
% ci=[-0.01 0];
% ci=[0 0.01];
ci=[0 0];

% solução via Runge-Kutta
[t,x]=ode45('leh',ts,ci);

% visualização da solução (estado)
subplot(2,1,1)
plot(t,x(:,1),'b-')
ylabel('deslocamento (m)');
xlabel('tempo (s)')
grid on
subplot(2,1,2)
plot(t,x(:,2),'r-')
ylabel('velocidade (m/s)');
xlabel('tempo (s)')
grid on

% casos de interesse
% a) ci: [-0.01 0], vibração livre (lse);
% b) ci: [0, 0.01], vibração livre (lse)
% c) ci: [0 0], vibração forçada harmônica (leh);
