A black hole described in SU(N) gauge theory consists of N D-branes. By separating one of the D-branes from others and studying the interaction between them, the black hole geometry can be probed. In order to obtain quantitative results, we employ the lattice Monte Carlo simulation. As a proof of the concept, we perform an explicit calculation in the matrix model dual to the black zero-brane in type IIA string theory. We demonstrate this method actually works in the high temperature region, where the stringy correction is large. We argue possible dual gravity interpretations.