Improving initial conditions for cosmological N-body simulations