# A Gift That Keeps On Giving

We see from “Seek-Lock-Strike!” Again that given the missile’s position $(x, y),$ $\frac{dx}{dy}=\frac{a+v_a t-x}{b-y}$

where $x$ and $y$ are themselves functions of time $t.$

It means $\frac{dx}{dy} = \frac{\frac{dx}{dt}}{\frac{dy}{dt}}=\frac{a+v_a t-x}{b-y}\implies \frac{dx}{dt} = \frac{a+v_a t-x}{b-y}\cdot\frac{dy}{dt}.$

That is, let $\kappa(x,y,t) = \frac{a+v_a t-x}{b-y},$ $\frac{dx}{dt}=\kappa\cdot \frac{dy}{dt}.\quad\quad\quad(1)$

We also have (see “Seek-Lock-Strike!”) $v_m = \sqrt{(\frac{dx}{dt})^2+(\frac{dy}{dt})^2}\implies v_m = \sqrt{1+\kappa^2}|\frac{dy}{dt}| .$

Since $\frac{dy}{dt} >0,$ $\frac{dy}{dt} = \frac{v_m}{\sqrt{1+\kappa^2}}.\quad\quad\quad(2)$

Substitute (2) into (1) yields $\frac{dx}{dt} = \frac{v_m\cdot \kappa}{\sqrt{1+\kappa^2}}.\quad\quad\quad(3)$

It follows that $(x(t), y(t))$, the position of the missile satisfies the initial-value problem $\begin{cases} \frac{dx}{dt} = \frac{\kappa\cdot v_m}{\sqrt{1+\kappa^2}} \\ \frac{dy}{dt} = \frac{v_m}{\sqrt{1+\kappa^2}} \\x(0)=0, y(0)=0\end{cases}\quad\quad\quad(4)$

To obtain the missile’s trajectory, we solve (4) numerically using the Runge-Kutta algorithm. It integrates (4) from $t=0$ to $t= t_*$ (see “Seek-Lock-Strike!”).

Fig. 1 $a=100\;m, b=3000\;m, v_a=1500\;ms^{-1},v_m=2000\;ms^{-1}$

The missile strike is illustrated in Fig. 1 and 2.

Fig. 2 $a=100\;m, b=3000\;m, v_a=1500\;ms^{-1},v_m=2000\;ms^{-1}$

Fig. 3 $a=-1200\;m, b=3000\;m, v_a=1500\;ms^{-1},v_m=2000\;ms^{-1}$

The trajectories shown are much smoother than those in “Seek-Lock-Strike!” Animated.

# “Seek-Lock-Strike!” Animated

In “Seek-Lock-Strike!” Again, we obtained the missile’s trajectory. Namely, $x = \frac{1}{2} \left(\frac{(b-y)^{r+1}}{b^r \cdot f \cdot (r+1)}-\frac{b^r \cdot f \cdot (b-y)^{1-r}}{1-r}\right) -k\quad\quad\quad(*)$

where $f = \frac{a}{b}+\sqrt{1+(\frac{a}{b})^2},$ $k = \frac{b}{2}\left(\frac{1}{f \cdot (r+1)}-\frac{f}{1-r}\right).$

Since the fighter jet maintains its altitude ( $y= b$), the missile must strike it at $(x_*, b)$. Setting $y=b$ in $(*)$ gives $x_* = -k.$

Hence, we can plot $(*):$

Fig. 1 $a = 100\;m, b=3000\;m, v_a=1500\;ms^{-1}, v_m=2000\;ms^{-1}$

We can also illustrate “Seek-Lock-Strike” in an animation:

Fig. 2 $a = 100\;m, b=3000\;m, v_a=1500\;ms^{-1}, v_m=2000\;ms^{-1}$

Fig. 3 $a = -2500\;m, b=3000\;m, v_a=1500\;ms^{-1}, v_m=2000\;ms^{-1}$