Solving the Bethe-Salpeter equation for positronium