A hybrid asymptotic-numerical method for calculating drag coefficients in 2-D low Reynolds number flows