Please, note an error in checking for radial velocity when computing true anomaly. You cannot use the Vr formula I mentioned, as theta is not corrected yet! You just check the sign of Vr by performing the dot product between the given V and r vectors. See also how it is done in uncollected hw 3. Will be corrected at beginning of next lecture.