Solution (source code)

= Solution

To first order, the radial and transverse comoving separations are
$$
\Delta\chi_\parallel=\frac{\Delta z}{H(z)}
=\frac{f(z,\Omega_{m,0})}{H_0}\Delta z,
\qquad
\Delta\chi_\perp=\chi(z)\Delta\theta
=\frac{I(z;\Omega_{m,0})}{H_0}\Delta\theta,
$$
where
$$
f(z,\Omega_{m,0})=
\frac1{\sqrt{\Omega_{m,0}(1+z)^3+1-\Omega_{m,0}}}.
$$
The local spatial geometry is that of <Euclidean space>, so the <Pythagorean theorem> gives
$$
\Delta\chi
=\frac{f}{H_0}
\sqrt{(\Delta z)^2+\left(\frac I f\right)^2(\Delta\theta)^2}.
$$
Consequently the requested function is
$$
\boxed{y(z,\Omega_{m,0})=
\frac{I(z;\Omega_{m,0})}{f(z,\Omega_{m,0})}
=H(z)\chi(z)}.
$$