plot (floor(x)+1/2) ln(x)-ln(floor(x)!),x-0.919023,x-0.91894,{x,1000.00009,1000.0002}