This formula is more accurate at predicting SNe Ia comoving distances than FLRW and since it's analytical, it will execute upwards of 10,000 faster and give you machine precision.

$D_C(z) = \frac{t_o\left(2V_0 - A t_o\right) z}{2 + z}$

That's it. That's all there is to it. Here's how it compares to FLRW:

Red is the analytical formula, blue is FLRW. Using the full covariance matrix of the Pantheon+ SH0ES project, the reduced $\chi^2$ of 1.68 for the analytical formula is significantly better than 2.49 for the numerical FLRW solution.

This is an outrageous claim, so it should be easy to debunk. The full paper is here: Acceleration Law Article and the notebook is here: Acceleration Law Notebook. I'd appreciate your feedback, positive or negative.