We're given the ODE $$2(x^2+1)y'(x) + (x-i)y(x)=0$$
Deriving $f$ in respect to $x$, we get the equality $$f'(x)=\frac{-1}{2(x+i)}f(x)$$ which lines up nicely as a solution to the above ODE.
All the solutions to this ODE hold the form $$y(x)=\lambda e^{A(t)}$$
with $\lambda\in\mathbb{R}$ fixed and $A(t)$ a fixed anti-derivative of, in this case, $x\mapsto\frac{-1}{2(x+i)}$, so we choose $A(t)=\frac{-1}{2}\ln|x+i|$.
This leaves us with
$$f(x)=\lambda e^{\frac{-1}{2}\ln|x+i|}$$
The problem is that it's nowhere near the expected answer of
$$f(x)=\sqrt{\pi} (x^2+1)^\frac{-1}{4}e^{\frac{i}{2}\arctan{x}}$$
I can also try to directly compute $f(x)$ since the paper suggests to use the Gaussian integral. A quick variable change of $t=u^2$ leaves us with
$$f(x)=2 \int_0^{+\infty} e^{-u^2} e^{i u^2 x} du$$
But I can't get through this integral in order to get $\int_0^{+\infty} e^{-u^2} du$ to appear, can you provide some hints?