Merokok merupakan salah satu ancaman kesehatan terbesar masyarakat yang saat ini dihadapi dunia dan dapat menyebabkan berbagai penyakit serius. Salah satu langkah untuk meningkatkan kesadaran berupa introspeksi diri di kalangan masyarakat dengan memberikan informasi tentang bahaya merokok melalui media. Penyebaran perokok dapat dimodelkan menggunakan pemodelan matematika yang dibentuk dalam suatu sistem persamaan diferensial biasa. Tujuan penelitian ini adalah menjelaskan terbentuknya model matematika penyebaran perokok dengan adanya media informasi dan introspeksi diri, menganalisis kestabilan model, menggunakan kontrol optimal penyebaran perokok dan membuat simulasi numerik. Hasil dari penelitian ini adalah terbentuknya model matematika  penyebaran perokok dengan adanya media informasi dan introspeksi diri dengan diperoleh dua titik ekuilibrium yaitu titik ekuilibrium bebas perokok dan titik ekuilibrium penyebaran perokok. Kemudian, diperoleh bilangan reproduksi dasar (R_0) dengan menggunakan metode Next Generation Matrix dan analisis kestabilan model menggunakan metode linearisasi dan Kriteria Routh-Hurwitz. Hasil analisis kestabilan model di titik ekuilibrium bebas perokok stabil asimtotik lokal dengan syarat (R_0)<1 sedangkan di titik ekuilibrium penyebaran perokok  stabil asimtotik lokal dengan syarat (R_0)>1. Selanjutnya, diperoleh kontrol optimal pada model matematika perokok dengan dua kontrol yang diberikan yaitu konseling dan terapi, di mana dapat meminimalkan subpopulasi Smokers dan memaksimalkan subpopulasi Smokers that quit smoking. Simulasi numerik diberikan untuk mendukung analisis kestabilan di titik ekuilibrium dan menunjukkan ketika dua kontrol diberikan, tingkat perokok menurun seiring waktu.

Kata kunci: Perokok, Media Informasi, Introspeksi Diri, Next Generation Matrix, Analisis Kestabilan, Kontrol Optimal.