plot (floor(x)+1/2) ln(x)-ln(floor(x)!),x-1/2 ln(2pi),{x,1,20}