plot (1-cos(x))/sin(x) for x=0.999*pi to 1.001*pi