Background This study establishes a mathematically consistent and computational framework for the simultaneous identification of two time-dependent coefficients in a one-dimensional second-order parabolic partial differential equation. The considered problem is governed by nonlocal initial, boundary, and integral overdetermination conditions. Methods The direct problem is solved using the Crank-Nicolson finite difference method (FDM), which ensures unconditional stability and second-order accuracy in both spatial and temporal discretizations. The corresponding inverse problem is reformulated as a nonlinear regularized least-squares optimization problem and efficiently solved used the MATLAB subroutine
