Plot[1 - e^(1000*f*ln(1000*f) + (1000-10)*ln(1000-10) - 1000*ln(1000) - (1000*f-10)*ln(1000*f-10)), {f, 0.5, 1}]