The present paper reports on the development of a relativistic CASPT2 program utilizing MO integrals formed by DIRAC software. This work is built upon previous research by one of the authors, which involved implementing a CASPT2 program utilizing MO integrals generated by the UTChem software. In addition to the previous work, the restricted active space secondorder perturbation theory (RASPT2) method was incorporated, and the program was parallelized.