A numerical method is investigated to design gradient-index (GRIN) fiber probes. The GRIN fiber probe is composed of a single mode fiber (SMF), a no-core fiber (NCF), and a GRIN fiber lens. The optical software GLAD is adopted to simulate the optical performance of the probe. The simulation results show that, given the length of the GRIN fiber lens 0.1mm and the length of the NCF 0.36mm, the working distance is 0.73mm and the focus spot size 33μm, which are well agreement with the experimental data. As a result, the proposed numerical method is validated to be effective to design such GRIN fiber probes.