Abstract
This article focuses on modeling bone formation process using a fractional differential approach, named bones remodeling process. The first goal of the work is to investigate existence and uniqueness of the proposed fractional differential model. The next goal is to investigate how similar is the proposed approach to the method based on system classical differential equations. The dynamical system of equations used is built upon three main parameters. These are chemical substances, namely, calcitonin secretion, osteoclastic and osteoblastic, which are involved in the bone’s formation process. We implement some numerical simulations to graphically show the impact of an arbitrary fractional order of derivative. We finally obtained that modeling bone formation process using fractional differential equations yielded comparable results with those obtained through a system of classical differential equations. Flexibility in the choice of the fractional order of derivative is an advantage as it helps in selecting the best fractional order of derivative.