A modified school bullying model with recovery compartment: stability analysis and optimal control

Iranian Journal of Numerical Analysis and Optimization

مقاله 6، دوره 17، Issue 1 - شماره پیاپی 40، بهار 2027، صفحه 145-176   

نوع مقاله: Research Article

شناسه دیجیتال (DOI): 10.22067/ijnao.2026.99035.1883

نویسندگان

Ayomi Sasmito^{* 1}؛ Camelia K H¹؛ Rosa Maeni²

¹Business Mathematics, School of STEM, Prasetiya Mulya University, Tangerang, Indonesia.

²Business Mathematics, School of STEM, Prasetiya Mulya University, Tangerang, Indonesia

چکیده

This study proposes a compartmental dynamical model of school bullying transmission incorporating an explicit recovery mechanism through counselling and institutional intervention. The population is divided into susceptible, bully, exposed, violent, and recovered subgroups. The model’s qualitative properties are established, including positivity, boundedness, equilibrium analysis, and stability conditions derived using the next-generation matrix approach. The bullying-free equilibrium is shown to be both locally and globally asymptotically stable whenever $\mathcal{R}_0 < 1$, and the endemic equilibrium is shown to exist and be locally asymptotically stable when $\mathcal{R}_0 > 1$ under explicit Routh--Hurwitz conditions. Sensitivity analysis reveals that the perpetrator-generation rate $\alpha$ is the primary driver of bullying persistence, whereas recovery-related interventions substantially reduce long-term prevalence. An optimal control framework based on Pontryagin's Maximum Principle is then formulated using disciplinary and victim-protection interventions. Numerical simulations show that combining both interventions yields the largest reductions in bullying and in overall cost, indicating that integrated institutional responses outperform isolated policies.

کلیدواژه‌ها

School bullying؛ Compartmental model؛ Basic reproduction number؛ Sensitivity analysis؛ Optimal control