Charge order in magnetite. An LDA+ U study