When you integrate the changing speed or the constant speed with the variable time step, you add up all these small, partial sections. Space is not made of sections of different lengths (that expanded at different rates) unlike the path calculated by the numerical integration of a motion of a body moving with variable speed. Can you see the difference? It's crucial.
I introduce parameter b(t)=1/a(t) and have the exact, same problem.
It occurred to me, that this integration would be suitable for the path of gravitationally redshifted photon, because gravitationally curved spacetime retains its distributions of time dilation and length contraction factors at all times at all distances:
This is not the case with expanding, intergalactic space.