We consider a reformulation of QED in which covariant Green functions are used to solve for the electromagnetic field in terms of the fermion fields. The resulting modified Hamiltonian contains the photon propagator directly. A simple Fock-state variational trial function is used to derive relativistic two-fermion equations variationally from the expectation value of the Hamiltonian of the field theory. The interaction kernel of the equation is shown to be, in essence, the invariant M matrix in lowest order. Solutions of the two-body equations are presented for muonium like system for small coupling strengths. The results compare well with the observed muonium spectrum, as well as that for hydrogen and muonic hydrogen. Anomalous magnetic moment effects are discussed.