Bone tissue remodeling in the human body occurs continuously in various parts of the skeleton. The process is characterized by a stage of resorption of existing mineralized tissue by osteoclasts and a stage of subsequent formation of new bone tissue by osteoblasts. The intensity of bone remodeling at each local site is regulated by numerous intracellular and extracellular parameters of cell regulation. Mathematical modeling provides an important scientific approach for integrating existing knowledge about the process of interaction of bone cells, as well as for predicting the balance between bone formation and its subsequent resorption. This study presents a mathematical model that describes changes in the populations of osteoblasts and osteoclasts, as well as subsequent changes in bone mass in the local area of bone remodeling, taking into account the diffusion movement of cells. Numerical modeling was used to study the influence of parameters of internal and external regulation of osteoclasts and osteoblasts on the process of dynamic bone formation. It is shown that the used mathematical model can describe some pathological bone diseases.