integrate (1/Gamma(e^(ix))-1/Gamma(e^(-ix)))/x from 0 to 5000 using the midpoint method with 20000 intervals