plot 2x (x - log_2 x + 1) - gamma ((log_2 x) + 1) for x=670..680