% solução numérica de sistema com 1 gdl
% amortecimento quadrático, excitações harmônica e impulsiva
clc
clearvars
close all

% ts=[0 100];
ts=linspace(0,100,500);
% ci=[-0.01 0];
ci=[0 0];

[t,x]=ode45('aqei',ts,ci);

plot(t,x(:,1))
ylabel('deslocamento (m)');
xlabel('tempo (s)')
grid on