Numerical solution of the inverse problem of reconstructing a distributed right-hand side of a parabolic equation