Numerical computation of the transport matrix in a tokamak plasma with electrostatic turbulence