% solução numérica de sistema com 1 gdl sob excitação harmônica
clc; clearvars; close all

ts=[0 15];
ci=[0.01 0.1];

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

plot(t,x(:,1),'r-',t,x(:,2),'b-')
xlabel('tempo (s)')
ylabel('deslocamento(m), velocidade(m/s)')
legend('deslocamento','velocidade')
grid on