Charge order in magnetite. An LDA+U study