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

ts=[0 40];
% ci=[-0.01 0];
ci=[0 0];

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

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