In order to derive the orbital formulas for first-order corrections to the induced quantity ∆X, the second quantization approach was used, and the final equations were implemented in Python utilizing the approach based on the Psi4NumPy 63 framework -an interactive quantum-chemical environment designed for method development, rapid prototyping, and educational purposes. The basis for the newly developed expressions is the open-source quantum-chemical code Psi4. It is used to perform integral and SCF calculations and solve the CPHF and CPKS equations. All the necessary quantities are then exported from Psi4 into NumPy arrays, and their contractions are optimized and performed using the opteinsum module. The first-order corrections to the induced property have been implemented using molecular orbitals. The contractions of corresponding orbital expressions have a numerical cost that scales as N , where N is the number of molecular orbitals (due to the cost of the (Ξ B ) r a intermediate calculations). However, the cost of necessary CPHF (or CPKS) calculations scales as N , and therefore, it constitutes the actual performance limiting factor of the proposed theory. Additionally, the density-fitting approximation may be employed to reduce the scaling of the propSAPT theory to N 3 and the CPHF (or CPKS) calculations to N 5 .