Hopf points of codimension two in a delay differential equation modeling leukemia

This paper continues the work contained in two previous papers, devoted to the study of the dynamical system generated by a delay differential equation that models leukemia. Here our aim is to identify degenerate Hopf bifurcation points. By using an approximation of the center manifold, we compute the first Lyapunov coefficient for Hopf bifurcation points. We find by direct computation, in some zones of the parameter space (of biological significance), points where the first Lyapunov coefficient equals zero. For these we compute the second Lyapunov coefficient, that determines the type of the degenerate Hopf bifurcation.