I'm doing a curious exercise. I have to find a method to interpolate the following ballistic trajectory in 3D space:
$$\left\{\begin{matrix} x(t) = \frac{v\cos(\phi)\cos(\beta)}{k}(1-e^{-kt})+d_{0}\cos(\alpha)\\ y(t) = \frac{v\sin(\phi)\cos(\beta)}{k}(1-e^{-kt})+d_{0}\sin(\alpha) \\ z(t) = \left ( \frac{v\sin(\beta)}{k}+\frac{g}{k^2}\right)(1-e^{-kt})-\frac{g}{k}t \end{matrix}\right.$$
provided a set of data $(x_n,y_n,z_n)$. Therefore, to do this I think that I should eliminate the dependence on the parameter 't' and write the curve in the form of two different functions of one variable: $z(x)$ and $z(y)$:
$$z(x)= \left( v\sin(\beta)+\frac{g}{k}\right )\frac{x-d_0\cos(\alpha)}{v\cos(\phi)\cos(\beta)}+\frac{g}{k^2}\log\left(1-k\frac{x-d_0\cos(\alpha)}{v\cos(\phi)\cos(\beta)}\right) $$
$$z(y)= \left( v\sin(\beta)+\frac{g}{k}\right )\frac{y-d_0\sin(\alpha)}{v\sin(\phi)\cos(\beta)}+\frac{g}{k^2}\log\left(1-k\frac{y-d_0\sin(\alpha)}{v\sin(\phi)\cos(\beta)}\right) $$
Is this actually the best way to set the problem before interpolating?