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

ts=[0 20]; % tempo
ci=[0.01 1]; % condições iniciais

% [t,x]=ode45('lse',ts,ci); % Runge-Kutta
% 
% plot(t,x(:,1),'r-',t,x(:,2),'b-')
% ylabel('variáveis de estado');
% xlabel('tempo (s)')
% grid on
% legend('deslocamento (m)','velocidade (m/s)')

% casos de interesse
% a) x0=0.01 m; xp0=0.1 m/s; f=0 (lse)
% b) x0=0 m; xp0=0 m/s; f=150sen(3t) N (leh)
% c) x0=0.01 m; xp0=0.1 m/s; f=150sen(6t) N (seh)
% d) x0=0 m; xp0=0 m/s; f=mc.g N (lec)
% e) x0=0 m; xp0=0 m/s; f=150.dl(t) N (lei)
% f) x0=0 m; xp0=0 m/s; f=1500sen[(omn/3)t] N (neh)
% g) x0=0.01 m; xp0=1 m/s; f=1500.dl(t) N (ner/ler)

% uso no caso (g), no lugar das linhas 9 a 15
% [t,x]=ode45('ner',ts,ci); % Runge-Kutta
% plot(t,x(:,1),'r-')
% hold on
% [t,x]=ode45('ler',ts,ci); % Runge-Kutta
% plot(t,x(:,1),'b-')
% hold off
% 
% ylabel('deslocamento (m)');
% xlabel('tempo (s)')
% grid on
% legend('caso não linear','caso linear')