Two-loop maximal unitarity with external masses

We extend the maximal-unitarity method at two loops to double-box basis integrals with up to three external massive legs. We use consistency equations based on the requirement that integrals of total derivatives vanish. We obtain unique formulas for the coefficients of the master double-box integrals. These formulas can be used either analytically or numerically.