void VectNormalize(Vect v, Vect v_prim) { double mag; mag=VectMag(v); v_prim[X]=v[X]/mag; v_prim[Y]=v[Y]/mag; v_prim[Z]=v[Z]/mag; }
double VectUnit(Vect v) { double mag; mag = VectMag(v); v[X] /= mag; v[Y] /= mag; v[Z] /= mag; return mag; }