code: def findbestparams_ho(xsamples):
      
      return np.abs(leastsq(deltaH_ho,np.array([10.,10.,10.]), Dfun = Jac_deltaH_ho, args=(xsamples,))[0])[:3]
docs: Minimize sum of square differences of H_sho-<H_sho> for timesamples