Dynamics and bifurcations in predator-prey type of systems

Research team:

  1. Dr. Eric (ITB)
  2. Marwan, MSi (PhD student, ITB)
  3. Livia Owen, MSi (PhD student, ITB)


  1. Prof. H.W. Broer (Univ Groningen),
  2. Prof F. Verhulst (Univ Utrecht),
  3. Kie Van Ivanky Saputra PhD., (Universitas Pelita Harapan).
  4. Kus Prihantoso M.Si., (Universitas Negeri Jogjakarta)

This research was funded by:

  1. ITB, by Riset KK B ITB 2019.
  2. ITB, by Riset P3MI 2018.
  3. ITB, by Riset P3MI 2017.
  4. ITB, by Riset KK ITB 2017.
  5. ITB, by Riset KK ITB 2015.
  6. ITB, by Penguatan Riset Institusi 2011.
  7. KNAW (Royal the Netherlands Society for Science), by Spin PostDoctoral Program 2007.


Two dimensional predator-prey type of dynamical system.

We investigate a predator-prey model with a non-monotonic response function subjected to a small periodic forcing.  This forcing term is chosen to model for example, the small variations in the natural growth or carrying capacity due to seasonal change in the environtment.  This is investigated by using the so-called stroboscopic map.

We begin with studying the dynamics and bifurcations in the system, especially the bifurcation of periodic solution as we turn on the time periodic perturbation. There we have observed some interesting co-dimension two bifurcation for equilibria, namely: Bogdanov-Takens bifurcation, cusp bifurcation and Bautin bifurcation. As for bifurcation of periodic solution, we have observed the occurrence of Swallowtail-like bifurcation.



  1. Marwan, J. M. Tuwankotta, and E. Harjanto, Application of Lagrange multiplier method for computing fold bifurcation point in a two-prey one-predator dynamical system, Journal of the Indonesian Mathematical Society, 24 (2018), pp. 7–19.
  2. Eric Harjanto, J.M. Tuwankotta,Bifurcation of periodic solution in a Predator-Prey type of systems with non-monotonic response function and periodic perturbation , International Journal of Nonlinear Mechanics 85 (2016) 188-196.
  3. Eric Harjanto, J.M. Tuwankotta, Hilangnya Dua Bifurkasi Fold Tanpa Melalui Bifurkasi Cusp pada Sistem Predator-Prey dengan Faktor Pertahanan Grup dan Gangguan Berkala, Prosiding Konferensi Nasional Matematika XVII-2014, 11-14 Juni 2014, ITS, Surabaya.

To appear:

  1. L. Owen and J. Tuwankotta, Computation of cusp bifurcation point in a two-prey one predator model using lagrange multiplier method, in Proceedings of the International Conference on Applied Physics and Mathematics 2019, Chulalongkorn University, Bangkok, Thailand, 2019.
  2. Marwan, J.M. Tuwankotta, Infinitely many equilibria and some codimension one bifurcations in a subsystem of a two-prey one-predator dynamical system, Proceedings of International Conference on Mathematics and Natural Sciences, 2018, Institut Teknologi Bandung.

Leave a Reply

Your email address will not be published. Required fields are marked *


You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>

To use reCAPTCHA you must get an API key from https://www.google.com/recaptcha/admin/create