Starting from the orthogonal dynamics of any given set of variables with respect to the projection variable used to derive the Mori–Zwanzig equation, a set of coupled Volterra equations is obtained that relate the projected time correlation functions between all the variables of interest. This set of equations can be solved using standard numerical inversion methods for Volterra equations, leading to a very convenient yet efficient strategy to obtain any projected time correlation function or contribution to the memory kernel entering a generalized Langevin equation. Using this strategy, the memory kernel related to the diffusion of tagged particles in a bulk Lennard–Jones fluid is investigated up to the long-term regime to show that the repulsive–attractive cross-contribution to memory effects represents a small but non-zero contribution to the self-diffusion coefficient.