  All three-loop on-shell QCD Feynman integrals with two masses can be reduced
to 27 master integrals. Here we calculate these master integrals, expanded in
epsilon, both exactly in the mass ratio and as series in limiting cases.
