oops realized huge bug - doesnt calculate phase for each point, just averages the phase at the center of the fit so output doesnt make any sense