I met an interesting function defined by series of functions, which is given as follows \begin{equation*} f(n,s)=\frac{s}{n(n+1)}-\sum_{k=0}^{+\infty}\frac{2s}{(sk+n)(sk+n+1)(sk+n+2)}. \end{equation*} Here, $n=1,2,\cdots$, and $s>0$. When I studied the difference of digamma functions, the function $f(n,s)$ was derived through the series representation of the digamma function.
It is not difficult to show $f(n,1)=0$, no matter what $n$ is. Through numerical experiments in Matlab, I guess for any given $n$, the function is strictly increasing with respect to $s$. However, I failed to give a rigorous proof. Can someone give me a hand? Thanks a lot. The Matlab code and numerical result are given below.
function y = frctseries(s,n,m) % m:number of terms in the partial sum; y = 0; for k=0:m y = y + 2*s./((s*k+n).*(s*k+n+1).*(s*k+n+2)); end y = s./(n.*(n+1)) - y; end % z1 = frctseries(s,n(1),5000); z2 = frctseries(s,n(2),5000); z3 = frctseries(s,n(3),5000); figure plot(s,z1,'g','LineWidth',1); grid on; hold on; plot(s,z2,'r','LineWidth',1); hold on; plot(s,z3,'k','LineWidth',1); hold on; L = legend('$n=1$','$n=2$','$n=5$','FontSize',14); set(L,'Interpreter','latex'); 