Bifurcation Analysis About a Mathematical Model of Somitogenesis Based on the Runge–Kutta Method

Abstract

In this paper, we investigated the modified two dimensional model which can explain somite patterning in embryos. It is suitable for exploring a design space of somitogenesis and can explain aspects of somitogenesis that previous models cannot. Here, we mainly studied the non-diffusing case. We have used the Hopf bifurcation theorem, the Center manifold theorem and Runge---Kutta method in our investigation. First, we investigate its dynamical behaviors and put forward a sufficient condition for the oscillation of the small network. Then, we give the mathematical simulation based on the Runge---Kutta method. In the process of solving ordinary differential equations, the four order Runge---Kutta method has the advantages of high accuracy, convergence and stability (under certain conditions), which can change the step size and do not need to calculate higher order derivatives. Therefore, it has become the most commonly used numerical solution. At the same time, we get the sufficient condition in which the bistable state of the system exists and give the numerical simulation. Because somitogenesis occupies an important position in the process of biological development, and as a pattern process can be used to study pattern formation and many aspects of embryogenesis. So our study have a great help for embryonic development, gene expression, cell differentiation. In addition, it is beneficial to study the clone animal variation problem of spinal bone number and is of great help to the treatment and prevention of defects of human spine disease.