% Equaçao do movimento do exercício 7.1
clear all 
close all
clc
%-----------------------------------------------------------------------
%% dados de entrada
% intervalo de tempo
%tf=input('intervalo de tempo (s): ');
tf=5;                                                                      % 25
tspan=[0; tf];    
% posição inicial
%y01=input('posição inicial (m): ');
y01=0;                                                                      % 50
% velocidade inicial
%y02=input('velocidade inicial (m/s): ');
y02=0;
y0=[y01;y02];       
%-----------------------------------------------------------------------
%% solução da equação do movimento
[t,y]=ode45('func_ex7_01_DINAMICA',tspan,y0);
q=y(:,1);
qp=y(:,2);
%-----------------------------------------------------------------------
%% gráfiicos da posição e velocidade de q
figure('Name',['    t = ',num2str(tf),' s      q_o = ',...
    num2str(y01),'m      qp_o = ',num2str(y02),'m/s'])
subplot(211), plot(t,q,'k');
xlabel('seg')
ylabel('q [m]')
title('Posiçao')
grid on
subplot(212), plot(t,qp,'k');
xlabel('seg')
ylabel('qp [rad/s]')
title('Velocidade angular')
grid on
%%
figure('Name',['    t = ',num2str(tf),' s      q_o = ',...
    num2str(y01),'m      qp_o = ',num2str(y02),'m/s'])
plot(q,qp,'k');
xlabel('q [m]')
ylabel('qp [m/s]')
title('Velocidade')
grid on
%% ---------------------------------------------------------------------
nomearq = 'dinamica.mat';
save(nomearq, 'y')
%-----------------------------------------------------------------------
%% FIM
%-----------------------------------------------------------------------

