integrate 2i*cos(cos(x))*sinh(sin(x))/x from x=0 to 5000 using the midpoint method with 20000 intervals