This is a special case of the solution by @Michal_Adamaszek at Writing x*log(1+x/y) .
rel_entr(x+1,x) + rel_entr(x,x+1)