For comparison, do harmonic sum with NumPy