Dust tail to be represented as a quadratic equation (vertex points along the curve) Then tail will be planar. Dust velocity equation will act as axis of symmetry, axisDV

1. First get quadratic equation kx**2, where k is calculated from dust ejection velocity, and magnitude from heliocentric distance (k is determined by comparing those two values of many comets, via testing)

2. Get a list of vertex points, vertexArray

3. Since that equation's axis of symmetry will be x=0, we can get the rotation to aDV, A

4. Next, get the translation to the point in space, B

5. Combine these to the final transformation matrix T(x) = B(A(x))

6. Put all points in vertexArray through T

**Finding k**, which represents the magnitude of the dust tail, from comparing the usual case of dust ejection velocity (which takes into account heliocentric distance) vector length.

Min dust ejection velocity vector length: 5.2853

Max DEV length: 23.2426

Usual DEV length: (10, 15)

**k = DEV_length * factor // TODO: need factor? Can only determine after drawing**

Still need to find projection from 3-D coordinates to 2-D window (screen).

**Comet data and parse code below**