Animações em Física com gnuplot: Lançamento de um foguete inicialmente acelerado

Caros leitores, uma pequena animação da cinemática do lançamento de um foguete inicialmente acelerado e os correspondentes gráficos das velocidades Vx(t), Vy(t) e V(t) é aqui construída segundo os comandos básicos do aplicativo gnuplot.  Os objetos presentes na animação (rectangle, arrow, circle e label) são implementados facilmente. As equações físicas envolvidas na animação são,

1. Primeiro Estágio: Acelerado por um intervalo de tempo, $\Delta t = 30\, s$, com aceleração de $a = 46\, m/s^{2}$, sob o ângulo de inclinação $\theta = 70^{o}$:

1.1 As funções horárias do movimento: $$X_{1}(t)=x_{o}+v_{ox}t + \frac{1}{2} a_{x} t^{2} \qquad \mbox{e} \qquad Y_{1}(t)=y_{o}+v_{oy}t + \frac{1}{2} a_{y} t^{2}$$
onde $v_{ox}=v_{o}\cos{\theta}$, $v_{oy}=v_{o}\sin{\theta}$,  $a_{x}=a \cos{\theta}$ e $a_{y}=a \sin{\theta}$.

1.2. As funções da velocidade: $$V_{x}(t)=v_{ox} + a_{x} t  \qquad \mbox{e} \qquad V_{y}(t)=v_{oy} + a_{y} t. $$  O módulo de $V(t)$ para qualquer instante $t$ é determinado por $$V(t) = \sqrt{(V_{x}(t))^{2}+(V_{y}(t))^{2}}.$$

2. Segundo Estágio: Para os instantes de tempo $t > T = 30 \,s$, agora, sob a ação da gravidade apenas, temos:

2.1 As funções horárias do movimento: $$X_{2}(t)=X_{1}(T)+v_{1x}(t-T)  \qquad \mbox{e} \qquad Y_{2}(t)=Y_{1}(T)+v_{1y}(t-T) + \frac{1}{2} g (t-T)^{2}$$
onde $v_{1x}=V(T) \cos{\theta}$, e $v_{1y}=V(T) \sin{\theta}$.

2.2. As funções da velocidade: $$V_{2x}(t)=v_{1x}  \qquad \mbox{e} \qquad V_{2y}(t)=v_{1y} - g (t-T)$$.

No aplicativo do gnuplot é possível alterar todos os parâmetros de entrada facilmente, ou seja, é bem geral.