Model Pertumbuhan Hidup Nyamuk Aedes Aegypti

Jurnal Konvergensi Vol. 4, No. 2, Oktober, 2014

Model Pertumbuhan Hidup Nyamuk Aedes Aegypti Wenny Wulandari Program Studi Matematika FMIPA UAD

Abstrak Penelitian ini dilakukan untuk mengetahui pemodelan matematika mengenai pertumbuhan siklus nyamuk aedes aegypti. Model matematika mengenai pertumbuhan nyamuk aedes aegypty pada penulisan ini terdapat 3 tahap (telur, larva dan nyamuk dewasa). Hasil dari analisis diperoleh dua titik equilibrium yaitu pada keadaan tidak ada pertumbuhan dan keadaaan yang terdapat pertumbuhan. Hasil simulasi yang diperoleh dengan menggunakan metode euler yaitu pada saat parameter batasannya kurang dari 1 ( < 1) diperoleh keadaan yang menunjukkan tidak ada pertumbuhan sehingga menyebabkan tidak ada keadaan endemik dan pada saat parameter batasannya lebih dari 1 ( > 1) diperoleh keadaan yang menunjukkan ada pertumbuhan sehingga menyebabkan keadaan endemik.

1. Pendahuluan Salah satu penyakit yang masih mewabah di Indonesia setiap tahunnya adalah penyakit Demam Berdarah Dangue (DBD). Sampai sekarang DBD masih merupakan masalah kesehatan masyarakat dan masih banyaknya laporan mengenai meninggalnya penderita karena kurang cepat ditangani (Gandahusada, dkk. 2000). Penyakit ini disebabkan oleh nyamuk Aedes Aegypti yang menularkan virusnya ke manusia. Cara penularannya adalah setiap kali nyamuk menghisap darah maka nyamuk segera mengekskresikan air liurnya yang mengandung zat anti pebekuan darah supaya darah bisa dihisap termasuk darah yang mengandung virus dangue, infeksi virus dangue dalam tubuh nyamuk dapat mengakibatkan perubahan prilaku pada peningkatan kompetensi nyamuk yaitu kemampuan nyamuk menyebarkan virus. Penyakit DBD ini tidak punya pengobatan ataupun vaksin,sehingga salah satu cara untuk mencegah atau mengurangi proses transmisi virus dangue dengan bergantung pada pengontrolan utama nyamuk (Moulay, 2012). Pengontrolan nyamuk Aedes aegypti dapat dilakukan dengan mengetahui dinamika pertumbuhan nyamuk tersebut. Karena dengan mengetahui dinamika nyamuk tersebut kita bisa menggambarkan atau memprediksi pertumbuhan atau penurunan populasi dari nyamuk tersebut sehingga berguna untuk pengontrolan nyamuk tersebut. Aplikasi matematika yang dianggap berkembang sangat pesat adalah aplikasi matematika dalam bidang biologi. Permasalahan biologi dapat ditransformasikan dalam model matematika dengan menggunakan beberapa asumsi. Salah satu aplikasi matematika dalam bidang biologi adalah model pertumbuhan yaitu model matematika. Sehingga pertumbuhan nyamuk dapat dimodelkan kedalam bentuk matematika.

Model Pertumbuhan Hidup Nyamuk Aedes Aegypti Wenny Wulandari

105


2. Perumusan Model Dalam penelitian ini tahapan larva dan pupa tidak akan dibedakan karena tahapan ini merupakan tahapan yang sama. Sehingga ada 3 tahapan pada pertumbuhan siklus nyamuk ini yaitu laju pertumbuhan telur yang dinotasikan dengan ( ) kemudian larva yang dinotasikan dengan ( ) dan nyamuk dewasa dinotasikan dengan ( ). a. Telur Pada tahapan telur laju pertumbuhannya dipengaruhi oleh jumlah nyamuk dewasa yang bertelur. Kemudian berkurang karena perubahan telur menjadi larva Pertumbuhan telur juga berkurang karena kematian alami dari telur itu sendiri. b. Larva Pada tahapan ini laju pertumbuhannya dipengaruhi oleh pertumbuhaan telur yang menjadi larva kemudian berkurang karena perubahan larva menjadi nyamuk aedes aegypti dewasa. Laju pertumbuhan larva juga berkurang karena kematian dari larva itu sendiri. c. Nyamuk Pada tahapan nyamuk laju pertumbuhannya dipengaruhi oleh pertumbuhan larva yang menjadi nyamuk dewasa kemudian berkurang karena kematian alami dari nyamuk. Nyamuk yang bertelur tidak akan menyebabkan berkurangnya laju pertumbuhan nyamuk tersebut. Dibawah ini merupakan diagram yang menggambarkan siklus hidup nyamuk aedes

106



aegypti Dari deskripsi di atas dapat diformulasikan siklus hidupnya ke dalam model matematika: 1) Laju pertumbuhan telur Parameter b menunjukkan rata-rata banyaknya telur dari nyamuk aedes aegypti. Jumlah telur akan berkurang karena adanya pertumbuhan telur yang menjadi larva yang dinyatakan dengan r. Selain Pertumbuhan telur yang menjadi larva jumlah telur akan berkurang karena kematian alami dari telur itu sendiri yang dinyatakan dalam d. Sehingga model matematika yang dapat dibentuk adalah sebagai berikut: =

−

−

2) Laju pertumbuhan larva Pada pertumbuhan larva, jumlah populasi larva akan bertambah dari telur yang menjadi larva yang dinyatakan dalam r. Kemudian larva akan menjadi nyamuk yang menyebabkan jumlah populasi larva berkurang dinyatakan dengan. Selain itu populasi larva juga berkurang karena kematian alami larva yang dinyatakan dengan parameter. Model laju pertumbuhan larva adalah sebagai berikut: =

−

−

3) Laju pertumbuhan nyamuk Jumlah populasi nyamuk bertambah dari pertumbuhan larva yang berkembang menjadi nyamuk yaitu dan populasinya berkurang karena kematian nyamuk yang dinyatakan dengan . Sehingga laju pertumbuhannya dinyatakan sebagai berikut: =

−

Namun model pertumbuhan pada siklus nyamuk ini didasarkan pada asumsi bahwa populasi bertambah dengan laju yang sebanding dengan besarnya populasi. Asumsi ini masuk akal untuk kondisi ideal yaitu lingkungan tak terbatas, nutrisi yang mencukupi tidak adanya pemangsa dan lain-lain. Dalam kenyataanya suatu lingkungan memiliki sumber daya yang terbatas. Dalam penelitian ini diasumsikan bahwa populasi dibatasi oleh ketersediaan makanan yang ada di air. Sehingga pada tahapan larva yang berada di air akan dibatasi. Sehingga model pertumbuhannya menjadi 1− Jadi model populasi dibentuk dari siklus nyamuk sebagai berikut: =


−

−

107


=

1−

=

−

−

−

Keterangan: E : Telur L : Larva/ Pupa A : Nyamuk dewasa : rata-rata banyaknya telur yang dihasilkan dari nyamuk betina : rata-rata banyaknya telur yang menjadi larva :rata-rata banyaknya larva yang menjadi nyamuk dewasa : tingkat kematian telur : tingkat kematian larva : tingkat kematian nyamuk dewasa : ketersediaan makanan (jumlah protozoa)

3. Analisis Model Salah satu masalah yang penting dalam dalam pemodelan adalah analisis dinamika sistem persamaan diferensial nonlinear. Pertama tentukan titik equilibrium model. Titik equilibrium memenuhi persamaan-persamaan 0= − − 0= 0=

1−

−

−

−

Diperoleh dua titik equilibrium,yaitu titik yang menunjukkan tidak adanya pertumbuhan yaitu = ( , , ) = (0,0,0)dan titik adanya pertumbuhan =

( ,

,

)

(

)

,

1−

,

Dalam model ini digunakan parameter batasan yang dibentuk dari perkembangan siklus nyamuk yaitu perkembangan dari telur, larva dan kemudian nyamuk. Perkembangan telur diperoleh dari rata-rata banyaknya telur yang dihasilkan dikalikan dengan usia hidup dari telur itu, perkembangan larva diperoleh dari rata-rata banyak nya larva dikalikan dengan usia hidup dari larva. Sedangkan perkembangan nyamuk diperoleh dari rata-rata banyaknya nyamuk dikalikan dengan usia hidup nyamuk. Sehingga parameter batasan yang dibentuk dari masa hidup nyamuk ini adalah sebagai berikut =

108

.



Selanjutnya dicari kestabilan dari titik keseimbangan.Kestabilan dari suatu titik equilibrium dapat dilihat dari nilai eigennya. Nilai eigen sendiri dapat dicari dari persamaan karakteristik yang merupakan determinan dari matriks Jacobi. Matriks Jacobi pada titik equilibrium yang menunjukkan tidak adanya pertumbuhan

+( + ) − 0 Maka didapat nilai eigen dari

0 +( + −

− 0 +

)

berupa persamaan karakteristiknya ( )= + + +

Dimana = ( + )+( + )+ = ( + )( + ) + ( + ) + ( + ) = ( + )( + )(1 − ) Dengan menggunakan kriteris Routh Hurwitz diperoleh nilai determinan yang positif Matriks Jacobi pada titik equilibrium yang menunjukkan adanya pertumbuhan −( + )

0

=

−( +

)

1−

1

+1

0 Maka didapat nilai eigen dari

0 −

berupa persamaan karakteristiknya ( )=

+

+

+

Dimana =

+

= ( + )+( + = = ( + )( +

)

1−

+

)

1−

+ 1

+1 +

+

+1 +( + )

= = ( + )( +

1

+ ( +

) 1−

1

− )

1−

1

+1

−

Dengan menggunakan kritetris Routh Hurwitz diperoleh nilai determinan yang positif.

4. Simulasi Numerik Berikut diberikan parameter untuk simulasi numerik dari model ini


109


1. Simulasi pada saat titik equilibriumnya = ( , , ) = ( , , ) yaitu tidak adanya pertumbuhan Pada simulasi ini menggunakan nilai awal (0) = 25, (0) = 10, (0) = 12

Gambar 4.4.1 : Grafik laju pertumbuhan Telur, Larva, Nyamuk dengan = . , = . , = . , = . , = . , = , = . , = . 110



Dari parameter tersebut diperoleh titik kesetimbangan telur, larva dan nyamuk berturut turut adalah Dari gambar di atas dapat dilihat: Laju pertumbuhan telur mengalami penurunan hal ini disebabkan karena laju kematian yang besar kemudian penurunan populasi telur disebabkan oleh telur yang berubah jadi larva. Populasi ini akan terus berkurang hingga menuju titik kesetimbangannya dan kemudian stabil pada → ∞. Pada populasi tersebut sudah tidak terjadi penambahan populasi telur. Laju pertumbuhan Larva mengalami penurunan karena tingkat kematian larva lebih besar dan perubahan telur yang menjadi larva lebih sedikit dari laju perubahan larva menjadi nyamuk dan terus menurun hingga menuju titik kesetimbangannya dan kemudian stabil pada → ∞. Laju pertumbuhan Nyamuk mengalami kenaikan dikarenakan laju perubahan larva menjadi nyamuk sangat besar kemudian mengalami penurunan dikarenakan kematian kemudian terus menurun menuju titik kestabilannya.

2. Simulasi pada saat titik equilibriumnya =(

,

,

)=

(

)

,

−

,

yaitu adanya

pertumbuhan Pada simulasi ini menggunakan nilai awal yaitu (0) = 50, (0) = 40, (0) = 5

Gambar 4.4.2 : Grafik Laju pertumbuhan Telur, Larva, Nyamuk dengan = . , = . , = . , = . , = . , = . , = . , = , = . Model Pertumbuhan Hidup Nyamuk Aedes Aegypti Wenny Wulandari

111


Dari parameter tersebut diperoleh titik kesetimbangan telur, larva dan nyamuk berturut turut adalah (128.6420, 41.6800, 231.5556) Dari gambar di atas dapat dilihat : Laju pertumbuhan telur mengalami penurunan pada awalnya karena laju kematian telur yang besar tapi kemudian mengalami kenaikkan karena populasi nyamuk yang meningkat. Populasi ini akan terus naik hingga menuju titik kesetimbangannya dan kemudian stabil pada → ∞. Pada populasi ini sudah tidak terjadi penambahan populasi telur. Laju pertumbuhan Larva mengalami penurunan pada awalnya dikarenakan tingkat kematian larva yang besar tapi kemudian mengalami kenaikkan karena jumlah telur yang berubah menajdi larva meningkat.Populasi iniakan terus meningkat sehingga mencapai titik kesetimbangannya dan populasi akan mengalami kestabilan. Laju pertumbuhan Nyamuk mengalami kenaikan dikarenakan laju perubahan larva yang besar sedangkan tingkat kematiannya kecil. Populasi ini mengalami kenaikkan hingga mencapai titik kestabilannya.

5. Kesimpulan Dalam model matematika pertumbuhan nyamuk ini terdapat 3 model yaitu model pertumbuhan untuk tahapan telur, larva dan nyamuk. Pada model digunakan parameter batasan yang diperoleh dari rata-rata pertumbuhan telur, larva dan nyamuk yaitu =

dimana kondisi parameter batasan ini menunjukkan

keberadaan endemik di suatu populasi. Analisis model menghasilkan dua titik equilibrium yaitu titik equilibrium yang menunjukkan tidak ada pertumbuhan = (0,0,0) dan ada pertumbuhan =

(

)

,

1−

,

.

Hasil simulasi yang diperoleh jika nilai parameter batasannya sangat kecil ( ≤ 1) maka didapat titik equilibrium yang menunjukkan tidak ada pertumbuhan. Artinya, pada model tersebut setelah beberapa waktu, populasi akan habis dan tidak ada endemik sehingga sangat kecil kemungkinan terkena DBD. Sedangkan untuk nilai parameter batasan yang besar ( > 1) didapat titik equilibrium yang menunjukkan adanya 14 pertumbuhan. Artinya, pada model tersebut populasi meningkat dan kemudian laju pertumbuhannya berhenti ketika mendekati titik kesetimbangannya. Pada kondisi ini terdapat keadaan endemik sehingga beresiko tinggi terkena DBD. Pengontrolan nyamuk dapat dilakukan dengan memperkecil nilai parameter batasan. Misalnya mengeringkan tempat bertelur yang dapat meningkatkan kematian telur dan larva sehingga pertumbuhan nyamuk berkurang dan tidak terjadi endemik.

112



6. Daftar Pustaka [1] [2] [3] [4]

[5] [6] [7] [8] [9]

[10]

[11] [12]

[13] [14] [15]

Anton, Howard dan Chris Rorres. 1988 .Penerapan Aljabar Linear.Erlangga. Jakarta. Burden,Richard dan J.Doughlas Faires .2001.Numerical Analysis seventh Edition.Wadsworth Group.United State of America. Christoper.S.R. 1960 .Aedes aegypty The Yellow Fever Mosquito.Cambridge University Press. London Finizio, N., Ladas, G., 1988. Persamaan Diferensial Biasa dengan Penerapan Modern, Edisi kedua, terjemahan Dra. Widiarti Santoso, Jakarta : Erlangga. Ginanjar,genis.2008.Demam Berdarah: A Survival Guide B-First. (PT. Bentang Pustaka) Anggota IKAPI. Gandahusada.S.Herry dan Pribadi W.2000.Parasitologi Kedokteran(ed 3). Fakultas Kedokteran Indonesia.Jakarta. Hamzah, amir.2010.Populasi nyamuk aedes aegypti. Institut Teknolgi Bandung. Huntley, ian & Jhonson Roy Michael. 1983. Linear and Nonlinear Diferential Equation. New York : Elis Horwood Ltd. Moulay,djamila, M.A. Aziz-Alaoui, dan M. Cadivel.2010. The chikungunyah disease: Modeling, vector and transmission global dynamics. Mathematical Biosciences 229 (2011) 50–63. Moulay,djamila, M.A. Aziz-Alaoui, dan Hee-Dae Kwon.2012. Optimal Control of Chikungunya Disease:Larvae Reduction, Treatment and Prevention. Mathematical Biosciences and Engineering Volume 9, Number 2. Mukhsar.2009.Model Pertumbuhan Nyamuk.Universitas Haluoleo: Kendari. Pinho, dll .2010.Modelling the dynamics of danguereal epidemics. Mathematical physical and Engineering science.The philosophical transaction of the royal society. Ross, SL. 1984. Diferential Equation. New York : J.Wiley Rahmawati, Dian. 2004.Perkembangan Nyamuk Aedes Aegypti. IPB: Bogor Soedarmo, SSP.1998.Demam Berdarah (Dangue) pada Anak.Universitas Indonesia: Jakarta


113

Model Pertumbuhan Hidup Nyamuk Aedes Aegypti

Recommend Documents