Abstract: Direct production of ⁹⁹Mo and ^99mTc with accelerator has a good development prospect in the field of nuclear medicine. Accurate and reliable nuclear reaction data is of great significance for isotope production. In this paper, the cross sections for ^natMo(p, x)^96m+gTc, ¹⁰⁰Mo(p, x)⁹⁹Mo and ¹⁰⁰Mo(p, 2n)^99mTc reactions being taken from EXFOR database have been analyzed, and the experimental data were fitted by polynomial fit method. The excitation functions of ^natMo(p, x)^96m+gTc、¹⁰⁰Mo(p, x)⁹⁹Mo and ¹⁰⁰Mo(p, 2n)^99mTc reactions were calculated up to 40 MeV by using TALYS-1.95 and EMPIRE-3.2.3 code with different energy level density and pre-equilibrium emission models. For the ^natMo(p, x)^96m+gTc reaction, the calculation results predicted by preeqmode 1, 2 models are in good agreement with the experimental data when the input parameter ldmodel 6 of level density is used in TALYS-1.95 code. For the ¹⁰⁰Mo(p, x)⁹⁹Mo reaction, the calculation results predicted by preeqmode 4 models are in good agreement with the experimental data when the input parameter ldmodel 6 of level density is used in TALYS-1.95 code. For the ¹⁰⁰Mo(p, 2n)^99mTc reaction, the calculated results predicted by LEVDEN 2 in EMPIRE-3.2.3 code are in good agreement with the experimental data.