## RESEARCH ARTICLE • OPEN ACCESS

## Analisis Dinamik Model Penyebaran COVID-19 dengan Vaksinasi
Resmawan, Lailany Yahya, Revandi S. Pakaya, Hasan S. Panigoro, and Agusyarif Rezka Nuha

## Volume 3, Issue 1, Pages 29–38, June 2022
Received 10 Januari 2022, Accepted 21 June 2022, Published Online 1 July 2022 To Cite this Article : Resmawan dkk., “Analisis Dinamik Model Penyebaran COVID-
19 dengan Vaksinasi”, Jambura J. Biomath , vol. 3, no. 1, pp. 29–38, 2022,
https://doi.org/10.34312/jjbm.v3i1.13176
© 2022 by author(s)

## JOURNAL INFO • JAMBURA JOURNAL OF BIOMATHEMATICS
S Homepage : http://ejurnal.ung.ac.id/index.php/JJBM/index ǚ Journal Abbreviation : Jambura J. Biomath. Ѓ Frequency : Biannual (June and December) Ŕ Publication Language : English (preferable), Indonesia  DOI : https://doi.org/10.34312/jjbm ʀ Online ISSN : 2723-0317 Ǜ Editor-in-Chief : Hasan S. Panigoro Ŗ Publisher : Department of Mathematics, Universitas Negeri Gorontalo ź Country : Indonesia ̎ OAI Address : http://ejurnal.ung.ac.id/index.php/jjbm/oai  Google Scholar ID : XzYgeKQAAAAJ ͐ Email : editorial.jjbm@ung.ac.id

## JAMBURA JOURNAL • FIND OUR OTHER JOURNALS
Jambura Journal of
Mathematics Jambura Journal of Mathematics Education Jambura Journal of Probability and Statistics EULER : Jurnal Ilmiah Matematika, Sains, dan Teknologi

## Analisis Dinamik Model Penyebaran COVID-19 dengan Vaksinasi
Resmawan 1 , ∗ , Lailany Yahya 2 , Revandi S. Pakaya 3 , Hasan S. Panigoro 4  , dan Agusyarif Rezka Nuha 5
1 , 2 , 3 , 4 , 5 Jurusan Matematika, Universitas Negeri Gorontalo, Bone Bolango 96119, Indonesia
ABSTRAK. Coronavirus Disease 2019 (COVID-19) adalah virus jenis baru dari keluarga besar virus yang ditu- larkan antara manusia dan hewan (ditularkan secara zoonosis) yang pertama kali ditemukan di Kota Wuhan, Provinsi Hubei, China pada akhir 2019 yang hingga saat ini masih menyebar luas dan menjadi ancaman dis- eluruh dunia termasuk Indonesia. Pada artikel ini dibahas tentang model matematika penyebaran COVID-19 dengan vaksinasi. Dalam hal ini, populasi manusia dibagi menjadi 5 kelas, yaitu kelas suspected, vaccine, ex- posed, infected dan recovered . Model yang dikontruksi membentuk model SVEIR yang mempunyai dua titik kesetimbangan, yaitu titik kesetimbangan bebas penyakit dan endemik. Analisis kestabilan menunjukkan bahwa titik ekuilibrium bebas penyakit stabil asimtotik lokal dan global jika R 0 < 1 dan tidak stabil jika R 0 > 1 . Kemudian dilakukan analisis sensitivitas untuk mengetahui parameter yang sangat berpengaruh terhadap model serta selanjutnya diberikan simulasi secara numerik untuk mengambarkan perilaku dari model yang telah diperoleh berdasarkan analisis sensitivitas bilangan reproduksi dasar, diperoleh beberapa parameter yang berpengaruh terhadap penyebaran COVID-19. Hasil simulasi numerik menunjukkan bahwa vaksinasi dapat menekan penambahan populasi terinfeksi dan bergantung pada tingkat efektivitas vaksinasi.
ABSTRACT. Coronavirus Disease 2019 (COVID-19) is a new type of virus from a large family of viruses transmitted between humans and animals (zoonotically transmitted) that was first discovered in Wuhan City, Hubei Province, China in late 2019 which is still widespread and threat throughout the world including Indonesia. This article discussed about the mathematical model of the spread of COVID-19 with vaccinations. In this case, the human population is divided into 5 classes, namely the suspected, vaccine, exposed, infected and recovered classes. The constructed model forms an SVEIR model that has two equilibrium points, namely disease-free and endemic equilibrium points. Stability analysis shows that the equilibrium point is stable local and global asymptotic if R 0 < 1 and unstable if R 0 > 1 . Then a sensitivity analysis was carried out to determine the parameters that greatly affect the model as well as furthermore, numerical simulations are given to describe the behavior of the model that has been obtained based on the analysis of the sensitivity of basic reproductive numbers, obtained several parameters that affect the spread of COVID-19. Numerical simulation results show that vaccination can suppress the addition of infected populations and depend on the level of effectiveness of vaccination.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution-NonComercial 4.0 International License. Editorial of JJBM: Department of Mathematics, Uni- versitas Negeri Gorontalo, Jln. Prof. Dr. Ing. B. J. Habibie, Bone Bolango 96119, Indonesia.

## 1. Pendahuluan
Coronavirus merupakan keluarga besar virus yang ditu- larkan antara manusia dan hewan (ditularkan secara zoonosis). Terdapat dua jenis coronavirus yang diketahui menyebabkan penyakit pada manusia, yaitu Middle East Respiratory Syndrome (MERS-CoV) dan Severe Acute Respiratory Syndrome (SARS-CoV) . Pada akhir tahun 2019 telah ditemukan jenis baru coronavirus yang da- pat menimbulkan gejala gangguan pernapasan akut seperti de- mam, batuk, dan sesak napas di Kota Wuhan, Provinsi Hubei, China yang mempunyai nama resmi Coronavirus Disease 2019 (COVID-19) [1] .
World Health Organization (WHO) atau Badan Kesehatan Dunia secara resmi menetapkan virus corona (COVID-19) seba- gai suatu pandemi pada tanggal 11 Maret 2020 [2] . Virus ini ke- mudian telah menyebar secara luas dan mendadak menjadi anca-
∗ Corresponding Author.
Check for updates
Research Article

## ARTICLE HISTORY
Received 10 Januari 2022
Accepted 21 June 2022
Published 1 July 2022

## KEYWORDS
Dynamic Analysis Lyapunov function Sensitivity Analysis COVID-19 Vaccinate
man di dunia, terutama setelah merenggut ribuan nyawa manu- sia dalam waktu yang relatif singkat [3] . Penyakit ini menjadi salah satu topik pembicaraan utama karena dampak yang diak- ibatkannya, sehingga membuat banyak pihak untuk terjun lang- sung mengambil langkah-langkah aktif dalam menemukan solusi penanganan kesehatan dan sarana pencegahan agar dampak yang ditimbulkan tidak akan lebih meluas. Dari sudut pandang mate- matis, fokus ini berkaitan erat dengan implementasi model mate- matika untuk mempermudah memahami dan menganalisa suatu permasalahan dalam kehidupan sehari-hari.
Model matematika adalah metode mendeskripsikan masalah yang kompleks dalam bentuk matematika. Selain itu, model matematika dapat diterapkan pada berbagai bidang seperti biologi, ekologi, epidemiologi dan kedokteran [4] . Salah satu permasalahan yang dapat diselesaikan dengan menggu- nakan model matematika yakni masalah penyebaran COVID-19.
Email : resmawan@ung.ac.id (R. Resmawan), lailany.math@gmail.com (L. Yahya), hspanigoro@ung.ac.id (H.S. Panigoro), and agusyarif@ung.ac.id (A.R.
Nuha) Homepage : http://ejurnal.ung.ac.id/index.php/JJBM/index / E-ISSN : 2723-0317 © 2022 by the Author(s).
Pemodelan matematika memudahkan seseorang untuk dapat memahami dan menemukan solusi yang tepat untuk penanganan COVID-19.
Beberapa penelitian telah membahas model matematika COVID-19. Seperti Khan dan Atangana [5] yang merancang model matematika dengan mengeksplorasi dinamika infeksi COVID-19 menggunakan operator turunan orde-fraksional ABC. Resmawan dkk [6] membahas model transisi COVID-19 dengan melibatkan intervensi karantika. Annas dkk [7] menggunakan model SEIR un- tuk mempelajari penyebaran COVID-19 di Indonesia. Resmawan dan Yahya [8] melakukan analisis sensitivitas untuk mengetahui parameter mana yang paling dominan untuk mempengaruhi en- demisitas COVID-19. Gebremeskel dkk. [9] mempelajari wabah pandemi COVID-19 dan menyelidiki lebih lanjut tren masa de- pannya di Ethiopia menggunakan data COVID yang dilaporkan harian. Pengembangan model matematika tentang penularan COVID-19 masih sangat terbuka dan terus dilakukan hingga saat ini. Terlebih, setelah ditemukannya vaksin untuk mencegah COVID-19.
Pada penelitian ini, diperkenalkan model baru yang men- gacu pada model [6] dengan melakukan modifikasi yaitu de- ngan menggabungkan kelas populasi manusia yang terinfeksi tanpa gejala dan populasi manusia yang dikarantina kedalam ke- las populasi manusia yang terpapar, hal ini dimaksudkan sesuai dengan keadaan dunia nyata bahwa tidak diberlakukannya lagi karantina. Modifikasi lainnya yaitu berupa penambahan kelas populasi vaksin yang mewakili individu yang telah divaksinasi. Vaksin ini dimaksudkan untuk memberi perlindungan tubuh agar tidak jatuh sakit akibat COVID-19 dengan cara menimbulkan atau menstimulasi kekebalan spesifik dalam tubuh. Model yang telah diperoleh akan dianalisis dengan memperhatikan titik kesetim- bangan bebas penyakit, selanjutnya menyelidiki stabilitas lokal dan global keseimbangan bebas penyakit. Pada bagian akhir akan diberikan simulasi numerik untuk mengonfirmasi hasil analisa di- namika penyebaran COVID-19 dengan vaksinasi.

## 2. Metode
Penelitian ini dilakukan menggunakan metode studi pus- taka dengan langkah-langkah sebagai berikut:
1. Kajian pustaka terkait pemodelan matematika penyebaran COVID-19
2. Memodifikasi model matematika penyebaran COVID-19 de- ngan vaksinasi dan melakukan penskalaan model
3. Menentukan titik kesetimbangan dan analisis kestabilan titik kesetimbangan modifikasi model matematika penye- baran COVID-19 dengan vaksinasi
4. Melakukan simulasi numerik terhadap modifikasi model ma- tematika penyebaran COVID-19 dengan vaksinasi
5. Menyelidiki efektivitas vaksinasi terhadap model yang di- modifikasi

## 3. Hasil dan Pembahasan

## 3.1. Model Matematika
Model pada penelitian ini membahas model matematika penyebaran COVID-19 tipe SVEIR. Model ini terbagi dalam enam kelas populasi, yaitu manusia yang rentan (S), manusia yang di- vaksin (V), manusia yang terpapar termasuk manusia terinfeksi tanpa gejala dalam masa inkubasi (E), manusia terinfeksi dis-
ertai gejala klinis (I) dan manusia yang pulih dari COVID-19 (R). Model yang dimodifikasi dengan pertimbangan atau batasan dalam merumuskan model matematika menggunakan asumsi se- bagai berikut:
1. Laju rekrutmen kelahiran hanya ditujukan pada populasi (S) dimana setiap individu yang baru lahir akan dimasukkan kedalam kelas manusia rentan
2. Tingkat kematian terdiri atas tingkat kematian alami yang ditujukan kepada seluruh kelas populasi dan tingkat kema- tian disebabkan oleh COVID-19 yang hanya ditujukan pada populasi (I)
3. Kelompok individu dalam kelas vaksin akan terpapar apabila tingkat kekebalan dari vaksin menurun
4. Kelompok individu vaksin merupakan individu yang telah melewati dua kali tahapan vaksinasi
5. Penularan COVID-19 terjadi ketika adanya kontak antara in- dividu pada populasi (S) dan populasi (V) dengan individu pada populasi (E) atau individu pada populasi (I)
6. Individu pada populasi (E) sudah mencakup individu yang terinfeksi tanpa gejala akan melewati masa inkubasi dan ma- suk kedalam populasi (I)
7. Individu pada populasi (I) akan mengalami pemulihan dan berpindah ke populasi (R)
8. Individu yang telah dinyatakan pulih dari COVID-19 tidak akan kembali rentan
Berdasarkan asumsi-asumsi diatas, diagram kompartmen model matematika penyebaran COVID-19 dengan vaksinasi diberikan oleh Gambar 1 sehingga diperoleh model dengan sistem per- samaan diferensial sebagai berikut.
dS dt = Π N − (1 − φ ) ηS ( ξ 1 E + ξ 2 I ) N − ( µ + φα ) S dV dt = φαS − ( µ + γ ( ξ 3 E + ξ 4 I ) N ) V dE dt = (1 − φ ) ηS ( ξ 1 E + ξ 2 I ) N + γV ( ξ 3 E + ξ 4 I ) N − ( ω + µ ) E dI dt = ωE − ( ρ + µ + δ ) I dR dt = ρI − µR (1)
dengan S + V + E + I + R = N . Untuk menyederhanakan model (1) dilakukan nondimensionalisasi sebagai berikut.
s = S N , v = V N , e = E N , i = I N , r = R N sehingga diperoleh
ds dt = Π − (1 − φ )( ξ 1 e + ξ 2 i ) ηs − ( µ + φα ) s dv dt = φαs − ( µ + γ ( ξ 3 e + ξ 4 i )) v de dt = (1 − φ )( ξ 1 e + ξ 2 i ) ηs + γv ( ξ 3 e + ξ 4 i ) − ( ω + µ ) e di dt = ωe − ( ρ + µ + δ ) i dr dt = ρi − µr (2)

## 3.2. Penentuan Titik Kesetimbangan
Untuk menentukan titik kesetimbangan pada sistem (2) da- pat ditentukan dengan menyelesaikan persamaan tersebut de- ngan ds dt = dv dt = de dt = di dt = dr dt = 0 . Dari penyelesaian persamaan (2) diperoleh dua titik kesetimbangan sebagai berikut.
1. Titik kesetimbangan bebas penyakit COVID-19 ( E 0 )
E 0 = ( s 0 , v 0 , 0 , 0 , 0) = ( Π µ + αφ , αφ Π µ 2 + αφµ , 0 , 0 , 0 ) (3)
2. Titik kesetimbangan endemik
E 1 = ( s ∗ , v ∗ , e ∗ , i ∗ , r ∗ )
dengan
s ∗ = ω Π (1 − φ )(( ρ + µ + δ ) ξ 1 + ωξ 2 ) ηi ∗ + ( αφ + µ ) ω , v ∗ = φαω 2 Π (( ξ 4 ω + ξ 3 ( ρ + µ + δ )) γi ∗ + µω ) ((( ρ + µ + δ ) ξ 1 + ξ 2 ω )(1 − φ ) ηi ∗ + ω ( αφ + µ )) , e ∗ = ( ρ + µ + δ ) i ∗ ω , i ∗ = (+ ω + µ − γξ 3 v ∗ − (1 − φ ) ηξ 1 s ∗ ) e ∗ (1 − φ ) ηξ 2 s ∗ + γξ 4 v ∗ , r ∗ = ρi ∗ µ .
3.3. Bilangan Reproduksi Dasar ( R 0 )
Bilangan reproduksi dasar yang dinotasikan R 0 adalah nilai harapan banyaknya infeksi tiap satu satuan waktu. Infeksi ini ter- jadi pada suatu populasi rentan yang dihasilkan oleh satu individu terinfeksi [10] . Dalam menentukan bilangan reproduksi dasar di- gunakan metode Next Generation Matrix (NGM) [11] . Next Gen- eration Matrix (K) didefinisikan sebagai berikut.
K = F V − 1
dengan
F = [ (1 − φ ) ηξ 1 s 0 + γξ 3 v 0 (1 − φ ) ηξ 2 s 0 + γξ 4 v 0 ω 0 ]
dan
V = [ ω + µ 0 0 ρ + µ + δ ]
sehingga diperoleh matriks K sebagai berikut
K =    (1 − φ ) ηξ 1 s 0 + γξ 3 v 0 ω + µ (1 − φ ) ηξ 2 s 0 + γξ 4 v 0 ρ + µ + δ ω ω + µ 0    .
R 0 diperoleh dari nilai eigen positif terbesar dari matriks K [8] yaitu
R 0 = R A + R S (4)
dengan
R A = (1 − φ ) ηξ 1 s 0 + γξ 3 v 0 2( ω + µ ) , R S = √ R A 2 + (1 − φ ) ηωξ 2 s 0 + γωξ 4 v 0 ( ω + µ )( ρ + µ + δ ) .

## 4. Kestabilan Titik Kesetimbangan
Kestabilan lokal dan global titik kesetimbangan E 0 diberikan oleh Teorema 1 dan 2 .
Teorema 1. Titik kesetimbangan bebas Penyakit COVID-19 E 0 stabil asimtotik lokal jika R 0 < 1 dan tidak stabil jika R 0 > 1 .
Bukti . Matriks Jacobi yang di evaluasi pada titik kesetimbangan E 0 diberikan oleh
J ( E 0 ) =       A 11 0 A 13 A 14 0 A 21 A 22 A 23 A 24 0 0 0 A 33 A 34 0 0 0 A 43 A 44 0 0 0 0 A 54 A 55       (5)
dengan
A 11 = − ( µ + αφ ) , A 13 = − (1 − φ ) ηξ 1 s 0 ,
A 14 = − (1 − φ ) ηξ 2 s 0 , A 21 = αφ, A 22 = − µ,
A 23 = − γξ 3 v 0 , A 24 = − γξ 4 v 0 , A 33 = (1 − φ ) ηξ 1 s 0 + γξ 3 v 0 − ( ω + µ ) = (2 R A − 1)( ω + µ ) A 34 = (1 − φ ) ηξ 2 s 0 + γξ 4 v 0 , A 43 = ω,
A 44 = − ( ρ + µ + δ ) , A 54 = ρ, A 55 = µ.
Dari matriks Jacobi J ( E 0 ) diperoleh persamaan karakteristik
( λ − A 11 )( λ − A 22 )( λ − A 55 ) ( λ 2 − ( A 33 + A 44 ) λ + A 33 A 44 − A 34 A 43 ) = 0 ,
sehingga diperoleh nilai eigen
λ 1 = A 11 = − µ − αφ, λ 2 = A 22 = − µ, λ 3 = A 55 = − µ, dan λ 4 , 5 adalah akar-akar dari
λ 2 − ( A 33 + A 44 ) λ + A 33 A 44 − A 34 A 43 = 0 , (6)
Jelas bahwa λ i < 0 , ∀ i = 1 , 2 , 3 dan Re ( λ 45 ) < 0 jika A 33 + A 44 < 0 dan A 33 A 44 − A 34 A 43 > 0 (Kriteria Routh-Hurwitz [4] ). Perhatikan bahwa
R S = √ R A 2 + (1 − φ ) ηωξ 2 s 0 + γωξ 4 v 0 ( ω + µ )( ρ + µ + δ ) > √ R A 2 = R A
sehingga jika R 0 < 1 maka diperoleh
1 > R A + R S > 2 R A 0 > 2 R A − 1
0 > (2 R A − 1)( ω + µ ) > (2 R A − 1)( ω + µ ) − ( δ + µ + ρ ) = A 33 + A 44
S V E I R Π N µS αφS µV ωE ρI µR µE µI µS γV ( ξ 3 E + ξ 4 I ) N R ηS (1 − φ )( ξ 1 E + ξ 2 I ) N
Gambar 1. Diagram kompartmen model matematika penyebaran COVID-19 dengan vaksinasi
Tabel 1. Deskripsi biologis variabel dan parameter pada model (1)
Parameter Deskripsi Π Tingkat perekrutan kelahiran µ Kematian alami manusia η Peluang infeksi setelah kontak antara individu terinfeksi dengan individu rentan tanpa vaksinasi ξ 1 Tingkat kontak antara individu rentan (S) dan individu terpapar (E)
ξ 2 Tingkat kontak antara individu rentan (S) dan individu terinfeksi dengan gejala klinis (I)
ξ 3 Tingkat kontak antara individu yang telah divaksinasi (V) dan individu terpapar (E) ξ 4 Tingkat kontak antara individu yang telah divaksinasi (V) dan individu terinfeksi dengan gejala klinis (I)
ω Tingkat tranmisi setelah menyelesaikan masa inkubasi dan berpindah ke (I)
ρ Tingkat pemulihan individu terinfeksi dengan gejala klinis berpindah ke kelas individu yang pulih φ Proporsi individu rentan yang diberikan vaksin δ Kematian karena terinfeksi COVID-19 γ
Peluang setelah kontak antara individu terinfeksi dengan individu yang telah divaksinasi
α Laju vaksinasi
dan
A 33 A 44 − A 34 A 43 = − (2 R A − 1)( ω + µ )( δ + µ + ρ ) − ((1 − φ ) ηωξ 2 s 0 + γωξ 4 v 0 ) = − (2 R A − 1)( ω + µ )( δ + µ + ρ ) − ( R S 2 − R A 2 ) ( ω + µ )( ρ + µ + δ ) = (2 R A − 1 − R 0 )) ( R 0 − 1)( ω + µ )( ρ + µ + δ ) = ( ( R 0 − 1) 2 − 2 R S ( R 0 − 1) ) ( ω + µ )( ρ + µ + δ ) > 2 R S (1 − R 0 )( ω + µ )( ρ + µ + δ ) > 0
Dengan demikian terbukti bahwa E 0 stabil asimtotik lokal jika R 0 < 1 . Dengan cara yang serupa dapat ditunjukkan sebaliknya untuk R 0 > 1 juga memiliki kestabilan yang berkebalikan.
Teorema 2. Jika 2 ω + µ < 2( ω + µ ) R S < ( µ + δ ) ω ( ρ + µ + δ ) R 0 dan
R 0 < 1 maka titik kesetimbangan bebas Penyakit COVID-19 E 0 stabil asimtotik global di Ω := { ( s, v, e, i, r ) | s 0 s + v v 0 < 1 }
Bukti . Diberikan fungsi Lyapunov definit positif yang merupakan kombinasi antara fungsi Volterra dan linear sebagai berikut.
V ( s, v, e, i, r ) = ( s − s 0 − s 0 ln s s 0 ) + ( v − v 0 − v 0 ln v v 0 ) + e + i + r.
Dengan demikian diperoleh turunan fungsi Lyapunov sebagai berikut:
dV dt = ( s − s 0 s ) ds dt + ( v − v 0 v ) dv dt + de dt + di dt + dr dt = ( s − s 0 s ) (Π − (1 − φ ) ηs ( ξ 1 e + ξ 2 i ) − ( µ + φα ) s ) + ( v − v 0 v ) ( φαs − ( µ + γ ( ξ 3 e + ξ 4 i )) v ) + (1 − φ ) ηs ( ξ 1 e + ξ 2 i ) + γv ( ξ 3 e + ξ 4 i )
− ( ω + µ ) e + ωe − ( ρ + µ + δ ) i + ρi − µr = Π ( 2 − s 0 s − s s 0 ) + αφs ( 2 − v 0 v − v v 0 ) − µ ( v + e + i + r ) − δi + γξ 3 v 0 e + γξ 4 v 0 i + αφs 0 + (1 − φ ) ηξ 1 s 0 e + (1 − φ ) ηξ 2 s 0 i + αφsv v 0 − αφs = Π ( 2 − s 0 s − s s 0 ) + αφs ( 2 − v 0 v − v v 0 ) + ( ω + µ )( ρ + µ + δ )( R S − R A ) R 0 i ω
+ 2( ω + µ ) R A e − µ ( v + e + i + r ) − δi + αφs ( s 0 s + v v 0 − 1 ) .
Karena R 0 < 1 dan R A = R 0 − R S maka
dV dt ≤ Π ( 2 − s 0 s − s s 0 ) + αφs ( 2 − v 0 v − v v 0 ) + (2( ω + µ )(1 − R S ) − µ ) e − µ ( v + r ) + 2( ω + µ )( ρ + µ + δ ) R S R 0 − ( µ + δ ) ω ω i + αφs ( s 0 s + v v 0 − 1 ) .
Karena 2 ω + µ < 2( ω + µ ) R S < ( µ + δ ) ω ( ρ + µ + δ ) R 0 maka di Ω diperoleh
dV dt ≤ Π ( 2 − s 0 s − s s 0 ) + αφs ( 2 − v 0 v − v v 0 )
Berdasarkan rata-rata aritmatik dan geometrik diperoleh 2 − s 0 s − s s 0 ≤ 0 dan 2 − v 0 v − v v 0 ≤ 0 sehingga dV dt < 0 . Akibatnya, titik kesetimbangan E 0 stabil asimtotik global di Ω .

## 5. Analisis Sensitivitas dan Simulasi Model

## 5.1. Nilai Parameter
Untuk melakukan simulasi maka diperlukan nilai nilai pa- rameter yang mendukung. Nilai parameter diambil berdasarkan data sekunder yang diperoleh dari berbagai referensi serta pemi- lihan parameter lainnya didasarkan pada asumsi. Nilai parameter yang digunakan dapat dilihat pada Tabel 2 dan nilai awal dari se- tiap populasi dapat dilihat pada Tabel 3 .

## 5.2. Analisis Sensitivitas
Analisis sensitivitas yang dilakukan pada penelitian ini adalah analisis sensitivitas lokal untuk mengetahui parameter- parameter yang berpengaruh terhadap angka reproduksi dasar. Angka reproduksi dasar yang dianalisis yaitu pada saat kondisi R 0 > 1 yang menunjukkan bahwa terjadinya penyebaran COVID- 19. Berdasarkan hal tersebut, maka perlu dilakukan analisis sen- sitivitas terhadap parameter-parameter untuk mengetahui upaya yang harus dilakukan agar menekan penyebaran COVID-19.
Analisis sensitivitas dilakukan dengan menentukan indeks sensitivitas dari parameter-parameter. Berdasarkan [14] indeks
Tabel 2. Estimasi nilai parameter model
Parameter R 0 < 1 Sumber R 0 > 1 Sumber Π 4 . 1 × 10 − 5 [12] 4 . 1 × 10 − 5 [12] µ 3 . 57 × 10 − 5 [6] 3 . 57 × 10 − 5 [6] η 0 . 2 [6] 0 . 2 [6] ξ 1 0 . 09 [6] 0 . 09 [6] ξ 2 0 . 07 [6] 0 . 07 [6] ξ 3 0 . 03 Asumsi 0 . 03 Asumsi ξ 4 0 . 02 Asumsi 0 . 02 Asumsi ω 0 . 005 [6] 0 . 005 [6] ρ 0 . 27 [13] 0 . 0027 Asumsi δ 0 . 001 [12] 0 . 001 [12] α 0 . 001 Asumsi 0 . 0003 Asumsi φ 0 . 2 Asumsi 0 . 2 Asumsi γ 0 . 01 [13] 0 . 01 [13] Tabel 3. Nilai awal model Populasi Nilai s (0) 267564000 270000000 v (0) 20534823 270000000 e (0) 2500000 270000000 i (0) 3372374 270000000 r (0) 2730720 270000000
sensitivitas normalisasi dari bilangan reproduksi dasar didefinisi- kan sebagai berikut.
C R 0 p = ∂R 0 ∂p × p R 0
Dimana p adalah parameter yang terdapat pada angka re- produksi dasar yaitu Π , η, ξ 1 , ξ 2 , ξ 3 , ξ 4 , µ, α, γ, ω, ρ, δ, φ Analisis sensitivitas untuk persamaan C R 0 p dilakukan de- ngan menginput nilai parameter pada Tabel 2 Selanjutnya diper- oleh R 0 = 2 , 5574 dan nilai indeks sensitivitas masing-masing parameter yang ditunjukkan pada Tabel 4 berikut.
Indeks sensitivitas dalam Tabel 4 diurutkan berdasarkan parameter dengan nilai indeks sensitivitas terbesar dan dapat disimpulkan bahwa tiap parameter memberikan pengaruh yang bervariasi terhadap R 0 . Indeks sensitivitas yang bernilai positif menunjukkan adanya pengaruh perubahan nilai parameter yang berbanding lurus dengan perubahan angka reproduksi dasarnya, artinya jika nilai parameter diperbesar (diperkecil) sementara pa- rameter lainnya dibuat konstan maka akan berkontribusi pada pe- ningkatan (penurunan) angka reproduksi dasar. Sebagai contoh, parameter η mempunyai nilai indeks sensitivitas positif. Hal ini menunjukkan bahwa jika η ditingkatkan maka nilai R 0 juga akan meningkat.
Akan tetapi, indeks sensitivitas yang bernilai negatif menunjukkan adanya pengaruh perubahan nilai parameter yang berbanding terbalik dengan perubahan angka reproduksi
Tabel 4. Nilai Indeks Sensitivitas Parameter
Parameter Indeks Sensitivitas Π 1 η 0 , 965498 α − 0 , 614466 φ − 0 , 614466 ω − 0 , 488833 ξ 2 0 , 487968 ξ 1 0 , 47753 µ − 0 , 397352 ρ − 0 , 357624 δ − 0 , 141725 γ 0 , 0345017 ξ 3 0 , 0183923 ξ 4 0 , 0161094
dasarnya, artinya jika nilai parameter diperbesar (diperkecil) se- mentara parameter lainnya dibuat konstan, maka akan berkon- tribusi pada penurunan (peningkatan) angka reproduksi dasar. Sebagai contoh, parameter α mempunyai nilai indeks sensitivi- tas negatif yang artinya bahwa jika nilai α ditingkatkan maka R 0 akan mengalami penurunan.
Analisis sensitivitas menunjukkan bahwa parameter pelu- ang infeksi setelah kontak antara individu terinfeksi dengan ind- vidu rentan tanpa vaksinasi dan parameter laju vaksinasi adalah parameter yang paling sensitif terhadap perubahan angka reprod- ksi dasar. indeks sensitivitas C R 0 η = 0 , 965498 menunjukkan bahwa dengan meningkatkan peluang infeksi setelah kontak an- tara individu terinfeksi dengan indvidu rentan tanpa vaksinasi sebesar 10 % akan meningkatkan nilai R 0 sebesar 9 , 65498 % dan C R 0 α = − 0 , 614466 menunjukkan bahwa dengan meningkatkan laju pemberian vaksinasi sebesar 10 % akan menurunkan nilai R 0 sebesar 6 , 14466 %

## 5.3. Simulasi Numerik
(i) Dinamika Populasi Untuk R − 0 < 1 . Berdasarkan per- samaan (2) diperoleh satu titik bebas penyakit COVID-19. Titik tetap bebas penyakit COVID-19 diperoleh dengan men- subtitusikan parameter parameter pada Tabel 2 dengan bi- langan reproduksi dasar R 0 = 0 , 563329 . Titik kesetimba- ngan bebas penyakit diberikan oleh.
E 0 = (0 . 17394 , 0 . 97450 , 0 , 0 , 0)
Selanjutnya untuk menganalisis titik kesetimbangan bebas penyakit COVID-19, dapat dilakukan dengan cara mensubs- titusikan nilai parameter-parameter pada Tabel 2 ke per- samaan (3) sehingga diperoleh nilai eigen untuk titik kese- timbangan bebas penyakit COVID-19 sebagai berikut.
λ E 0 = − 0 . 271146 , − 0 . 002198 ,
− 0 . 000235 , − 0 . 000035 , − 0 . 00035
Karena semua nilai eigennya bernilai negatif maka dapat disimpulkan bahwa titik kesetimbangan bebas penyakit COVID-19 stabil. Dinamika populasi penyebaran penyakit COVID-19 dengan vaksinasi untuk R 0 < 1 dapat dilihat pada Gambar 2 . Gambar 2 menunjukkan bahwa setiap populasi stabil menuju titik kesetimbangan bebas penyakit COVID-
19. Dimana populasi manusia yang rentan terhadap COVID- 19 mengalami penurunan jumlah populasi dari nilai awal hingga mencapai kondisi stabil dititik S = 0 . 174 atau se- kitar 46 . 980 . 000 orang. Sedangkan populasi manusia yang telah diberi vaksin mengalami peningkatan jumlah popu- lasi dari nilai awal hingga mencapai kondisi stabil dititik V = 0 . 975 atau sekitar 263 . 250 . 000 orang. Adapun po- pulasi manusia yang terinfeksi COVID-19 mengalami penu- runan jumlah populasi dari nilai awal kemudian mengalami fluktuasi dan mencapai kondisi stabil dititik I = 0 . Sedan- gkan populasi manusia yang terpapar COVID-19, dan po- pulasi manusia yang telah sembuh dari COVID-19 masing- masing mengalami peningkatan jumlah populasi dari nilai awal dan kemudian mengalami penurunan jumlah populasi hingga mencapai kondisi stabil E = R = 0 .
(ii) Dinamika Populasi R 0 > 1 Berdasarkan persamaan (2) diperoleh titik penyakit COVID- 19. Titik kesetimbangan endemik penyakit COVID-19 diper- oleh dengan cara mensubstitusikan parameter-parameter pada Tabel 2 pada persamaan (2) dengan nilai bilangan re- produksi dasar R 0 = 2 . 55746 titik kesetimbangan endemik penyakit COVID-19 diberikan oleh.
E 1 = (0 . 167911 , 0 . 261206 , 0 . 00509969 ,
0 . 00670007 , 0 . 506728)
Selanjutnya dilakukan pelinearan dan perhitungan pada sis- tem persamaan (2) disekitar titik tetap endemik. Sehingga diperoleh nilai eigen sebagai berikut.
λ E 1 = ( − 0 . 00635747 + 000001 i, − 0 . 00011526 + 0 . 00064901 i, − 0 . 00011526 − 0 . 00064901 i, − 0 . 0000398727 + 0 . 000001 i, − 0 . 0000357 + 0 . 000001 i )
Berdasarkan Gambar 3 dapat dilihat bahwa populasi menuju titik kesetimbangan endemik atau stabil disekitar titik tetap penyakit. Populasi manusia yang rentan terhadap COVID-19 mengalami penurunan jumlah populasi dari nilai awal dan mengalami fluktuasi hingga mencapai kondisi stabil pada saat S = 0 . 168 atau sekitar 45 . 360 . 000 orang. Popu- lasi manusia yang telah diberi vaksin mengalami pening- katan jumlah populasi dari nilai awal hingga stabil pada titik V = 0 . 261 atau sekitar 70 . 470 . 000 orang. Populasi manusia yang telah sembuh dari COVID-19 mengalami pe- ningkatan jumlah populasi manusia dan mengalami fluktuasi hingga stabil pada titik R = 0 . 506 atau sekitar 136 . 620 . 000 orang. Populasi manusia yang terpapar COVID-19 dan popu- lasi manusia yang terinfeksi COVID-19 masing-masing men- galami peningkatan jumlah populasi dari nilai awal dan ke- mudian mengalami penurunan jumlah populasi hingga men- capai stabil dititik E = 0 . 005 , I = 0 . 006 .

## 5.4. Simulasi dengan Variasi Nilai Parameter
a. Simulasi laju individu yang divaksinasi terhadap penyebaran penyakit COVID-19 pada dinamika populasi.
Simulasi ini dilakukan untuk mengetahui pengaruh α yaitu laju vaksinasi terhadap penyebaran penyakit COVID-19. Dalam hal ini
Gambar 2. Dinamika Populasi R 0 < 1
Gambar 3. Dinamika Populasi R 0 > 1
akan ditunjukkan bahwa peningkatan atau penurunan nilai para- meter α dapat mengubah nilai bilangan reproduksi dasar ( R 0 ) yang didefinisikan pada persamaan (4) . Adapun nilai-nilai para- meter α yang dapat mempengaruhi nilai R 0 dapat dilihat pada Tabel 4 .
Berdasarkan Tabel 4 dapat diperhatikan bahwa dengan meningkatnya nilai α dapat menyebabkan turunnya nilai R 0 . Se- hingga upaya yang harus dilakukan untuk menekan laju penye-
baran penyakit COVID-19 yaitu dengan meningkatkan laju vaksi- nasi terhadap populasi rentan. Selanjutnya hasil simulasi nu- merik disajikan pada Gambar 4 .
Perubahan dari nilai laju individu yang divaksinasi mempe- ngaruhi jumlah individu pada populasi manusia rentan, populasi manusia yang telah divaksin, jumlah populasi manusia terpapar dan terinfeksi serta manusia yang telah sembuh dari COVID-19. Jika laju dari individu yang telah divaksin diperbesar dan nilai

## Gambar 4. Simulasi Laju Individu yang divaksinasi
Tabel 5. Hasil Simulasi Variasi Nilai Parameter α terhadap R 0 Parameter α R 0 Π 0.0003 2.5574 µ 0.00033 2.4141 η 0.00033 2.4141 γ 0.00039 2.1729
parameter lain tidak berubah, maka hal ini akan menyebabkan semakin berkurangnya jumlah manusia pada populasi rentan, manusia yang terpapar dan manusia yang terinfeksi serta manu- sia yang sembuh dari COVID-19. Sedangkan jumlah manusia pada populasi manusia yang telah divaksinasi meningkat. Hal ini me- nunjukkan bahwa jumlah manusia rentan, manusia yang terpa- par dan manusia yang terinfeksi serta manusia yang sembuh dari COVID-19 dapat berkurang dengan adanya upaya vaksinasi yang diberikan pada individu rentan.
Gambar 4 menunjukkan bahwa pemberian vaksinasi pada individu yang rentan terhadap COVID-19 mempunyai peran cukup penting dalam menekan penyebaran COVID-19 yang di-
tandai dengan berkurangnya jumlah populasi rentan, manusia terpapar dan manusia yang terinfeksi COVID-19 dari jumlah po- pulasi sebelumnya.
a. Simulasi peluang infeksi setelah kontak antara individu terinfeksi dengan individu rentan tanpa vaksinasi.
Simulasi ini dilakukan untuk mengetahui pengaruh η yaitu peluang infeksi setelah kontak antara individu terinfeksi dengan individu rentan tanpa vaksinasi. Berdasarkan hasil analisis sen- sitivitas, η merupakan parameter yang berpengaruh terhadap angka reproduksi dasar R 0 . Oleh karena itu, akan ditunjukkan bahwa perubahan nilai parameter η dapat mengubah nilai R 0 . Adapun nilai-nilai parameter η yang dapat mempengaruhi nilai R 0 dapat dilihat pada Tabel 5 .
Berdasarkan Tabel 5 dapat diperhatikan bahwa dengan meningkatkan nilai η dapat menyebabkan naiknya R 0 . hal ini sesuai dengan hasil analisis sensitivitas bahwa parameter η berbanding lurus dengan R 0 . Selain itu, dari Tabel 5 bahwa jika nilai η ditingkatkan akan menyebabkan makin menyebarnya COVID-19 sehingga populasi orang yang rentan akan berkurang. Selanjutnya, hasil simulasi numerik disajikan pada Gambar 5 .
Perubahan dari nilai peluang infeksi setelah kontak antara
Gambar 5. Simulasi Peluang Infeksi Setelah Kontak Antara Individu Terinfeksi dengan Individu Rentan Tanpa Vaksinasi
Tabel 6. Hasil Simulasi Variasi Nilai Parameter η terhadap R 0 Parameter α R 0 Π 0.21 2.6813 µ 0.23 2.9290 η 0.25 3.1767 γ 0.27 3.4244
individu terinfeksi dengan individu rentan tanpa vaksinasi mem- pengaruhi jumlah individu pada populasi manusia yang rentan, populasi manusia yang diberikan vaksin dan populasi manusia yang sembuh. Jika peluang infeksi setelah kontak antara individu terinfeksi dengan individu rentan tanpa vaksinasi diperbesar dan parameter lain dibuat konstan, maka hal ini akan menyebabkan semakin berkurangnya jumlah populasi manusia yang rentan dan manusia yang diberi vaksinasi. Sedangkan populasi manusia yang telah sembuh dari COVID-19 mengalami peningkatan jumlah po- pulasi manusia.
Gambar 5 menunjukkan bahwa jika peluang infeksi setelah
kontak antara individu terinfeksi dengan individu rentan tanpa vaksinasi semakin besar maka akan berdampak pada menyebarlu- asnya COVID-19 dikalangan masyarakat. Sebaliknya jika peluang infeksi setelah kontak antara individu terinfeksi dengan individu rentan tanpa vaksinasi dikurangi ini akan membantu mencegah terjadinya penyebaran COVID-19.

## 6. Kesimpulan
Model yang diperoleh memiliki dua titik kesetimbangan yaitu titik kesetimbangan tanpa penyakit dan titik kesetimbangan endemik. Titik kesetimbangan tanpa penyakit bersifat stabil asimtotik lokal dan global jika R 0 < 1 . Sedangkan titik kesetim- bangan endemik bersifat asimtotik lokal dan global jika R 0 > 1 . Analisis sensitivitas menunjukkan bahwa parameter yang paling sensitif terhadap perubahan angka reproduksi dasar adalah laju vaksinasi. Hasil simulasi numerik menunjukkan bahwa pada kon- disi R 0 < 1 jumlah manusia yang terinfeksi COVID-19 berku- rang dan pada waktu tertentu akan menghilang. Sedangkan pada kondisi R 0 > 1 jumlah manusia yang terinfeksi COVID-19 akan bertambah dan dapat terjadi penyebaran COVID-19. Peningka-
tan laju vaksinasi menyebabkan bilangan reproduksi dasar berku- rang. Hal tersebut menunjukkan bahwa dengan adanya pening- katan laju vaksinasi terhadap individu rentan dapat mengurangi jumlah penyebaran COVID-19 sehingga tidak terjadi penyebaran COVID-19.

## References
[1] Kementerian Kesehatan, “Pertanyaan dan Jawaban Terkait COVID-19,” 2020. [Online]. Available: https://infeksiemerging.kemkes.go.id/ uncategorized/qna-pertanyaan-dan-jawaban-terkait-covid-19 (Accessed 20 November 2021).
[2] I. Darti, A. Suryanto, H. S. Panigoro, and H. Susanto, “Forecasting COVID- 19 Epidemic in Spain and Italy Using A Generalized Richards Model with Quantified Uncertainty,” Communication in Biomathematical Sciences , vol. 3, no. 2, pp. 90–100, 2020. DOI: 10.5614/cbms.2020.3.2.1
[3] D. Aldila, M. Z. Ndii, and B. M. Samiadji, “Optimal control on COVID-19 erad- ication program in Indonesia under the effect of community awareness,” Mathematical Biosciences and Engineering , vol. 17, no. 6, pp. 6355–6389, 2020. DOI: 10.3934/mbe.202033
[4] M. Z. Ndii, Pemodelan Matematika Dinamika Populasi Dan Penyebaran Penyakit Teori, Aplikasi, Dan Numerik . Deepublish, 2018. ISBN 9786232090354.
[5] M. A. Khan and A. Atangana, “Modeling the dynamics of novel coron- avirus (2019-nCov) with fractional derivative,” Alexandria Engineering Journal , vol. 59, no. 4, pp. 2379–2389, 2020. DOI: 10.1016/j.aej.2020.02.033
[6] R. Resmawan, A. R. Nuha, and L. Yahya, “Analisis Dinamik Model Transmisi COVID-19 dengan Melibatkan Intervensi Karantina,” Jambura Journal of Math-
ematics , vol. 3, no. 1, pp. 66–79, 2021. DOI: 10.34312/jjom.v3i1.8699
[7] S. Annas, M. Isbar Pratama, M. Rifandi, W. Sanusi, and S. Side, “Stability analysis and numerical simulation of SEIR model for pandemic COVID-19 spread in Indonesia,” Chaos, Solitons and Fractals , vol. 139, p. 110072, 2020. DOI: 10.1016/j.chaos.2020.110072
[8] R. Resmawan and L. Yahya, “Sensitifity Analysis of Mathematical Model of Coronavirus Disease (COVID-19) Transmission,” CAUCHY , vol. 6, no. 2, p. 91, 2020. DOI: 10.18860/ca.v6i2.9165
[9] A. A. Gebremeskel, H. W. Berhe, and H. A. Atsbaha, “Mathematical mod- elling and analysis of COVID-19 epidemic and predicting its future sit- uation in Ethiopia,” Results in Physics , vol. 22, p. 103853, 2021. DOI: 10.1016/j.rinp.2021.103853
[10] R. Resmawan, “Model Epidemik SEIRS-SEI Penyebaran Penyakit Malaria de- ngan Vaksinasi dan Pengobatan,” in Prosiding Seminar Nasional Matematika IndoMS Wilayah Sulawesi . IndoMS, 2017, pp. 128–40.
[11] J. Giesecke, Modern Infectious Disease Epidemiology , third edit ed. London: CRC Press, 2017. ISBN 9781315222714. DOI: 10.1201/9781315222714
[12] B. Rinaldi, “Pemodelan Matematika Penyebaran Penyakit Corona Viruses Dis- ease 2019 (COVID-19) pada Kasus Penggunaan Masker Kesehatan,” Skripsi, Universitas Jambi, 2021.
[13] M. Azizah, “Model Matematika Penyebaran Penyakit Coronavirus Disease 2019 (COVID-19) Dengan Vaksinasi, Isolasi mandiri, dan Karantina di Rumah Sakit,” Ph.D. dissertation, UIN Syarif Hidayatullah, 2021.
[14] N. Chitnis, J. M. Hyman, and J. M. Cushing, “Determining Important Parame- ters in the Spread of Malaria Through the Sensitivity Analysis of a Mathemat- ical Model,” Bulletin of Mathematical Biology , vol. 70, no. 5, pp. 1272–1296, 2008. DOI: 10.1007/s11538-008-9299-0