Tuberkulosis (TB) merupakan salah satu penyakit menular yang masih menjadi permasalahan kesehatan masyarakat di berbagai negara, termasuk Indonesia. Untuk memahami dinamika penyebarannya, penelitian ini membangun model matematika tipe VSEIT yang membagi populasi ke dalam lima kompartemen, yaitu divaksinasi (Vaccinated) (V), rentan (Susceptible) (S), terpapar (Exposed) (E), infektif (Infectious) (I), dan menjalani pengobatan (Treatment) (T). Model tersebut menggambarkan proses penularan TB melalui interaksi antara individu rentan dan infektif, perkembangan individu dari fase laten menjadi infektif, serta proses pengobatan yang berperan dalam menekan penyebaran penyakit. Analisis matematis dilakukan untuk menentukan titik ekuilibrium bebas penyakit dan titik ekuilibrium endemik beserta kestabilannya. Hasil analisis menunjukkan bahwa diperoleh bilangan reproduksi dasar sehingga titik ekuilibrium bebas penyakit bersifat stabil, sedangkan titik ekuilibrium endemik tidak memenuhi syarat eksistensi biologis. Selanjutnya, simulasi numerik pada titik ekuilibrium bebas penyakit dilakukan pada subpopulasi rentan (Susceptible) menggunakan Metode Euler, Runge–Kutta Orde Empat (RK4), dan Nonstandard Finite Difference (NSFD). Perbandingan ketiga metode dilakukan pada beberapa ukuran langkah waktu untuk mengevaluasi akurasi dan kestabilan numerik solusi. Hasil simulasi menunjukkan bahwa ketiga metode menghasilkan pola dinamika yang serupa dan konvergen menuju titik ekuilibrium bebas penyakit. Pada ukuran langkah yang kecil, solusi yang diperoleh hampir berimpit, sedangkan pada ukuran langkah yang lebih besar mulai terlihat perbedaan antar metode. Dengan demikian, model VSEIT mampu menggambarkan dinamika penyebaran tuberkulosis dengan baik, sedangkan metode Euler, RK4, dan NSFD memberikan solusi numerik yang konsisten dalam menyelesaikan model tersebut.