plot(1-(365!)/(365^x * (365-x)!)), 1<=x<=40