code: def compute_resids(xy,uv,fit):
      
      print('FIT coeffs: ',fit['coeffs'])
      xn,yn = apply_fit(uv,fit['coeffs'])
      resids = xy - np.transpose([xn,yn])
      return resids
docs: Compute the residuals based on fit and input arrays to the fit