A quasidegenerate perturbation theory (QDPT) with Epstein–Nesbet paritioning based on MCSCF reference functions is presented. The formulas for CSF-based matrix operations are derived. They are applied to the systems Be–H2 and CO, and compared with other methods. The potential energy curves agree well with those of the full-CI or MR-CI methods.