An efficient approach for evaluating effective core potential integrals which involve projection operators has been implemented in the MAGIC quantum chemistry program. The methodology is presented and its performance is examined through illustrative calculations on transition metal and actinide compounds.