%% Aula 10 - Exemplo de filme: queda livre
clear; 
% Inicializacao
tempo0=0;
tempoF=4.5;
Deltat=0.1;
tempoArray=tempo0:Deltat:tempoF;
yarray=tempoArray;
Ntempos=length(tempoArray);
y0=100;    % altura inicial
yarray(1)=y0; % y em t=0
g=9.8;     % aceleração da gravidade
% Ajusta os parâmetros da figura: fundo branco, tamanho da janela
figure1 = figure(1);
set(figure1,'color','white');
winsize = get(figure1,'Position');
winsize(1:2) = [0 0];
clear M;
count=1;
p1=plot(0,yarray(1),'bo'); % plota um ponto
ylim([0 100]);             % mantem a escala do gráfico constante
M(count)=getframe(figure1,winsize); % armazena a frame em "M"
delete(p1);   % deleta o gráfico (já armazenou em "M")
count=count+1;  % atualiza o contador    
%% Loop nos tempos
for n=2:Ntempos
    yarray(n)=y0-0.5*g*(tempoArray(n))^2; % atualiza y
    p1=plot(0,yarray(n),'bo'); % plota um ponto
    ylim([0 100]);             % mantem a escala do gráfico constante
    M(count)=getframe(figure1,winsize); % armazena a frame em "M"
    delete(p1);               % deleta o gráfico (já armazenou em "M")
    count=count+1;            % atualiza o contador 
end
movie2avi(M,'QuedaLivre.avi'); % gera o avi 
% (É possível gerar em outros formatos)