This work promotes the room temperature energy storage properties of the multiferroics. In this approach, impacts of PrFeO₃ doping on PT-based solid solutions (Pb1−xPrxTi1−xFexO₃, x = 0.21, 0.22, 0.23, 0.24, 0.25 and 0.26) have been explored. X-ray diffraction (XRD) patterns were used to estimate the crystallographic parameters, confirming the single phase tetragonal structure. The ferroelectric Curie temperature (TcFE) is observed to drop from 410 K to below room temperature as the Pr concentration increases. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The sample x = 0.24 achieved the maximum value of energy storage density of 362.25 mJ/cm³ with the efficiency of 40.5%. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The validity of magnetoelectric coupling in all samples was confirmed by magneto-dielectric studies and found that the sample x = 0.24 shows the maximum response with the coupling coefficient (γ) = 15.54 g²/emu².This work promotes the room temperature energy storage properties of the multiferroics. In this approach, impacts of PrFeO₃ doping on PT-based solid solutions (Pb1−xPrxTi1−xFexO₃, x = 0.21, 0.22, 0.23, 0.24, 0.25 and 0.26) have been explored. X-ray diffraction (XRD) patterns were used to estimate the crystallographic parameters, confirming the single phase tetragonal structure. The ferroelectric Curie temperature (TcFE) is observed to drop from 410 K to below room temperature as the Pr concentration increases. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The sample x = 0.24 achieved the maximum value of energy storage density of 362.25 mJ/cm³ with the efficiency of 40.5%. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The validity of magnetoelectric coupling in all samples was confirmed by magneto-dielectric studies and found that the sample x = 0.24 shows the maximum response with the coupling coefficient (γ) = 15.54 g²/emu².