You can first substitute $u=x^2$, use a useful property of the Laplace Transform and the Beta function and then utilize the reflection formula for the Gamma function
\begin{align}
\int_0^\infty \sin(x^2)\,dx &= \frac{1}{2}\int_0^\infty \frac{\sin u}{\sqrt{u}}\,du \\
&=\frac{1}{2\sqrt{\pi}} \int_0^\infty \frac{ds}{\sqrt{s}(s^2+1)}\, \\
&=\frac{1}{2\sqrt{\pi}} \int_0^\infty \frac{1}{\nu^{1/4}(1+\nu)}\frac{d\nu}{2\sqrt{\nu}} \\
&=\frac{1}{4\sqrt{\pi}} \int_0^\infty \frac{d\nu}{\nu^{3/4}(1+\nu)} \\
&=\frac{1}{4\sqrt{\pi}} B(1/4,3/4) \\
&=\frac{1}{2\sqrt{\pi}} \frac{\Gamma(1/4)\Gamma(3/4)}{\Gamma(1)} \\
&=\frac{1}{2\sqrt{\pi}} \frac{\pi}{\frac{1}{\sqrt{2}}} \\
&=\sqrt{\frac{\pi}{8}}
\end{align}
My approach:
Let
$$I(t) = \int_{0}^{\infty} \frac{xt - \sin(xt)}{x^3\left(x^2 + 4\right)} \:dx$$
Where $I = I(1)$
Taking the first derivative:
$$ \frac{dI}{dt} = \int_{0}^{\infty} \frac{x - x\cos(xt)}{x^3\left(x^2 + 4\right)} \:dx = \int_{0}^{\infty} \frac{1 - \cos(xt)}{x^2\left(x^2 + 4\right)} \:dx$$
Taking the second derivative:
$$ \frac{d^2I}{dt^2} = \int_{0}^{\infty} \frac{x\sin(xt)}{x^2\left(x^2 + 4\right)} \:dx = \int_{0}^{\infty} \frac{\sin(xt)}{x\left(x^2 + 4\right)} \:dx$$
Now, take the Laplace Transform w.r.t $t$:
\begin{align}
\mathscr{L}\left[ \frac{d^2I}{dt^2}\right] &= \int_{0}^{\infty} \frac{\mathscr{L}\left[\sin(xt)\right]}{x\left(x^2 + 4\right)} \:dx \\
&= \int_{0}^{\infty} \frac{x}{\left(s^2 + x^2\right)x\left(x^2 + 4\right)}\:dx \\
&= \int_{0}^{\infty} \frac{1}{\left(s^2 + x^2\right)\left(x^2 + 4\right)}\:dx
\end{align}
Applying the Partial Fraction Decomposition we may find the integral
\begin{align}
\mathscr{L}\left[ \frac{d^2I}{dt^2}\right] &= \int_{0}^{\infty} \frac{1}{\left(s^2 + x^2\right)\left(x^2 + 4\right)}\:dx \\
&= \frac{1}{s^2 - 4} \int_{0}^{\infty} \left[\frac{1}{x^2 + 4} - \frac{1}{x^2 + s^2} \right]\:dx \\
&= \frac{1}{s^2 - 4} \left[\frac{1}{2}\arctan\left(\frac{x}{2}\right) - \frac{1}{s}\arctan\left(\frac{x}{s}\right)\right]_{0}^{\infty} \\
&= \frac{1}{s^2 - 4} \left[\frac{1}{2}\frac{\pi}{2} - \frac{1}{s}\frac{\pi}{2} \right] \\
&= \frac{\pi}{4s\left(s + 2\right)}
\end{align}
We now take the inverse Laplace Transform:
$$ \frac{d^2I}{dt^2} = \mathscr{L}^{-1}\left[\frac{\pi}{4s\left(s + 2\right)} \right] = \frac{\pi}{8}\left(1 - e^{-2t} \right) $$
We now integrate with respect to $t$:
$$ \frac{dI}{dt} = \int \frac{\pi}{8}\left(1 - e^{-2t} \right)\:dt = \frac{\pi}{8}\left(t + \frac{e^{-2t}}{2} \right) + C_1$$
Now
$$ \frac{dI}{dt}(0) = \int_{0}^{\infty} \frac{1 - \cos(x\cdot 0)}{x^2\left(x^2 + 4\right)} \:dx = 0 = \frac{\pi}{8}\left(\frac{1}{2} \right) + C_1 \rightarrow C_1 = -\frac{\pi}{16}$$
Thus,
$$ \frac{dI}{dt} = \int \frac{\pi}{8}\left(1 - e^{-2t} \right)\:dt = \frac{\pi}{8}\left(t + \frac{e^{-2t}}{2} \right) - \frac{\pi}{16}$$
We now integrate again w.r.t $t$
$$ I(t) = \int \left[\frac{\pi}{8}\left(t + \frac{e^{-2t}}{2} \right) - \frac{\pi}{16} \right] \:dt = \frac{\pi}{8}\left(\frac{t^2}{2} - \frac{e^{-2t}}{4} \right) - \frac{\pi}{16}t + C_2 $$
Now
$$I(0) = \int_{0}^{\infty} \frac{x\cdot0 - \sin(x\cdot0)}{x^3\left(x^2 + 4\right)} \:dx = 0 = \frac{\pi}{8}\left( -\frac{1}{4} \right) + C_2 \rightarrow C_2 = \frac{\pi}{32}$$
And so we arrive at our expression for $I(t)$
$$I(t)= \frac{\pi}{8}\left(\frac{t^2}{2} - \frac{e^{-2t}}{4} \right) - \frac{\pi}{16}t + \frac{\pi}{32}$$
Thus,
$$I = I(1) = \frac{\pi}{8}\left(\frac{1}{2} - \frac{e^{-2}}{4} \right) - \frac{\pi}{16} + \frac{\pi}{32} = \frac{\pi}{32}\left(1 - e^{-2}\right)$$
Best Answer
The method I took was:
Let
$$ I(t) = \int_{0}^{\infty} \frac{\sin(kxt)}{x\left(x^2 + 1\right)} \:dx$$
Take the Laplace Transform
\begin{align} \mathscr{L} \left[I(t) \right]&= \int_{0}^{\infty} \frac{\mathscr{L} \left[\sin(kxt)\right]}{x\left(x^2 + 1\right)} \:dx \\ &= \int_{0}^{\infty} \frac{kx}{\left(k^2x^2 + s^2\right)x\left(x^2 + 1\right)} \:dx \\ &= \int_{0}^{\infty} \frac{k}{\left(k^2x^2 + s^2\right)\left(x^2 + 1\right)} \:dx \\ &= \int_{0}^{\infty} \left[\frac{k^3}{\left(k^2 - s^2\right)\left(k^2x^2 + s^2 \right)} - \frac{k}{\left(k^2 - s^2\right)\left(x^2 + 1\right)}\right] \:dx \\ &= \left[\frac{k^3}{\left(k^2 - s^2\right)}\frac{\arctan\left(kx\right)}{ks} - \frac{k}{\left(k^2 - s^2\right)}\arctan(x)\right]_{0}^{\infty} \\ &= \frac{k^2}{s\left(k^2 - s^2\right)}\frac{\pi}{2} - \frac{k}{\left(k^2 - s^2\right)}\frac{\pi}{2} \\ &= \frac{k}{k^2 - s^2}\left[ \frac{k}{s} - 1 \right]\frac{\pi}{2} \\ &= \frac{k}{s\left(k + s\right)}\frac{\pi}{2} \end{align}
And thus,
$$ I(t) = \mathscr{L}^{-1}\left[\frac{k}{s\left(k + s\right)}\frac{\pi}{2} \right] = \left[1 - e^{-kt} \right]\frac{\pi}{2}$$
Lastly,
$$ I = I(1) = \left[1 - e^{-k} \right]\frac{\pi}{2}$$