% 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)
ts=[0 15];

% condições iniciais
% ci=[0 1.5];
% ci=[0 0];
ci=[1e-3 20e-3];

% 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 1.5], 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).