plot(d(sum 1/(h-n) n=1 to (h-1))/dh)