A three-step method due to Nijhoff and Bobenko & Suris to derive a Lax pair for scalar partial difference equations (PΔEs) is reviewed. The method assumes that the PΔEs are defined on a quadrilateral, and consistent around the cube. Next, the method is extended to systems of PΔEs where one has to carefully account for equations defined on edges of the quadrilateral. Lax pairs are presented for scalar integrable PΔEs classified by Adler, Bobenko, and Suris and systems of PΔEs including the integrable two-component potential Korteweg–de Vries lattice system, as well as nonlinear Schrödinger and Boussinesq-type lattice systems. Previously unknown Lax pairs are presented for PΔEs recently derived by Hietarinta (J. Phys. A, Math. Theor. 44:165204, 2011). The method is algorithmic and is being implemented in Mathematica.