There is actually no need to express the result in terms of digamma function. A neat way to solve the problem is to choose points in the segments so that they form a geometric progression.
So we are going to find the following integral:
$$
I = \int_1^4 {1\over x^3}\mathop{dx}
$$
Consider the segment $[1, 4]$. Let's split it into $n$ parts with points $x_0, x_1,\dots x_n$ so that the points form a geometric progression. Let $q$ be the denominator of the geometric progression, then the segment in split by the points: $\{1, q, q^2, \dots, q^n\}$
and therefore:
$$
q^n = 4 \iff q =\sqrt[n]{4}
$$
By choosing those point in such a way we can see that the length of each interval is as follows:
$$
\begin{align}
\Delta x_1 &= q - 1\\
\Delta x_2 &= q^2 - q\\
&\dots \\
\Delta x_n &= q^n - q^{n-1}
\end{align}
$$
Now we have to choose points $\zeta_k$ in every subsegment. Let stick to the right side of each segment for simplicity. This gives:
$$
\zeta_k = \{q, q^2, q^3, \dots, q^n\}
$$
Let's calculate the value of the function in those points:
$$
f(\zeta_k) = \left\{{1\over q^3}, {1\over q^6}, \dots, {1\over q^{3k}} \right\}
$$
We are now ready to build the Riemann's sum:
$$
\begin{align}
S_n &= \sum_{k=1}^n f(\zeta_k)\Delta x_k\\
&= \sum_{k=1}^n {1\over q^{3k}}(q^k - q^{k-1})\\
&= {1\over q^3}(q-1) + {1\over q^6}(q^2 - q) + \cdots + {1\over q^{3n}}(q^n - q^{n-1})\\
&= {1\over q^3}(q-1) + {1\over q^5}(q - 1) + \cdots + {1\over q^{2n-1}}(q - 1)\\
&= (q - 1)\left( {1\over q^3} + {1\over q^5} + \cdots + {1\over q^{2n+1}} \right)\\
&= (q-1)\sum_{k=1}^n {1\over q^{2k+1}}
\end{align}
$$
But the sum in last expression is nothing but a geometric series and we know how to find its sum. So the expression above becomes:
$$
S_n = (q-1) \cdot \frac{{1\over q} - {1\over q^{2n+1}}}{q^2 - 1} = \frac{{1\over q} - {1\over q^{2n+1}}}{q+1}
$$
Remember our $q = \sqrt[n]{4}$:
$$
I = \lim_{n\to\infty} \frac{{1\over \sqrt[n]{4}} - {1\over (\sqrt[n]{4})^{2n+1}}}{\sqrt[n]{4}+1} = \boxed{{15\over 32}}
$$