An analytic energy gradient method for second-order quasidegenerate perturbation theory with multiconfigurational self-consistent field reference functions (MC-QDPT) is derived along the lines of the response function formalism (RFF). According to the RFF, the gradients are calculated without solving coupled perturbed equations. Instead, it is necessary to solve seven sets of linear equations in order to determine Lagrangian multipliers, corresponding to four sets of parameter constraining conditions and three sets of additional parameter defining conditions in the Lagrangian. Just one of these linear equations is a large scale linear equation; the others are reducible to just partial differentiations or simple equations solvable by straightforward subroutines.