plot digamma(1+x/2)-digamma(1/2+x/2) from 0 to 1000