$\newcommand{\bbx}[1]{\,\bbox[15px,border:1px groove navy]{\displaystyle{#1}}\,}
\newcommand{\braces}[1]{\left\lbrace\,{#1}\,\right\rbrace}
\newcommand{\bracks}[1]{\left\lbrack\,{#1}\,\right\rbrack}
\newcommand{\dd}{\mathrm{d}}
\newcommand{\ds}[1]{\displaystyle{#1}}
\newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,}
\newcommand{\ic}{\mathrm{i}}
\newcommand{\mc}[1]{\mathcal{#1}}
\newcommand{\mrm}[1]{\mathrm{#1}}
\newcommand{\pars}[1]{\left(\,{#1}\,\right)}
\newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}}
\newcommand{\root}[2][]{\,\sqrt[#1]{\,{#2}\,}\,}
\newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}}
\newcommand{\verts}[1]{\left\vert\,{#1}\,\right\vert}$
\begin{align}
&\bbox[10px,#ffd]{\int_{0}^{\pi/2}{x^{3}\cos\pars{x} \over
3\sin\pars{x} - \sin\pars{3x}}\,\dd x} =
{1 \over 4}\int_{0}^{\pi/2}x^{3}\,{\cos\pars{x} \over
\sin^{3}\pars{x}}\,\dd x
\\[5mm] = &\
-\,{1 \over 8}\int_{x\ =\ 0}^{x\ =\ \pi/2}x^{3}\,\dd\bracks{1 \over \sin^{2}\pars{x}} =
-\,{\pi^{3} \over 64} + {3 \over 8}
\color{#00f}{\bf\int_{0}^{\pi/2}{x^{2} \over \sin^{2}\pars{x}}\,\dd x}
\label{1}\tag{1}
\end{align}
\begin{align}
&\color{#00f}{\bf\int_{0}^{\pi/2}{x^{2} \over \sin^{2}\pars{x}}
\,\dd x} =
\left.\Re\int_{x\ =\ 0}^{x\ =\ \pi/2}{\bracks{-\ic\ln\pars{z}}^{\, 2}
\over \bracks{\pars{z - 1/z}/\pars{2\ic}}^{\, 2}}
{\dd z \over \ic z}\,\right\vert_{\ z\ =\ \exp\pars{\ic x}}
\\[5mm] = &\
\left.4\,\Im\int_{x\ =\ 0}^{x\ =\ \pi/2}{z\ln^{2}\pars{z}
\over \pars{z^{2} - 1}^{\, 2}}\,\dd z\,\right\vert_{\ z\ =\ \exp\pars{\ic x}}
\\[5mm] = &\
-\,4\,\Im\int_{1}^{0}{\ic y\,\bracks{\ln\pars{y} + \ic\pi/2}^{\, 2}
\over \pars{\ic^{2}y^{2} - 1}^{\, 2}}\,\ic\,\dd y
\\[5mm] = &\
-\,4\pi\int_{0}^{1}{y\ln\pars{y}\over \pars{1 + y^{2}}^{\, 2}}\,\dd y =
-\,\pi\int_{0}^{1}{\ln\pars{y}\over \pars{1 + y}^{\, 2}}\,\dd y
= \pi\ln\pars{2}
\end{align}
(\ref{1}) becomes
$\ds{\ \bbox[#ffd,10px,border:1px solid navy]{{3 \over 8}\,\pi\ln\pars{2} - {\pi^{3} \over 64}}}$ $\ds{\ \approx\ 0.3321}$
\begin{align}
\int_{0}^{1}{\ln\pars{y}\over \pars{1 + y}^{\, 2}}\,\dd y & =
\lim_{\epsilon \to 0^{+}}\bracks{\ln\pars{\epsilon} +
\int_{\epsilon}^{1}{\dd y \over y\pars{1 + y}}}
\\[5mm] & =
\lim_{\epsilon \to 0^{+}}\bracks{\ln\pars{\epsilon} +
\int_{\epsilon}^{1}{\dd y \over y} - \int_{\epsilon}^{1}{\dd y \over 1 + y}} = -\ln\pars{2}
\end{align}