We present a generalized nonlinear susceptibility retrieval method for metamaterials based on transfer matrices and valid in the nondepleted pump approximation. We construct a general formalism to describe the transfer matrix method for nonlinear media and apply it to the processes of three- and four-wave mixing. The accuracy of this approach is verified via finite element simulations. The method is then reversed to give a set of equations for retrieving the nonlinear susceptibility. Finally, we apply the proposed retrieval operation to a three-wave mixing transmission experiment performed on a varactor loaded split ring resonator metamaterial sample and find quantitative agreement with an analytical effective medium theory model.