Identity with Pochhammer and harmonic numbers This came out of some work on the digamma function. 
Let $(x)_k=x(x+1)\cdots(x+k-1)$ denote the Pochhammer symbol. Then,

Question. Can you prove/disprove this identity?
  $$\pmb{\frac{(\frac12)_j^2}{j!^2}}\sum_{i=0}^{j-1}\frac4{2i+1}
=\sum_{i=0}^{j-1}\pmb{\frac{(\frac12)_i^2}{i!^2}}\frac1{j-i}.$$

I found this fascinating in view of fact that the factors in bold are able "go in and out" of the sum.
 A: Here's a sketch of a proof using "creative telescoping."
Let 
$$T(i,j) = \frac{j!^2}{(\tfrac12)_j^2}\cdot \frac{(\tfrac12)_i^2}{i!^2}\frac1{j-i}.$$
 Since the identity holds for $j=1$, it suffices to show that
$$\sum_{i=0}^{j} T(i,j+1) -\sum_{i=0}^{j-1}T(i,j)=\frac{4}{2j+1};$$
i.e., that
$$T(j,j+1) +\sum_{i=0}^{j-1}\bigl(T(i,j+1) - T(i,j)\bigr) = \frac{4}{2j+1}.\tag{1}$$
But it is easy to verify the indefinite summation
$$\sum_{i=0}^{k}\bigl(T(i,j+1) - T(i,j)\bigr)
  =\frac{j!^2}{(\tfrac32)_j^2}\cdot \frac{(\tfrac32)_k^2}{k!^2}\frac{1}{k-j}$$
from which $(1)$ follows by a straightforward computation.
A: The identity under question can be found in the paper S. Boettner, V.H. Moll The integrals in Gradshteyn and Ryzhik. Part 16: Complete elliptic integrals, pages 11-12 https://arxiv.org/abs/1005.2941

The proof given by OP is miraculously similar to the one in this paper, not only in the method used, but also in the choice of words and formatting of equations. By the way, one of the authors V. Moll, is a frequent collaborator of the OP's, and they have written numerous papers that start in the same way The integrals in Gradshteyn and Ryzhik. Part ..., e.g. https://arxiv.org/abs/1004.2440. Coincidences happen, and with some people they happen more often than with others, this is just law of probabilities. No foul play is suspected here. But it is really intriguing to know, how is the identity under question related to digamma function, as OP claims? Unfortunately his proof does not contain any digamma functions. But in the paper the identity naturally came out of some work on the elliptic integral. Please, @T.Amdeberhan, would you share your insight? It is really intriguing to know, how is digamma function related here?
Here is a screenshot of OP's own answer, just in case:

