Rokok merupakan hasil olahan tembakau yang terbungkus, dihasilkan dari tanaman Nicotiana Tabacum, Nicotiana Rustica dan spesies lainnya atau sintetisnya yang mengandung nikotin dan tar dengan atau tanpa bahan tambahan. Banyak faktor yang menjadi latar belakang dan mempengaruhi perilaku merokok. Model matematika dapat digunakan untuk memahami penyebaran perokok, Model dibagi menjadi lima kelas yaitu: perokok potensial, perokok sesekali, perokok berat, perokok yang berhenti sementara dan perokok yang berhenti selamanya.Tujuannya menjelaskan terbentuknya model matematika perokok terhadap adanya pengaruh perokok berat, menentukan ekuilibrium model, menentukan kestabilan lokal model, dan menginterpretasikan model ke dalam simulasi numerik. Untuk pembentukan model digunakan asumsi-asumsi yang sesuai sehingga terbentuklah kompartemen dan model matematika perokok terhadap adanya pengaruh perokok berat, selanjutnya dengan menentukan perilaku model yaitu menentukan ekuilibrium model dan kestabilan lokal melalui linearisasi mengunakan matriks jacobian dan kriteria Routh-Hurwitz, serta menginterpretasikan model melalui hasil simulasi. Hasil yang diporeh adalah terbentuknya model matematika perokok terhadap adanya pengaruh perokok berat, diperoleh dua ekulibrium dari model yaitu ekulibrium bebas perokok stabil asimtotik lokal dan ekuilibrium endemik juga stabil asimtotik lokal.