The AP1000 adopts the advanced Mechanical Shim(MSHIM) core control strategy to control the reactor power and axial power distribution. This paper presents the development of a dynamic simulation program of AP1000 reactor core in Matlab/Simulink. Firstly, an optimum nodal core model is developed, which both has satisfactory accuracy and fast simulation speed. Then, the MSHIM control system is established. Based on the simulation program, the typical load follow simulation is performed. The correctness of the program is preliminary validated by comparing the simulation results with the references.