logplot (365!/(365-x)!/365^x) from 1 to 365