Update do_phase_fit.c

Correcting a coding error by inserting a number of good lags, nphi, in the expression for the velocity variance (for more detail see Task #1 on Pasha's website).