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

% limites de tempo
ts=[0 25]; % (lse,leh,lep,lei)
% ts=[0 5]; % (mqs)

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

% solução via Runge-Kutta
[t,x]=ode45('lse',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);
% c) ci: [0 0], excitação harmônica (leh);
% d) ci: [0 0], excitação periódica (lep);
% e) ci: [0 0], excitação impulsiva (lei);
% máquina de serraria (mqs).