Bullying menjadi masalah yang mengakar dari generasi sebelumnya ke generasi berikutnya. Kasus bullying terbanyak terjadi pada pelajar di tingkat SMP atau sederajat dengan rentang usia 13-15 tahun. Bullying dapat berdampak pada kesehatan mental pelajar. Emosi negatif akibat gangguan kesehatan mental dapat menyebar melalui interaksi sosial. Pada penelitian ini, dibahas analisis kestabilan model matematika kesehatan mental pelajar akibat bullying. Berdasarkan hasil analisis, didapatkan dua titik ekuilibrium, yaitu titik ekuilibrium bebas gangguan kesehatan mental pelajar akibat bullying (P_1^*) dan titik ekuilibrium endemik (P_2^*). Selanjutnya, menentukan bilangan reproduksi dasar dengan menggunakan metode generasi matriks selanjutnya. Berdasarkan hasil simulasi numerik dengan menggunakan software Matlab R2013a, diperoleh jika R_0<1 maka titik ekuilibrium bebas gangguan kesehatan mental pelajar akibat bullying bersifat stabil asimtotik lokal dan gangguan kesehatan mental akibat bullying tidak dapat menyebar dalam populasi pelajar. Jika R_0>1 maka titik ekuilibrium endemik bersifat stabil asimtotik lokal dan gangguan kesehatan mental akibat bullying dapat menyebar dalam populasi pelajar. Laju interaksi antara subpopulasi rentan mengalami gangguan kesehatan mental akibat bullying dengan subpopulasi yang mengalami gangguan kesehatan mental akibat bullying (α) memiliki pengaruh yang signifikan terhadap jumlah individu tiap subpopulasi.