func = lambda x: sc.digamma(x) - y