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

ts=[0 30]; % tempo
ci=[0; -0.1; 0; 0]; % condições iniciais

[t,x]=ode45('lei2',ts,ci); % Runge-Kutta

subplot(2,1,1)
plot(t,x(:,1),'r-',t,x(:,2),'b-')
ylabel('deslocamentos (m)');
xlabel('tempo (s)')
grid on
legend('massa 1','massa 2')
subplot(2,1,2)
plot(t,x(:,3),'r-',t,x(:,4),'b-')
ylabel('velocidades (m/s)');
xlabel('tempo (s)')
grid on
legend('massa 1','massa 2')

% casos de interesse
% a) ci: [0; 0.1; 1; 0], excitação harmônica (leh2);
% b) ci: [0; -0.1; 0; 0], excitação tipo pulso (lei2)