<s>[INST] Evaluation of $\sum_{n=1}^{\infty}\frac{1}{n(2n+1)}=2-2\ln(2)$ I came across the following statements
$$\sum_{n=1}^{\infty} \frac{1}{n(2 n+1)}=2-2\ln 2 \qquad \tag{1}$$
$$\sum_{n=1}^{\infty} \frac{1}{n(3 n+1)}=3-\frac{3 \ln 3}{2}-\frac{\pi}{2 \sqrt{3}} \qquad \tag{2}$$
$$\sum_{n=1}^{\infty} \frac{1}{n(4 n+1)}=4-\frac{\pi}{2}-3 \ln 2 \qquad \tag{3}$$
$$\sum_{n=1}^{\infty} \frac{1}{n(6 n+1)}=6-\frac{\sqrt{3} \pi}{2}-\frac{3 \ln 3}{2}-2 \ln 2 \qquad \tag{4}$$
The (1) by partial fractions
$$\sum_{n=1}^{\infty}\frac{1}{n(2n+1)}=\sum_{n=1}^{\infty}\frac{1}{n}-\frac{2}{2n+1}$$
$$=\sum_{n=1}^{\infty}\frac{1}{n}-\frac{1}{n+\frac{1}{2}}$$
Recall the Digamma function
$$\psi(x+1)=\gamma+\sum_{n=1}^{\infty}\frac{1}{n}-\frac{1}{n+x}$$
Therefore
$$\sum_{n=1}^{\infty}\frac{1}{n}-\frac{1}{n+\frac{1}{2}}=\psi(1+\frac{1}{2})-\gamma$$
$$\sum_{n=1}^{\infty}\frac{1}{n(2n+1)}=\psi\left(\frac{3}{2}\right)-\gamma$$
In the same token we can derive the relation for the other three ralations. My Question is: can we calculate the values of the digamma function for those values without resorting in the Gauss´s Digamma formula?
$$\psi\left(\frac{r}{m}\right)=-\gamma-\ln (2 m)-\frac{\pi}{2} \cot \left(\frac{r \pi}{m}\right)+2 \sum_{n=1}^{\left\lfloor\frac{m-1}{2}\right\rfloor} \cos \left(\frac{2 \pi n r}{m}\right) \ln \sin \left(\frac{\pi n}{m}\right)$$

I tried this approach also, but I think the resulting integral is divergent $$\sum_{n=1}^{\infty} \frac{1}{n(2 n+1)}=\sum_{n=1}^{\infty} \frac{1}{n}\int_{0}^{1}x^{2n}dx=\int_{0}^{1}\sum_{n=1}^{\infty} \frac{x^{2n}}{n}=-\int_{0}^{1}\ln(1-x^2)dx $$
 [/INST] Here is another method to add to the list.  Consider the summation
$$\sum\limits_{n=1}^{\infty }{\frac{1}{n\left( 2n+1 \right)}}$$
Note that
$$\psi \left( -z \right)+\gamma \underset{z\to n}{\mathop{=}}\,\frac{1}{z-n}+{{H}_{n}}+\sum\limits_{k=1}^{\infty }{\left( {{\left( -1 \right)}^{k}}H_{n}^{k+1}-\zeta \left( k+1 \right) \right){{\left( z-n \right)}^{k}}}$$
for $n\ge 0$. So if we have a function $f$ that goes like $1/{{z}^{2}}$then taking  $\left( \psi \left( -z \right)+\gamma  \right)f\left( z \right)$over an infinitely large closed contour about $z=0$ just sums all the residues the total of which is zero.  We have then
$$\sum\limits_{n=1}^{\infty }{f\left( n \right)}=-\sum\limits_{f\left( z \right)}^{{}}{res\left( \psi \left( -z \right)+\gamma  \right)f\left( z \right)}$$
where the sum on the right is over all residues due to $f$.  In this case its very simple: there are two simple poles at $z=-1/2,0$.  The residue at $z=0$ is easy (use the expansion above).  We have
$$\sum\limits_{n=1}^{\infty }{\frac{1}{n\left( 2n+1 \right)}}=2+\psi \left( \tfrac{1}{2} \right)+\gamma $$
At this point I’m just going to assume its well known that $\psi \left( \tfrac{1}{2} \right)=-\gamma -2\log \left( 2 \right)$ 
</s>