g ISSN 0216-2393 GRADIEN · PDF fileKimia 4. Penentuan Adanya ... Abstrak - Telah dilakukan penelitian mengenai pemodelan fisis aplikasi metode geolistrik tahanan jenis 2 ... ciri

ISSN 0216-2393

GGRRAADDIIEENN Vol. 2 No. 1 Januari 2006 JURNAL MIPA

FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM UNIVERSITAS BENGKULU

Gradien Vol. 2 No. 1 Hal. 105-151 Bengkulu, Januari 2006

ISSN 0216-2393

g

ISSN 0216-2393

GGRRAADDIIEENN Vol. 2 No. 1 Januari 2006 JURNAL MIPA

DAFTAR ISI

Fisika 1. Pencitraan Konduktivitas Bawah Permukaan dan Aplikasinya untuk Identifikasi

Penyebaran Limbah Cair Dengan Menggunakan Metode Geolistrik Tahanan Jenis 2D (Suhendra) 105-108

2. Analisis Citra Digital Dengan Menggunakan Teknik Penajaman Citra (Arif Ismul Hadi) 109-112

3. Uji Kualitas Akustik Auditorium (V. Sozi Karnefi) 113-115 Kimia 4. Penentuan Adanya Senyawa Triterpenoid Dan Uji Aktivitas Biologis Pada

Beberapa Spesies Tanaman Obat Tradisional Masyarakat Pedesaan Bengkulu (Eni Widiyati). 116-122

5. Studi Kandungan Etanol Dalam Tapai Hasil Fermentasi Beras Ketan Hitam Dan Putih (Teja Dwi Susanto) 123-125

6. Efek Perendaman Larutan Alkali Terhadap Prilaku Film Kertas Dari Nata De Coco Yang Dimodifikasi (Irfan Gustian) 126-129

7. Efek Kenaikan pH Pada Mekanisme Ekstraksi Cair-Cair Terhadap Asam Asam Karboksilat (Agus Martono) 130-133

Matematika 8. Karakterisasi Matrik Leslie Ordo Tiga (Mudin Simanihuruk) 134-137 9. Bentuk Umum Perluasan Teorema Phythagoras (Mulia Astuti) 138-143 10. Metode Diagonally Implicit Runge_Kutta Dalam penyelesaian Model Matematika

Sistem Cardiovascular (Fanani) 144-151

g

Jurnal Gradien Vol.2 No.1 Januari 2006 : 105-108

Pencitraan Konduktivitas Bawah Permukaan dan Aplikasinya untuk Identifikasi Penyebaran Limbah Cair

Dengan Menggunakan Metode Geolistrik Tahanan Jenis 2 D

Suhendra

Jurusan Fisika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Bengkulu, Indonesia

Diterima 5 Desember 2005; disetujui 20 Desember 2005

Abstrak - Telah dilakukan penelitian mengenai pemodelan fisis aplikasi metode geolistrik tahanan jenis 2-D untuk pemantauan masalah lingkungan yaitu pencemaran limbah cair dengan menggunakan metoda geolistrik tahanan jenis memakai konfigurasi elektroda Wenner-Schlumberger. Tujuan penelitian ini untuk mengidentifikasi penyebaran polutan cair (oli bekas) dengan menggunakan metode geolistrik tahanan jenis 2D dan menentukan batas daerah yang terkontaminasi. Pengolahan data diolah dengan menggunakan software RES2Dinv ver 3.3 dan diperoleh nilai tahanan jenis medium yang terkontaminasi oleh limbah cair dari data geolistrik tahanan jenis adalah 2.09 – 4.36 Ωm serta dari penampang (pseudosection) geolistrik terlihat batas atau daerah yang terkontaminasi oleh limbah cair. Kata kunci : Tahanan Jenis; Pencemaran limbah; pemodelan fisis

1. Pendahuluan Metoda geolistrik dapat dimanfaatkan untuk studi masalah lingkungan, yaitu untuk mendeteksi kontras resistivitas medium akibat penyebaran kontaminan (rembesan limbah) di bawah permukaan yang sering diasosiasikan sebagai fluida konduktif [3]. Beberapa studi telah dilakukan diantaranya untuk identifikasi intrusi air laut, kebocoran limbah hasil industri [4] sehingga metoda geolistrik dapat digunakan secara efektif untuk mengidentifikasi distribusi polutan baik secara spasial maupun temporal. Namun untuk keperluan tersebut diperlukan perangkat dan teknik pengukuran yang dapat menghasilkan citra konduktivitas (atau resistivitas) bawah permukaan dengan resolusi tinggi secara tepat dan cepat. Di bawah permukaan sulit untuk mengamati secara langsung rembesan limbah, sejauh mana limbah tersebut akan mengalir dan berapa cepat dia mengalir sehingga perlu dilakukan studi awal untuk penyelidikan masalah lingkungan ini sehingga nantinya akan mudah untuk mencari metode dan cara penangulangan rembesan limbah cair yang akhirnya dapat mengangu keberadaan air tanah yang dikonsumsi secara langsung oleh masyarakat.

Pencemaran limbah dapat dilakukan dengan metoda geolistrik, ini dilakukan dengan menggunakan arus listrik searah yang diinjeksikan melalui dua buah elektroda arus ke dalam bumi, lalu mengamati potensial yang terbentuk melalui dua buah elektroda potensial yang berada di tempat lain. Perbedaan potensial yang terukur merefleksikan keadaan di bawah permukaan bumi. Pada dasarnya metoda ini didekati dengan menggunakan konsep perambatan arus listrik di dalam medium homogen isotropis, dimana arus listrik bergerak kesegala arah dengan nilai yang sama besar. Berdasarkan asumsi tersebut, maka bila terdapat anomali yang membedakan jumlah rapat arus yang mengalir diasumsikan diakibatkan oleh adanya perbedaan akibat anomali tahanan jenis. Anomali ini nantinya digunakan untuk merekonstruksi keadaan geologi di bawah permukaan. Perbedaan konfigurasi elektroda, variasi tahanan jenis spesifik yang diselidiki, prosedur memperoleh data sangat menentukan dalam pemakaian metoda ini.

Suhendra / Jurnal Gradien Vol. 2 No. 1 Januari 2006 : 105-108

106

Metoda tahanan jenis mempunyai dua macam pendekatan, yaitu pendekatan horisontal dan pendekatan vertikal, kedua pendekatan ini mempunyai prosedur kerja dan interpretasi yang berbeda antara satu sama lainnya. Metoda tahanan jenis pendekatan horisontal dimaksudkan sebagai eksplorasi metoda tahanan jenis untuk mendeteksi lapisan atau formasi batuan yang mempunyai kedudukan stratigrafi bidang lapisan yang membentang secara horisontal. Sedang eksplorasi dilakukan untuk mempelajari urutan stratigrafi batas lapisan secara vertikal dari atas sampai bawah.

2. Metode Penelitian Lokasi Penelitian Penelitian ini merupakan penelitian model fisis yang dilakukan di Laboratorium Jurusan Fisika Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Bengkulu, untuk mengamati rembesan limbah di bawah permukaan medium berdasarkan anomali nilai resistivitas medium. Tempat pengambilan data merupakan model fisis berupa bak kaca berukuran 2 m x 1m x 1 m. Teknik Pengambilan Data Data yang diperlukan dalam penelitian ini berupa data tahanan jenis medium berupa tanah lempung (clay sand) dan tahanan jenis medium setelah diinjeksikan dengan limbah cair (Oli bekas), pengukuran di lapangan sebelum dan sesudah diinjeksi dengan limbah untuk mensimulasikan penyebaran polutan cair pada suatu model fisis [2][4] analisis kualitatif dan kuantitatif data lapangan untuk memperoleh penampang tahanan jenis semu melalui pemodelan inversi dengan menggunakan metode yang dikemukan oleh Loke dan Barker [1], untuk mengamati rembesan limbah yang diperoleh malalui pengukuran arus dan potensial listrik dengan menggunakan peralatan resistivitimeter. Dengan adanya data resistivitas dari medium pasir dan medium pasir yang telah diinjeksi kita dapat membandingkan antara medium yang belum terkontaminasi dengan medium yang telah terkontaminasi oleh limbah polutan cair.

Akuisisi data resistivitas dilakukan dengan menggunakan konfigurasi elektroda Wenner-Schlumberger dengan jarak spasi elektroda 5 cm dengan lintasan pengukuran sebanyak satu lintasan yaitu lintasan U-S. Dengan pertimbangan akuisisi ini dapat mengambarkan pola arah aliran secara Vertikal dan horisontal dalam bentuk 2D. secara tepat dan akurat. Alat Ukur Geolistrik Tahanan Jenis Alat yang digunakan untuk penelitian bawah permukaan yaitu Resismeter tipe Naniura Nrd22r beserta kelengkapanya yang terdiri dari accumulator (aki) sebagai sumber arus DC, kabel penghubung, elektroda arus dan potensial serta kelengkapan pendukung seperti meteran. Pengolahan Data Pemantauan variasi terhadap waktu sangat diperlukan dalam penelitian ini karena data tersebut akan menunjukan pola dan arah aliran polutan, pengambaran pola dan arah aliran polutan bawah permukaan dilakukan dengan pemodelan numerik Metoda Finite Elemen. Sebelum data diolah dengan pemodelan numerik metoda Finite elemen terlebih dahulu datum nilai beda potensial dan arus yang terukur di lapangan diolah terlebih dahulu untuk memperoleh nilai tahanan jenis semu (ρa) untuk setiap titik pengukuran dapat dihitung dengan menggunakan persamaan berikut:

ρa = K ∆V / I (1) ρa = Tahanan jenis semu (Ohm. m), K = Koefesien geometri untuk konfigurasi elektroda Wenner-Schlumberger, ∆V = Beda potensial (volt) dan I = Arus listrik (Ampere)

3. Hasil Dan Pembahasan Penelitian ini menggunakan konfigurasi elektroda Wenner-Schlumberger sehingga dapat menjelaskan penyebaran limbah secara lateral dan vertikal sehingga dianalogikan.penyebaran limbah dalam bentuk 2-D. Perhitungan data dengan persamaan 2, diperoleh nilai apparent resistivity pada masing-masing titik ukur. Selanjutnya dilakukan inversi dengan menggunakan


107

program RES2DINV ver 3.3 untuk melihat penampang tahanan jenis dari pasir serta anomaly yang mengidentifikasikan terdapatnya rembesan polutan di bawah permukaan. Berdasarkan hasil inversi dengan menggunakan program RES2DINV ver 3.3 dari konfigurasi elektroda Wenner-Schlumberger dengan jarak spasi elektroda potensial 5 cm, serta panjang lintasan 2 m diperoleh penampang tahanan jeni seperti terlihat pada Gambar 1 berikut ini.

Gambar 1. Penampang hasil inversi 2-D dari Tanah

lempung (clay sand) Tahanan jenis pasir yang diinjeksikan polutan cair Nilai tahanan jenis tanah lempung tanpa limbah yang merupakan pembanding antara medium yang terkontaminasi dengan yang belum terkontaminasi yang merupakan hasil pengukuran arus dan potensial listrik dapat terlihat pada gambar 1. Dari hasil penampang geolistrik pasir tanpa limbah terdapat nilai tahanan jenis berkisar antara 28.6 – 46 Ωm yang diinterpretasikan berupa air, ini dapat diyakini karena bagian tengah-tengah bak kaca pasir terlihat masih basah. Tahanan jenis pasir yang diinjeksikan polutan cair. Nilai tahanan jenis pasir yang diinjeksikan polutan cair yang merupakan hasil pengukuran arus dan potensial listrik dapat dilihat pada Gambar 2. Polutan cair (oli) yang diinjeksikan sebanyak 1000 mL yang ditumpahkan ditengah-tengah bak kaca lalu dibiarkan selama kurang lebih 30 menit untuk membiarkan penyebaran oli benar-benar sudah berhenti sehingga dapat

diukur sejauh mana oli tersebut akan menyebar pada medium yang telah disiapkan.

Gambar 2 Penampang hasil inversi 2-D dari pasir yang

terkontaminasi Dari hasil pengolahan data melalui inversi terlihat bahwa daerah yang sudah terkontaminasi dengan limbah oli yang memiliki nilai resistivitas berkisar dari 2.09 Ω.m sampai dengan 4.36 Ωm, dimana terlihat dari Gambar .2 bahwa limbah terdapat pada kedalaman 5 cm – 20 cm dari permukaan tanah (ground surface). Dari hasil pengolahan data terdapat nilai tahanan jenis yang cukup tinggi yaitu 172 – 359 Ωm ini di interpretasikan merupakan pasir besi yang terdapat pada medium tanah lempung (clay sand) yang digunakan.

4. Kesimpulan Berdasarkan penampang geolistrik tahanan jenis 2-D terlihat adanya polutan limbah cair yang ditunjukkan dengan nilai tahanan jenis yang rendah ini dikarenakan oli bekas bersifat konduktif. Hasil pengamatan menunjukan bahwa nilai resistansi oli memiliki nilai yang cukup rendah yaitu berkisar dari 2.09 Ωm sampai dengan 4.36 Ωm.Aplikasi metode geolistrik dapat dipakai sebagai alat monitoring rembesan limbah, namun untuk kondisi lapangan yang sebenarnya harus masih perlu pengkajian yang lebih mendalam lagi. Serta mutlak harus ditunjang dengan data-data skunder.


108

Daftar Pustaka [1] Loke, M.H., R.D. Barker, Least-square dekonvolusi of

apparent resistivity pseudosection, 1995, Geophysics, Vol 60 pp 1682-1690.

[2] Park, S.K., Fluid migrasi in the vadose zone from 3D inversion of resistivity monitoring data, 1998, Geophysics, 60, 1682-1690.

[3] Telford,W.M., L.P. Geldart and R.E. Sheriff, Applied Geophysic : Second Edition, 1990, Cambrige University Press, USA, 522-538.

[4] Van, G.P., Park, S.K., Hamilton, P., Monitoring leaks from stroge ponds using resistivity methods, 1991, Geophysics, 56, 1267-1270.


Analisis Citra Digital Dengan Menggunakan Teknik Penajaman Citra

Arif Ismul Hadi


Diterima 8 September 2005; disetujui 15 Desember 2005

Abstrak - Analisis citra digital telah dilakukan dengan menggunakan teknik penajaman citra. Proses pengolahan data menggunakan software ILWIS (Integrated Land and Water Information System) versi 1.4 yang merupakan paket program yang memadukan sistem pengolahan citra dengan sistem informasi geografi. Hasil yang diperoleh menunjukkan citra komposit dengan kombinasi FCC 432 adalah yang paling mudah dibedakan. Kata kunci : citra digital; landsat; citra komposit.

1. Pendahuluan Suatu teknik pengumpulan informasi tentang objek dan lingkungannya dari jarak jauh tanpa sentuhan fisik (penginderaan jauh) dapat menghasilkan beberapa bentuk citra untuk selanjutnya diproses dan diinterpretasi, sehingga menghasilkan data yang bermanfaat untuk aplikasi dibidang pertanian, kehutanan, geografi, geologi, dan bidang-bidang lainnya (Lo, C.P., 1986 dalam [5]). Informasi tentang objek disampaikan ke pengamat melalui gelombang elektromagnet. Pada penginderaan jauh, sensor merekam tenaga yang dipantulkan atau dipancarkan oleh objek di permukaan bumi. Objek yang terekam pada citra dikenali berdasarkan ciri-cirinya. Menurut Lillesand dan Kiefer (1983) dalam [2] ciri pengenalan objek ini disebut unsur interpretasi citra (element of image interpretation). Ciri pengenalan benda meliputi ciri spektral, ciri spasial, dan ciri temporal. Di antara tiga ciri ini, maka ciri spektral merupakan ciri utama. Ciri spasial seperti bentuk, ukuran, bayangan, tekstur, pola, situs, dan asosiasi baru tampak kemudian berdasarkan ciri spektralnya. Data penginderaan jauh yang berupa data digital atau data numerik dapat dianalisis dengan menggunakan komputer, sedangkan untuk data visual pada umumnya dianalisis secara manual [3]. Kecerahan rona objek pada citra bergantung atas besarnya tenaga pantulan yang datang dari objek. Berdasarkan

warna yang digunakan, foto berwarna dapat dibedakan atas [3]: foto berwarna semu (false color) atau foto infra merah berwarna dan foto warna asli (true color). Pada foto berwarna semu, warna objek tidak sama dengan warna foto. Objek seperti vegetasi yang berwarna hijau dan banyak memantulkan spektrum infra merah, akan tampak merah (reddish) pada foto, air akan tampak biru (blueish), dan tanah akan tampak hijau (greenish). Sedangkan untuk foto warna asli hasilnya akan tampak sesuai dengan aslinya. Keunggulan foto berwarna semu justru terletak pada warnanya yang tidak serupa dengan warna aslinya. Dengan warna semu ini banyak objek yang pengenalannya akan menjadi lebih mudah. Karena vegetasi yang lebat memantulkan saluran infra merah dekat jauh lebih besar bila dibandingkan dengan pantulannya pada saluran hijau [4], maka ia akan tampak dengan pelbagai tingkat kemerahan pada foto infra merah berwarna dan jika vegetasinya telah mati, maka pada foto infra merah berwarna tidak akan tampak merah. Cara paling umum untuk menonjolkan masing-masing keunggulan secara serentak adalah dengan membuat citra komposit warna [1]. Citra ini menghasilkan paduan saluran 1, 2, dan 3, sehingga menghasilkan informasi seperti pada foto udara infra merah warna semu. Prinsipnya, masing-masing saluran (yang aslinya hanya tersaji dalam rona hitam-putih) diberi warna dasar. Pada

Arif Ismul Hadi / Jurnal Gradien Vol. 2 No. 1 Januari 2006 :109-112

110

umumnya, saluran 3 (infra merah dekat) diberi warna merah, saluran 2 (merah) diberi warna hijau, dan saluran 1 (hijau) diberi warna biru. Kombinasi warna ketiga saluran ini menghasilkan citra berwarna semu, seperti halnya foto udara infra merah berwarna. Cara interpretasi visual untuk komposit ini pun sama dengan cara interpretasi untuk foto udara infra merah berwarna. Di samping teknik penyusunan citra komposit, dikenal pula teknik penajaman citra, yang menghasilkan citra masing-masing saluran dengan kenampakan yang lebih tajam untuk aspek-aspek tertentu. Teknik-teknik ini dapat dikelompokkan menjadi dua, yaitu penajaman kontras (contras enhancement) dan pemfilteran (filtering). Penajaman kontras merupakan teknik manipulasi citra yang mengubah nilai spektral asli menjadi nilai baru, sehingga kekontrasan antar objek menjadi lebih jelas. Pemfilteran merupakan teknik penonjolan (dan sekaligus penghilangan) variasi spektral tertentu, sehingga menghasilkan citra baru yang lebih ekspresif dalam menonjolkan pola-pola tertentu, misalnya kelurusan, jaringan sungai, ataupun jaringan jalan [1]. Dengan penajaman warna ini maka objek pada citra hitam-putih yang semula belum dapat dikenali, kemudian dapat dikenali karena diujudkan dengan warna yang bedanya terhadap objek lain dipertajam. Hasil akhirnya berupa citra paduan warna.

2. Metode Penelitian Dalam penelitian ini, metode yang digunakan adalah metode eksperimen dan visual yaitu dengan teknik penajaman citra. Citra ini merupakan paduan saluran 1, 2, dan 3, sehingga menghasilkan informasi seperti pada foto udara inframerah warna semu. Prinsipnya adalah masing-masing saluran (yang aslinya hanya tersaji dalam rona hitam-putih) di “beri” warna dasar. Pada umumnya, saluran 3 (inframerah dekat) diberi warna merah, saluran 2 (merah) diberi warna hijau. Dan saluran 1 (hijau) diberi warna biru. Kombinasi warna ketiga saluran ini menghasilkan citra berwarna semu.Vegerasi berwarna merah, air berwarna biru gelap-biru kehijauan, dan lahan terbuka berwarna biru kehijauan-putih cerah [1].

3. Hasil Dan Pembahasan Pengenalan objek merupakan bagian vital dalam interpretasi citra. Tanpa dikenali identitas dan jenis objek yang tergambar pada citra, tidak mungkin dilakukan analisis untuk memecahkan masalah yang sedang dihadapi. Prinsip pengenalan objek pada citra mendasarkan atas penyelidikan karakteristiknya atau atributnya pada citra. Karakteristik objek yang tergambar pada citra dan digunakan untuk mengenali objek disebut interpretasi citra [3].

FCC 432

FCC 532

FCC 732

TCC 342

TCC 372


111

TCC 352

Tabel 1. Interpetasi

Citra Landsat Komposit Objek FCC

432 FCC

532 FCC

732 TCC

342 TCC

352 TCC

372 Air *** *** *** * * *

Tanah *** ** ** * * *

Vegetasi *** *** ** *** *** **

Sungai *** *** *** *** *** ***

Tambak *** *** *** *** *** ***

Tanaman pantai

*** ** * *** ** *

Kerapatan vegetasi

*** ** * *** *** ***

Tk. Kekeringan tanah

*** ** * * * *

Keterangan : *** Sangat jelas (mudah dibedakan) ** Sedang

* Tidak jelas (sulit dibedakan)

Dari tabel 1 di atas dapat dijelaskan bahwa untuk objek air akan tampak sangat jelas atau mudah dibedakan pada citra komposit FCC (False Color Composite) 432, FCC 532, dan FCC 732 yaitu berwarna biru. Sedangkan pada TCC (True Color Composite) 342, TCC 352, dan TCC 372 tidak jelas atau sulit dibedakan. Hal ini terjadi karena antara rona air dan rona tanah hampir mempunyai warna yang sama (pink gelap dan pink terang).

Untuk objek tanah, citra komposit FCC 432 tampak sangat jelas/mudah dibedakan yaitu berwarna hijau yang tampak jelas. Pada FCC 532 dan FCC 732 tanah tampak sedang/kurang jelas walaupun masih tampak hijau sedikit. Sedangkan tanah pada TCC 342, TCC, 352, dan TCC 372 sulit dibedakan, karena warnanya hampir sama dengan warna air yaitu warna pink.

Untuk objek vegetasi, citra komposit FCC 432 dan FCC 532 tampak sangat jelas dengan rona merah dan masih tampak antara kerapatan vegetasinya (rona merah muda dan merah tua), untuk citra komposit FCC 732 agak kurang tampak, karena warnanya hampir sama dengan warna tanah yaitu warna coklat tua dan coklat muda. Sedangkan vegetasi pada TCC 342 dan TCC 352 tampak sangat jelas dengan rona hijau dan pada TCC 372 tampak kurang jelas, karena rona tanah juga tampak hijau muda. Untuk objek sungai pada semua citra komposit yang ada tampak sangat jelas/mudah dibedakan. Hal ini terjadi karena aliran sungai tampak memanjang munuju ke arah laut. Untuk objek tambak juga tampak sangat jelas/mudah dibedakan pada semua citra komposit yang ada seperti tampak pada sebelah kanan atas (pinggir laut) dari gambar-gambar tersebut. Untuk objek tanaman pantai, citra komposit FCC 432 dan TCC 342 tampak sangat jelas/mudah dibedakan dengan rona merah pada FCC 432 dan hijau pada TCC 342 di pinggir laut/pantai. Pada FCC 532 dan TCC 352 tampak kurang jelas, walaupun masih tampak ada rona merah pada FCC 532 dan hijau untuk TCC 352 di pantai. Dan untuk FCC 732 dan TCC 372 sulit dibedakan. Untuk kerapatan vegetasi tampak jelas pada citra komposit FCC 432, TCC 342, TCC 352, dan TCC 372. Hal ini karena masing-masing citra komposit tersebut masih terlihat rona merah tua dan mudanya untuk FCC 432, dan rona hijau tua dan muda untuk semua citra komposit TCC. Sedangkan untuk FCC 532 tampak kurang jelas, walaupun masih sedikit terlihat perbedaan antara rona merah tua dan muda, sedangkan untuk FCC 732 sulit dibedakan. Untuk tingkat kekeringan tanah, pada FCC 432 tampak sangat jelas dengan rona hijau muda. Untuk FCC 532 tampak kurang jelas perbedaan warna hijaunya, walaupun masih tampak sedikit, sedangkan untuk FCC 732, TCC 342, TCC 352, dan TCC 372 sulit dibedakan.


112

4. Kesimpulan Dari hasil penelitian dan pembahasan di atas dapat disimpulkan : 1. Pada citra komposit FCC (False Color Composite)

objek vegetasi akan berwarna merah, objek air akan berwarna biru dan objek tanah akan berwarna hijau.

2. Dari semua citra komposit yang telah ditampilkan, maka citra komposit dengan kombinasi FCC 432 adalah yang paling mudah untuk diinterpretasi. Karena rona yang ditampilkan mudah dibedakan (seperti terlihat pada tabel interpretasi).

Ucapan Terima Kasih

Pada kesempatan ini penulis ucapkan banyak terima kasih kepada Bapak Drs. Heru, M.Sc atas diskusinya dalam pengolahan citra digital di Lab. Penginderaan Jauh, UGM.

Daftar Pustaka [1] Danoedoro, Petunjuk Praktikum Pengolahan Citra Digital

dengan ILWIS versi 1.4, 2001, Fakultas Geografi UGM, Yogyakarta.

[2] Sutanto, Penginderaan Jauh Dasar, 1995, Fakultas Geografi UGM, Yogyakarta.

[3] Sutanto, Penginderaan Jauh Jilid 1, 1986, Gadjah Mada University Press, Yogyakarta.

[4] Sutanto, Penginderaan Jauh Jilid 2, 1987, Gadjah Mada University Press, Yogyakarta.

[5] Purbowaseso, B., Penginderaan Jauh Terapan, 1996, UI-Press. Jakarta.


Uji Kualitas Akustik Auditorium

V. Sozi Karnefi



Abstrak. - Pentingnya ilmu akustik dalam disain arsitektur masa kini merupakan suatu yang tidak dapat diingkari lagi. Rancang bangun arsitektur terasa kurang sempurna tanpa memperhatikan pentingnya keberadaan ilmu akustik, dalam menciptakan keadaan mendengar yang ideal. Oleh karena itu pengendalian bunyi secara arsitektural dibagi dalam dua sasaran utama yaitu akustik ruang dan pengendalian bising. Dengan demikian keberadaan akustik ruang dan pengendalian bising akan memungkinkan terciptanya kesempurnaan rancang bangun, baik secara akustik maupun arsitektural. Kata kunci : Waktu Reverbrasi; Akustik Ruang.

1. Pendahuluan

Pengembangan arsitektur masa kini berfungsi lebih dari sekedar menyediakan ruang yang memerlukan aspek keindahan, fungsional teknik, maupun ekonomi tetapi juga sebagai perlindungan bagi penghuninya terhadap pengaruh termal, cahaya dan bunyi yang berasal dari lingkungan untuk mengharapkan kenyamanan, keamanan, lingkungan yang menyenangkan, penerangan yang cukup, pemandangan yang mengesankan serta bunyi yang laik didengar disamping kualitas acara itu sendiri [4]. Hal-hal yang menentukan kondisi mendengar dalam suatu ruang – dan merupakan pertimbangan arsitektur murni – seperti bentuk ruang, dimensi dan volume, letak batas-batas permukaan, pengaturan tempat duduk, lapisan permukaan dan bahan-bahan untuk dekorasi interior, sedikit banyak akan menentukan penampilan akustik ruang tersebut. Oleh karena itu pengendalian bunyi secara arsitektural mempunyai dua sasaran yaitu 1. Akustik Ruang

Akustik ruang bertujuan menyediakan keadaan yang paling diinginkan untuk menghasilkan perambatan dan

penerimaan bunyi dalam ruang yang digunakan untuk berbagai tujuan mendengar seperti auditorium.

2. Pengendalian Bising Ruang lingkup pengendalian menyangkut peniadaan atau pengurangan bising dalam jumlah yang memadai.

Suatu ruang dengan waktu reverbrasi yang tidak memenuhi tuntutan kualitas akustik akan merusak kondisi akustik ruang, padahal waktu reverbrasi merupakan parameter akustik ruang yang sangat sering digunakan dalam mendisain raung tertutup, disamping pentingnya menghilangkan cacat-cacat akustik yang terjadi dalam ruang dengan jalan menciptakan penyebaran bunyi. Persamaan Waktu Reverbrasi Sabine Pentingnya pengedalian dengung dalam rancangan akustik auditorium telah mengharuskan masuknya besaran standar yang relevan yaitu waktu reverbrasi T, yaitu lamanya waktu yang dibutuhkan untuk bunyi meluruh 60 dB. Sabine adalah orang pertama yang membentuk hubungan kuantitatif antara waktu reverbrasi, volume ruang dan jumlah penyerapan total ruang.

V. Sozi Karnefi / Jurnal Gradien Vol. 2 No. 1 Januari 2006 : 113-115

114

Persamaan Sabine menyediakan pendekatan yang beralasan untuk perhitungan waktu reverbrasi yang memungkinkan arsitek melakukan pengendalian dengan mudah dan cepat untuk rancangan akustik suatu auditorium. Waktu reverbrasi akan berubah dengan perubahan frekuensi, dikarenakan penyerapan bermacam-macam bahan dan lapisan penutup yang digunakan dalam rancangan auditorium kadang-kadang berubah cukup banyak dengan perubahan frekuensi, maka sangat penting untuk menyatakan atau menghitung waktu reverbrasi untuk sejumlah frekuensi wakil pada jangkauan frekuensi audio misalnya frekuensi 125, 250, 500, 1000, 2000, dan 4000. Jumlah frekuensi wakil ini tergantung pada pentingnya hal-hal yang berkaitan. Bila suatu waktu reverbrasi disebutkan tanpa menyertakan suatu frekuensi tertentu, maka reverbrasinya adalah untuk frekuensi tengah 500 Hz. Secara matematis persamaan yang memberikan peluruhan bunyi difuse uniform diperoleh dengan mengambil harga

0 =π . Jika sumber dimatikan pada saat t = 0, tekanan bunyi kemudian menjadi [5],

E -t /2r

2r e (0) P P τ=

Sehingga perubahan tingkat tekanan bunyi – sound preassure level yang disingkat SPL – menjadi,

E t/4,35 SPL τ=∆ Sedangkan waktu reverbrasi didefenisikan sebagai waktu yang dibutuhkan untuk menurunkan tingkat tekanan bunyi 60 dB, maka didaptkan harga T,

cA V 55,2 E 13,8 T == τ

Diketahui suhu rata-rata untuk daerah Pekanbaru adalah 30oC maka kecepatan bunyi untuk daerah Pekanbaru dapat diperoleh dengan menggunakan persamaan konversi berikut,

c = 331,6 + 0,6 T Dengan c adalah keepatan bunyi dan T adalah suhu dalam oC. Maka didapat harga kecepatan bunyi untuk daerah Pekanbaru adalah c = 349,6 m/dt, sehingga rumus waktu reverbrasi untuk daerah Pekanbaru menjadi,

α SV 0,158 T =

Dengan T adalah waktu reverbrasi, α adalah koefisien penyerapan, V adalah volume ruang dan S adalah luas permukaan

2. Metode Penelitian 1. Atur dan susunlah peralatan. 2. Setelah generator attenuator pada – 40 dB. 3. Putarlah counter clockwise pada posisi maksimum [2]. 4. Setelah power menjadi on yang akan ditunjukan dengan

menyalakan lampu CAL 10 detik kemudian. 5. Biarkan posisi recorder push off (jika ingin

menggunakan recorder tekan push on) tombol 16. 6. Tekanlah saklar SYNC pada INT untuk pengaturan

peralatan, nomor tombil 14. 7. Tekanlah tombol MIC untuk penggunaan input jack

nomor tombol 12 8. Aturlah saklar SCALE pada posisi 50 dB. Jika

menggunakan recorder setelah kertas grafik dengan meletakan IC pen pada titik awal. Bila REC MODE diset pada GRAPH ujung pen harus diletakan pada tanda ∆ (45 Hz).

9. Aturlah saklar frekuensi pada posisi 500 Hz tombol 21. 10. Warble swing diset pada lebar 1 oktaf. Tombol 26. 11. FS selector switch diset pada posisi 1 sec, untuk

pembacaan skala bagian atas pada reverbrasi meter, tombol 30.

12. FUNCTION diset pada REV, tombol 31. 13. REC MODE diset pada switch yang lain, lalu waktu

reverbrasi dapat dibaca pada meter reverbrasi tombol 34.

14. Aturlah control level (tombol 8) pada INDIKATOR LEVEL dan KONTROL LEVEL (tombol 24) pada SIGNAL GENERATOR sehingga METER REVERBRASI menunjukan titik CAL, setelah itu LEVEL SIGNAL GENERATOR menjadi minimum. Jika pada saat itu, METER REVERBRASI gagal menunjukan daerah noise, maka putarlah SWITCH LEVEL pada LEVEL INDIKATOR COUNTER CLOCKWISE untuk CAL pada keadaan yang cocok (tingkat noise di luar daerah noise yang akan diukur akibatnya akan terjadi salah pengukuran).

V. Sozi Karnefi / Jurnal Gradien Vol. 2 No. 1 Januari 2006 : 113-115

115

15. Periksa kembali apakah lampu CAL telah hidup, maka tekan tombol START, signal akan keluar yang ditunjukan dengan hidupnya lampu MEASURE, bila signal telah keluar sebanyak lima kali maka lampu MEASURE akan mati, METER REVERBRASI menunjukan periode yang telah ditentukan oleh HOLD time, kemudian lampu REV akan hidup bila pengukuran tidak sesuai dengan noise ruang.

16. Jika lampu CAL hidup kembali, setelah frekuensi menjadi 125 Hz maka lakukan lagi langkah ke 15. putarlah frekuensi yang dibutuhkan yaitu 500 Hz

3. Hasil dan Pembahasan

Hasil yang diperoleh untuk beberapa macam auditorium disajikan pada tabel 1 di bawah ini:

Tabel 1. Data auditorium yang diamati

Waktu Reverbrasi No Auditorium Teoritis

(detik) Pengukuran

(detik) Ideal

(detik) 1. 2. 3. 4.

Ruang Pertemuan PKM Mesjid Ruang Kuliah

1,73 2,87 1,00 0,44

1,6 2,9 0,95 0,41

1,43 0,9 0,6 0,1

Dari hasil pengukuran pertimbangan arsitektur murni dan manipulasi matematis diperoleh bahwa waktu reverbrasi untuk ruang PKM adalah 1,72 detik untuk Mesjid Akramunas 2,87 detik untuk ruang kuliah 1 detik dan untuk studio 0,44 detik. Sedangkan waktu reverbrasi ideal untuk ruang PKM adalah 1,43 detik untuk Mesjid Akramunas 0,9 detik, untuk ruang kuliah 0,6 detik dan untuk studio 0,1 detik. Dari hasil pengukuran diperoleh juga harga waktu reverbrasi ruang kuliah 0,95 detik, untuk studio 0,41 detik dan untuk PKM 1,6 detik, * maka didapat kesalahan antara harga waktu reverbrasi teoritis terhadap pengukuran untuk PKM 7,5%, untuk ruang kuliah 7,3% dan untuk studio 5,3%, sehingga harga waktu reverbrasi yang dihitung dapat dipercaya.

Sedangkan penyimpangan harga waktu reverbrasi teoritis terhadap harga ideal untuk PKM 20%, untuk Mesjid Akramunas 309%, untuk ruang kuliah 100% dan untuk ruang studio 340%. Diperoleh juga penyimpangan harga waktu reverbrasi pengukuran terhadap harga ideal untuk PKM 12%, untuk ruang kuliah 50% dan untuk ruang stuio 310%. Dapat dilihat penyimpangan harga waktu reverbrasi ruang-ruang diatas sangat besar dikarenakan kesalahan desain ruang hal ini menunjukan bahwa disain arsitekturnya tidak memperhatiakn pentingnya keberadaan akustik.

4. Kesimpulan

1. Satu parameter yang penting dalam desain ruang adalah waktu reverbrasi.

2. Akustik ruang bertujuan menyediakan keadaan yang paling diinginkan untuk menghasilkan rambatan dan penerimaan bunyi dalam ruang yang digunakan untuk berbagai tujaun mendengar.

3. Ruang lingkup pengendalian bising menyangkut peniadaan atau pengurangan bising.

Daftar Pustaka

[1] Anonim, Acoustic Measurement, 1962, John Wiley & Sons,

Inc, new York,. [2] Anonim, Instruction Manual, Sound and Acoustic

Meansurement, Equipment, SE 3000, Trio Ken Wood. [3] Baranek Leo.L, Noise and Reduction, 1960, McGraw-Hill

Book Company, Inc, New York. [4] Doelle.L, Environmental Acoustic, Terjemahan Lirik Prasetio,

1972, McGraw-Hill, Inc. [5] Kinsler.L.E., Fundamental of acoustic, 1982, John Welley

and Sons, Third edition. [6] M.Morse.P, Vibration and Sound, Second edition, 1948,

McGraw-Hill Book Company, Inc, Tokyo. [7] Thomson.W.T, Theory of Vibration With applicatons, Second

edition, 1981, Prenctice-Hall, Inc.


Penentuan Adanya Senyawa Triterpenoid Dan Uji Aktivitas Biologis Pada Beberapa Spesies

Tanaman Obat Tradisional Masyarakat Pedesaan Bengkulu

Eni Widiyati

Jurusan Kimia, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Bengkulu, Indonesia


Abstrak - Telah dilakukan penelitian tentang penentuan adanya senyawa Triterpenoid dan uji aktivitas biologis pada beberapa spesies tanaman obat tradisional masyarakat pedesaan Bengkulu. Penelitian ini terdiri dari dua tahap yaitu tahap satu penentuan adanya senyawa Triterpenoid pada bagian tanaman dengan menggunakan pereaksi Liebermann-Burchard dan tahap dua, uji aktivitas biologis (bioassay) dengan uji Brine Shrimp menggunakan hewan uji Arthemia salina Leach terhadap ekstrak bagian tanaman yang mengandung banyak dan sangat banyak senyawa Triterpenoid. Dari penelitian yang dilakukan telah diketahui sebanyak 34 spesies tanaman obat tradisional masyarakat pedesaan Bengkulu mengandung senyawa Triterpenoid dan tidak semua bagian tanaman mengandung senyawa Triterpenoid. Hasil uji aktifitas biologis menunjukkan bahwa dari 8 ekstrak bagian tanaman yang diuji, diketahui semuanya mempunyai harga LC 50 < 1000 ppm, berarti ekstrak kasar yang mengandung Triterpenoid tersebut mempunyai aktivitas biologis. Kata kunci: Triterpenoid; Liebermann-Burchard; bioassy

1. Pendahuluan Saat ini dunia farmasi dan kedokteran telah berkembang pesat, sehingga sudah banyak dibuat dan dipakai berbagai jenis obat-obatan yang diproduksi oleh pabrik-pabrik farmasi. Oleh karena sebagian besar bahan baku untuk pembuatan obat-obatan tersebut masih diimport dari luar negeri, maka mengakibatkan harganya menjadi mahal dan kadang-kadang tidak terjangkau oleh sebagian masyarakat, khususnya masyarakat pedesaan. Adanya krisis moneter dan krisis ekonomi yang berkepanjangan di Indonesia akhir-akhir ini membuat harga obat-obatan produksi pabrik menjadi semakin mahal dan semakin tidak terjangkau lagi oleh masyarakat kecil. Untuk mengatasi permasalahan tersebut maka perlu digalakkan penggunaan obat-obatan tradisional, khususnya yang dibuat dari ramuan tanaman. Tanaman, khususnya tanaman obat tradisional mudah diperoleh karena dapat ditanam sendiri di pekarangan rumah, selain itu tanaman (tumbuhan) merupakan sumber daya alam yang dapat diperbaharui sehingga tidak akan mengalami kepunahan apabila dilestarikan.

Tanaman dikenal banyak mengandung senyawa-senyawa kimia khususnya senyawa metabolit sekunder. Salah satu senyawa metabolit sekunder yang terkandung di dalam tanaman adalah senyawa Triterpenoid. Senyawa tersebut dapat dijumpai pada bagian akar, batang, daun, buah maupun biji tanaman.

Triterpenoid adalah senyawa metabolid sekunder yang kerangka karbonnya berasal dari enam satuan isoprena dan diturunkan dari hidrokarbon C 30 asiklik , yaitu skualena. Senyawa ini berbentuk siklik atau asiklik dan sering memiliki gugus alkohol, aldehida, atau asam karboksilat [1]. Sebagian besar senyawa Triterpenoid mempunyai kegiatan fisiologi yang menonjol sehingga dalam kehidupan sehari-hari banyak dipergunakan sebagai obat seperti untuk pengobatan penyakit diabetes, gangguan menstuasi, patukan ular, gangguan kulit, kerusakan hati dan malaria. Sedang bagi tumbuhan yang mengandung senyawa Triterpenoid terdapat nilai ekologi karena senyawa ini bekerja sebagai anti fungus, insektisida, anti pemangsa, anti bakteri dan anti virus [2]. Uji kimia yang dapat dilakukan untuk mengetahui adanya senyawa

Eni Widiyati / Jurnal Gradien Vol. 2 No. 1 Januari 2006 :116-122

117

Triterpenoid dalam bagian tumbuhan adalah dengan menggunakan pereaksi Liebermann- Burchard [1], sedangkan untuk mengetahui adanya keaktifan biologis dari ekstrak bagian tanaman yang mengandung senyawa Triterpenoid dapat dilakukan dengan uji Brine Shrimp menggunakan hewan uji Arthemia Salina Leach [3]. Sejak dahulu, masyarakat pedesaan Propinsi Bengkulu sudah menggunakan obat-obatan tradisional untuk mengatasi berbagai penyakit. Meskipun perkembangan sistem dan cara pengobatan moderen telah maju pesat, namun pengobatan tradisional masih tetap mendapat tempat di dalam masyarakat pedesaan. Obat-obatan tradisional sebagai sarana penyembuhan berbagai penyakit, telah diwarisi secara turun-temurun dan telah menunjukkan bukti keampuhannya yang dapat diakui keberadaan dan kebenarannya sampai sekarang [4]. Sampai saat ini, belum banyak informasi tentang adanya kandungan senyawa metabolit sekunder khususnya Triterpenoid pada tanaman yang dipergunakan sebagai ramuan obat tradisional masyarakat pedesaan Bengkulu. Oleh karena itu perlu ditingkatkan penelitian-penelitian di bidang kimia, khususnya kimia organik bahan alam, sehingga akan memperkaya informasi di bidang tersebut. Informasi tentang adanya tanaman yang mengandung senyawa Triterpenoid serta adanya keaktifan biologis dari beberapa tanaman obat tradisional masyarakat pedesaan Bengkulu dapat memberikan sumbangan yang positif di bidang dunia kesehatan. Selain itu informasi tersebut dapat memberikan data dasar bagi peneliti selanjutnya yang tertarik untuk melakukan penelitian di bidang senyawa Triterpenoid. Berdasarkan pada latar belakang di atas maka penelitian ini bertujuan untuk mengetahui adanya senyawa Triterpenoid dan menentukan adanya aktivitas biologis beberapa tanaman yang mengandung banyak senyawa Triterpenoid pada beberapa spesies tanaman obat tradisional masyarakat pedesaan Propinsi Bengkulu.

3. Metode Penelitian Dalam penelitian ini digunakan alat dan bahan sebagai berikut:

1. Alat-alat yang digunakan pada penelitian ini meliputi pisau potong, lumpang porselin dan penumbuk, pipet tetes, pelat tetes, tabung reaksi, gelas ukur, tempat penetasan Arthemia salina, seperangkat alat ekstraksi sokhlet dan labu takar.

2. Bahan-bahan yang digunakan pada penelitian ini meliputi sampel tanaman obat, kloroform, asam sulfat pekat, asam asetat anhidrat, Brine shrimp egg, DMSO, metanol, akuades, dan garam.

Sedangkan metode penelitian meliputi: a. Penentuan Adanya Senyawa Triterpenoid Pada

Bagian Tanaman Obat Tradisional. Pada penelitian ini dikumpulkan beberapa spesies tanaman obat tradisional yang biasa dipergunakan oleh masyarakat pedesaan Propinsi Bengkulu. Bagian yang diambil meliputi akar, batang , daun, bunga dan buah. Masing-masing spesies dimasukkan ke dalam kantong plastik yang berbeda, kemudian dilakukan identifikasi adanya senyawa Triterpenoid sebagai berikut : Sebanyak 4 gram sampel segar (bagian tanaman dapat berupa akar, daun, kulit batang, bunga atau buah) digerus dalam lumpang porselin dan dilarutkan dengan kloroform sebanyak 10 ml. Filtrat yang dihasilkan dipindahkan ke dalam pelat tetes dan dibiarkan sampai menguap pelarutnya. Kemudian ke dalamnya ditambahkan pereaksi Liebermann-Burchard (1 tetes asam sulfat pekat dan 3 tetes asam asetat anhidrat). Uji positif jika timbul warna merah, merah jambu atau ungu. Untuk mengetahui banyak sedikitnya kandungan Triterpenoid dalam bagian tanaman, digunakan biji mahoni sebagai pembanding yang telah diketahui banyak mengandung senyawa riterpenoid (+++). Apabila sampel menunjukkan warna yang lebih pekat dari pembanding diberi tanda (++++), bila intensitas warnanya lebih kecil dari pembanding diberi tanda (++) atau (+) dan sampel yang tidak menunjukkan adanya Triterpenoid diberi tanda (-) [5]. b. Uji Aktivitas biologis (bioassy) dengan uji Brine

Shrimp Untuk keperluan bioassay, maka bagian tanaman yang banyak mengandung Triterpenoid diekstrak dengan pelarut metanol. Ekstrak yang diperoleh diuapkan


118

pelarutnya sehingga diperoleh ekstrak kental yang siap untuk bioassay. Adapun pelaksanaan bioassay adalah sebagai berikut : sebanyak 500 ml larutan garam dimasukkan ke dalam wadah sebagai tempat penetasan yang dibuat terdiri dari dua bagian yaitu separo bagian tertutup dan separo bagian lain dibiarkan terbuka. Dimasukkan ke dalamnya Brine shrimp eggs secukupnya, kemudian ditempatkan di bawah sinar lampu yang diberi airasi. Setelah telur menetas, larva akan bergerak bebas. Ditunggu selama 2 hari sampai tumbuh menjadi larva dewasa. Kemudian disiapkan larutan sampel yaitu larutan ekstrak bagian tanaman yang mengandung Triterpenoid, dengan konsentrasi 10 ppm, 100 ppm, dan 1000 ppm. Ke dalam 3 buah tabung reaksi masing-masing dimasukkan larutan sampel dengan konsentrasi yang berbeda yaitu 10 ppm, 100 ppm dan 1000 ppm serta satu tabung (tabung ke empat) diisi larutan blanko sebagai kontrol. Pada setiap tabung yang telah diisi larutan sampel dan larutan blanko di atas dimasukkan sebanyak 10 individu larva shrimp. Pengamatan dilakukan setelah 24 jam untuk mengetahui adanya larva shrimp yang mati. Percobaan tersebut diulangi sebanyak tiga kali. Data yang diperoleh kemudian dihitung harga LC50 nya dengan menggunakan program Finney’s Probit Analysis [6] atau analisis probit dan persamaan regresi linier [7] . Apabila harga LC 50 ≤ 1000 ppm, berarti ekstrak sampel yang dianalisis mempunyai aktivitas biologis [2].

4. Hasil Dan Pembahasan 1. Penentuan Adanya Senyawa Triterpenoid Pada

beberapa Spesies Tanaman Obat Tradisional Masyarakat Pedesaan Bengkulu

Untuk mengetahui adanya senyawa Triterpenoid pada tanaman dapat dilakukan dengan menambahkan pereaksi Lieberman-Burchard yang terdiri dari asam sulfat pekat dan asam asetat anhidrat. Jika pada bagian tanaman yang dianalisis mengandung senyawa Triterpenoid, maka ekstrak bagian tanaman yang diuji menunjukkan terjadi perubahan warna yaitu warna merah, merah jambu atau ungu. Untuk mengetahui seberapa banyak kandungan Triterpenoid yang terdapat pada bagian tanaman, maka intensitas warna yang ditimbulkan pada penambahan

pereaksi Liebermann-Burchard dibandingkan dengan intensitas warna yang terjadi pada ekstrak biji mahoni yang ditambah pereaksi Liebermann-Burchard. Digunakan biji mahoni sebagai pembanding karena di dalamnya terkandung banyak senyawa Triterpenoid yang diberi tanda (+++). Hasil penelitian menunjukkan, intensitas warna yang dihasilkan oleh beberapa bagian tanaman berbeda-beda, hal ini berarti kandungan Triterpenoid dalam bagian tanaman juga bervariasi, mulai dari (+) berarti sedikit, (++) berarti cukup banyak, (+++) berarti banyak dan (++++) sangat banyak mengandung senyawa Triterpenoid.

Dari hasil penelitian tersebut juga dapat diketahui ada beberapa bagian tanaman yang memberikan hasil negatif, yang berarti tidak semua bagian tanaman mengandung senyawa Triterpenoid. Hal ini sesuai dengan teori yang mengatakan bahwa senyawa metabolit sekunder, salah satunya adalah senyawa Triterpenoid tidak disebarkan secara universal ke seluruh bagian tanaman [1]. Sampai saat ini telah dilakukan analisis untuk menentukan adanya senyawa Triterpenoid terhadap beberapa spesies tanaman obat tradisional masyarakat pedesaan Bengkulu. Hasil penelitian menunjukkan 34 spesies tanaman positif mengandung senyawa Triterpenoid. Daftar beberapa spesies tanaman yang mengandung senyawa Triterpenoid terdapat di tabel 1. Dari tabel 1 terlihat sebagian besar bagian tanaman yang biasa digunakan oleh masyarakat pedesaan Bengkulu untuk mengobati jenis penyakit tertentu ternyata mengandung senyawa Triterpenoid. Kandungan Triterpenoid yang terdapat pada bagian tanaman obat tradisional yang digunakan juga bervariasi yaitu ada yang sedikit (+), cukup banyak (++), banyak (+++) dan sangat banyak (++++). Namun ada juga bagian tanaman yang tidak mengandung senyawa Triterpenoid tetapi sering digunakan oleh masyarakat untuk pengobatan tradisional. Kemungkinan, bagian tanaman tersebut mengandung senyawa metabolit sekunder selain Triterpenoid seperti alkaloid, steroid atau flavonoid. Ada juga bagian tanaman yang banyak mengandung senyawa Triterpenoid namun belum dimanfaatkan untuk pengobatan tradisional. Data tabel 1 tersebut dapat berubah dengan bertambahnya penelitian-penelitian tentang penentuan adanya senyawa


119

Triterpenoid pada tanaman obat tradisional masyarakat pedesaan Bengkulu.

2. Uji Aktivitas Biologis Terhadap Ekstrak Bagian

Tanaman Yang Mengandung Banyak Senyawa Triterpenoid

Pada penelitian ini, bioassay (uji aktivitas biologis) dilakukan pada bagian tanaman yang mengandung kadar Triterpenoid banyak (+++) dan sangat banyak (++++). Alasan pemilihan bagian tanaman yang mengandung Triterpenoid banyak dan sangat banyak untuk bioassay adalah agar setelah bagian tanaman tersebut diekstrak, maka diharapkan kandungan terbesar dari ekstrak kasar tersebut adalah Triterpenoid, sehingga pada saat bioassy, diharapkan senyawa yang mempunyai aktivitas biologis adalah senyawa Triterpenoid tersebut. Namun dugaan tersebut masih harus dibuktikan dengan jalan mengisolasi Triterpenoid yang terdapat pada ekstrak bagian tanaman tersebut kemudian dilakukan bioassay terhadap Triterpenoid murni yang dihasilkan. Bagian tanaman yang mengandung banyak dan sangat banyak senyawa Triterpenoid tersebut diekstrak dengan alat ekstraksi sokhlet menggunakan pelarut metanol. Dipilihnya alat ekstraksi sokhlet karena alat tersebut mempunyai beberapa keuntungan antara lain dapat mengekstrak sampel dalam jumlah banyak (skala gram), penyarian dapat dilakukan berulang-ulang sehingga semua ekstrak dapat terambil, dan pelarut yang digunakan dapat digunakan kembali karena tidak menguap. Pada penelitian ini untuk mengetahui adanya aktivitas biologis dari ekstrak tanaman yang mengandung Triterpenoid digunakan uji Brine Shrimp. Metode ini dipilih karena dapat digunakan untuk mengetahui adanya efek sitotoksik dan juga untuk memperoleh hewan uji lebih mudah, harganya murah, telurnya dapat tahan beberapa tahun bila disimpan ditempat kering, mengerjakannya sederhana dan lebih cepat. Disamping itu metode ini telah diuji dan mempunyai korelasi positif dengan metode yang telah biasa digunakan untuk penapisan senyawa anti kanker. Jika ekstrak bagian tanaman yang diuji mempunyai harga LC50 < 1000 ppm maka ekstrak tersebut dinyatakan mempunyai aktivitas biologis yaitu mengandung senyawa yang bersifat sitotoksik [2]. Hasil bioassay pada penelitian ini terdapat di tabel 2.

Data tabel 2 menunjukkan bahwa dari 8 sampel yang mengandung banyak dan sangat banyak Triterpenoid yang dianalisis semua mempunyai harga LC 50 < 1000. Berarti pada konsentrasi tersebut terdapat 50 % individu Artemia salina Leach ( Brine shrimp) yang mati. Hal ini berarti ke 8 ekstrak bagian tanaman tersebut mempunyai aktivitas biologis terhadap Arthemia salina Leach atau mengandung senyawa yang bersifat sitotoksik. Namun untuk mengetahui apakah senyawa yang mempunyai aktivitas tersebut adalah Triterpenoid, perlu dilakukan penelitian lanjutan, yaitu dengan mengisolasi Triterpenoid yang terdapat pada ekstrak tersebut. Setelah dihasilkan Triterpenoid murni maka diuji kembali dengan uji Brine shrimp. Perlu juga dilakukan ekstraksi menggunakan pelarut selain metanol seperti kloroform atau karbon tetra klorida untuk mengetahui kemungkinan adanya senyawa Triterpenoid yang larut dalam pelarut yang non polar atau kurang polar. Dengan diketahui adanya aktivitas biologis pada bagian tanaman yang mengandung banyak senyawa Triterpenoid maka penggunakan bagian tanaman tersebut sebagai obat tradisional perlu lebih digalakkan.

4. Kesimpulan

Dari penelitian yang telah dilakukan dapat diambil kesimpulan sebagai berikut : a. Dari beberapa spesies tanaman obat tradisional yang

dianalisis diperoleh 34 spesies tanaman mengandung senyawa Triterpenoid dan tidak semua bagian tanaman yang dianalisis mengandung Triterpenoid.

b. Dari 8 ekstrak bagian tanaman yang mengandung banyak dan sangat banyak senyawa Triterpenoid diketahui semuanya mempunyai aktivitas terhadap Brine shrimp berarti ekstrak kasar ke 8 bagian tanaman tersebut mempunyai aktivitas biologis.

Mengingat masih banyak tanaman obat tradisional yang belum dianalisis maka perlu dilakukan penelitian tentang identifikasi adanya Triterpenoid dan bioassay dengan sampel atau spesies tanaman obat tradisional masyarakat pedesaan Bengkulu yang berbeda sehingga akan menambah jumlah data yang telah ada. Juga perlu dilakukan identifikasi adanya senyawa metabolit sekunder


120

yang lain selain Triterpenoid seperti alkaloid, steroid dan flavanoid. Perlu dilakukan penelitian lanjutan untuk memastikan apakah senyawa aktif yang terdapat pada ekstrak bagian tumbuhan pada penelitian ini adalah benar-benar senyawa Triterpenoid yaitu dengan mengisolasi senyawa tersebut.

Daftar Pustaka

[1] Harborne, J.B., Metode Fitokimia Tumbuh-tumbuhan,

(Penterjemah Kosasih Padmawinata dan Iwang Soediro), 1987, terbitan kedua, Penerbit ITB, Bandung.

[2] Robinson, T., Kandungan Organik Tumbuhan Tinggi, (Penterjemah : Prof. Dr. Kosasih Padmawinata), Edisi keenam, Institut Teknologi Bandung, Bandung

[3] Mayer B. N., Ferrigni N. R., Putnam J. E., Jacabsen l. B., Nichols D.E., Mc Laughin J.L., A Convenient General Bioassay for Active Plant Constituent, 1982, Journal of Medicines Plant Research, Vol. 452

[4] DepDikBud, Pengobatan Tradisional Pada Masyarakat Pedesaan Daerah Bengkulu, 1995, Proyek Pengkajian dan Pembinaan Nilai-nilai Budaya Daerah Bengkulu, Bengkulu.

[5] Majang, Yunazar, Isolasi Karakterisasi Senyawa Terpenoid dan Steroid, 2002, Proyek Peningkatan Sumber Daya Manusia, Universitas Andalas, Padang.

[6] Abdi, D., Uji Bioaktivitas Metabolit Sekunder, 2001, Universitas Andalas, Padang

[7] Zulkhaidir, Uji kandungan senyawa steroid dan Bioassay pada beberapa spesies tanaman obat tradisional di desa Kelilik Kabupaten Kepahiang, 2005, skripsi, Universitas Bengkulu, Bengkulu.

[8] Sari, L. N., Uji Pendahuluan Penentuan Adanya Senyawa Triterpenoid dan Bioassay Pada Beberapa Spesies Tanaman Obat Tradisional di Taman Hutan Raya (Tahura) Rajo Lelo Propinsi Bengkulu, 2002, skripsi, Universitas Bengkulu, Bengkulu.

[9] Loisyana, Identifikasi, Inventarisasi Senyawa Triterpenoid dan Uji Brine Shrimp pada Beberapa Tanaman Obat Tradisional Masyarakat Bengkulu, 2005, skripsi, Universiats Bengkulu, Bengkulu.

[10] Yanti, M., Studi Senyawa Terpenoid Pada Beberapa Spesies Tanaman Obat Tradisional di Taman Hutan Raya Rajo Lelo Bengkulu, 2001, Skripsi, Universitas Bengkulu, Bengkulu


121

Tabel 1. Daftar beberapa spesies tanaman obat tradisional masyarakat pedesaan Bengkulu yang mengandung senyawa Triterpenoid

No Nama daearah/Latin/famili Bagian yang mengandung Triterpenoid / kadar

Bagian yang digunakan sebagai obat /khasiatnya

1 Kamboja / Plumeria acuminata /apocynaceae

Kulit batang (+), bunga (+) Getah untuk bengkak

2 Kapuk/randu / Caiba petrandra /bombaceae

Daun (+++), kulit batang (+) Daun untuk obat bengkak

3 Selasih / ocinum bacilicum Akar (+++), daun (+), batang ( +), bunga (+++), akar (+++)

Daun untuk bengkak, campak

4 Mengkudu/ Morinda citrifolis / Daun (++++), kulit batang (+++) , buah (++++)

Daun untuk mencret Buah untuk amandel, cacingan

5 Kunyit/ Curcuma domestica / zingiberaceae Daun ( +) , umbi (+++) Daun untuk bisul, Rimpang untuk kurang darah

6 Kangkung/ Ipomea reptans Akar (++), daun (+), batang (+) Akar untuk wasir 7 Sirsak / Anona muricata / annonaceae Kulit buah,(+), buah (++) Daun untuk asma 8 Sirih / piper betle / piperaceae Daun (+), akar (++) Daun untuk bisul, hidung

berdarah, Batang untuk batuk

9 Serai / Andropogon nardus / gramineae Daun (+), akar (++) Batang untuk batuk 10 Jeruk nipis/ Citrus aurantifolia/ rutaceae Daun ( +), buah (+) Daun untuk batuk 11 Kemuning / Murraya paniculata Daun (++), kulit batang (+++),

buah (+++) Daun untuk batuk, encok

12 Jarak pagar/ Jatropa curcar Daun (+), kulit batang (+) Daun untuk mencret 13 Melati/ Jasminum sambac/ Oleaceae Bunga (++) Daun, bunga untuk campak 14 Daun anjuang/ Cordyline terminalis/ Daun (+++), kulit batang (+) Daun untuk campak

15 Jagung/ Zea mais, Linn/ graminae Kulit batang (++), rambut (+) Rambut untuk darah tinggi 16 Jambu biji/ Psidium guajava, Linn / Daun (+), buah (+) Daun untuk desentri 17 Kecubung/ Datura metel / solanaceae Bunga (++) Daun untuk encok 18 Singkong/ Manihot utilissima, Pohl/ Daun (+++), kulit batang (+) Daun untuk encok 19 Kumis kucing / Orthosiphon stamineus/ Daun (+), batang (++), bunga

(+++) Daun untuk kencing batu

20 Kemiri / Aleurites moluccana, Willd/ Daun (+), kulit batang (+++) Untuk kudis, borok, terkilir 21 Keji beling / Strobilanthes crispus/

acanthaceae Kulit batang (+) Untuk kencing batu

22 Temu lawak / curcuma xanthorrhiza, Roxb / Daun (+), rimpang (+++) Rimpang untuk kencing batu

23 Alang-alang / Imperata cylindrica, beav / gramineae

Daun (+) Akar untuk obat kencing manis

24 Lengkuas / Alpiniaofficinarum hance / zingiberaceae

Daun (+), batang (+) , rimpang (++)

Rimpang untuk panu

25 Alpokat / Parsea americana, Mill/ lauraceae Kulit batang (+) Daun untuk sakit buah pinggang, darah tinggi

26

Dadap / Erythrina orientalis / Daun (++), kulit batang (+) Daun untuk sakit buah pinggang

27 Lamtoro / Leucaena glauca / leguminosae Daun (+), kulit batang (+), bunga (+), buah (++), kulit buah (+)

Untuk cacingan

28 Jahe / Zingiber officinale / zingiberaceae Daun (+), umbi (+++), akar (+) Rimpang untuk masuk angin

29 Jeruk purut / Citrus hystrix, D.C. / rutaceae Daun (+++) Untuk haid tidak teratur 30 Kencur / Kaemphera galanga, L / Daun (+) Untuk kencing manis 31 Pule (pelawi)/ Alstonia scholaris/

apocynaceae Kulit batang (++++), daun (++++), getah (+++)

Kulit batang untuk kencing manis


122

32 Pepaya / Carica papaya / caricaceae Daun (+), batang (+), akar (+), buah (+) bunga (++)

Daun untuk malaria, Getah untuk kulit

33 Jarak (ginje)/ thevetio peruviana / apocynaceae

Daun (+), kulit batang (+++), bunga (+++), kulit buah (+++), biji (++++)

Daun untuk penurun panas

34 Asam jawa / tamarindus indica , L, / Kulit batang (+) Penambah nafsu makan, Keterangan : (+) = sedikit Triterpenoid (++) = cukup banyak Triterpenoid (+++) = mengandung banyak Triterpenoid (++++) = mengandung sangat banyak Sumber : [4], [8], [9] dan [10].

Tabel 2. Hasil Uji Aktivitas Biologis terhadap ekstrak bagian tanaman obat tradisional yang banyak dan sangat banyak mengandung senyawa Triterpenoid

Prosentasi kematian Arthemia salina (%) No Nama Tumbuhan (daerah/latin), bagian yang

diekstrak. 10 ppm 100 ppm 1000 ppm LC -50

1 Anjuang /Cordyline terminalis / daun 0 30 100 442,70 2 Pule / Alstonia Scholaris /daun 3,33 50 100 356,22 3 Serai / Andropogon nardus / daun 100 100 100 <10 4 Jeruk purut / Citrus hystrix, D.C. / daun 50 100 100 < 10 5 Sirih / Piper betle / daun 0 23,33 100 429,74 6 Jambu biji / Psidium guajava, Linn / daun 16,66 26,66 50 988,69 7 Kemiri /Aleurites moluccana, Willd / kulit batang 76,66 96,66 96,66 < 10 8 Kapok randu/ daun 23,33 20 83,33 490,50

LC-50 = Lethal Concentrasion 50 Sumber : [8 ], [9 ], dan [10 ]


Studi Kandungan Etanol Dalam Tapai Hasil Fermentasi Beras Ketan Hitam Dan Putih

Teja Dwi Sutanto, Agus Martono Hp



Abstrak - Etanol merupakan senyawa alkohol yang yang dihasilkan dalam fermentasi bahan makanan yang mengandung karbohidrat termasuk beras ketan. Karena adanya kandungan etanol inilah sehingga bahan makanan hasil fermentasi menjadi lebih enak rasanya. Penelitian ini bertujuan untuk menentukan kandungan etanol yang terdapat dalam tapai hasil fermentasi beras ketan hitam dan putih yang difermentasi dengan menggunakan 3 macam pembungkus yang berbeda yaitu daun, plastik dan kaca. Pada penelitian ini beras ketan difermentasi dengan menggunakan ragi atau yeast selama 7 hari dengan menggunakan 3 macam pembungkus yang berbeda setelah lebih dahulu dimasak dengan cara dikukus, sehingga seluruh etanol hasil fermentasi sudah diubah menjadi asam asetat. Selanjutnya kandungan etanol dalam tapai ditentukan sebagai asam asetat dengan metode titrasi alkalimetri menggunakan larutan standar NaOH 0,1 N. Berdasar hasil pengukuran menunjukkan bahwa kandungan etanol dalam hasil fermentasi ketan putih dengan pembungkus daun, kaca dan plastik berturut-turut adalah 0,0751 %; 0,0599 % dan 0,0338 % sedangkan kandungan etanol dalam hasil fermentasi beras ketan hitam dengan pembungkus daun, kaca dan plastik berturut-turut adalah 0,0407 %; 0,0403 % dan 0,0388 %.

Kata Kunci : Etanol; Fermentasi; Beras Ketan

1. Pendahuluan

Tapai merupakan hasil dari proses fermentasi dari bahan-bahan yang mengandung karbohidrat seperti beras ketan dan ubi kayu. Dalam proses fermentasi yang melibatkan aktifitas mikroorganisme ini terjadi proses pengubahan karbohidrat menjadi etanol, sehingga bahan makanan hasil fermentasi menjadi lebih enak rasanya. Pada proses pembuatan tapai, masyarakat biasanya menggunakan daun pisang sebagai pembungkusnya. Dengan semakin sulitnya mendapatkan daun terutama di daerah perkotaan maka masyarakat beralih ke pembungkus atau wadah alternatif yang lebih mudah diperoleh untuk proses pembuatan tapai seperti plastik, gelas dan wadah yang lain. Namun demikian belum diketahui apakah wadah fermentasi dalam pembuatan tapai ini berpengaruh terhadap kandungan etanol dalam tapai yang dihasilkan. Diantara bahan dasar pembuatan tapai yang sering digunakan adalah beras ketan baik yang hitam maupun yang putih. Tapai dari beras ketan inilah yang paling

banyak dijumpai di toko-toko makanan maupun di Super Market. Namun demikian juga belum diketahui apakah ada perbedaan kandungan etanol dalam tapai yang dihasilkan pada pembuatan dengan kedua macam beras ketan tersebut, sehingga dalam penelitian ini akan dipelajari perbedaan kandungan etanol dalam tapai yang dihasilkan dari bahan dasar beras ketan hitam dan putih dan dengan menggunakan beberapa macam pembungkus yang berbeda yaitu dengan daun, plastik dan gelas.

2. Metode Penelitian Sampel beras ketan hitam dan putih diambil secara acak dari 3 pasar yang ada di kotamadya Bengkulu yaitu Pasar Barukoto, Pasar Minggu dan Pasar Panorama. Tidak ada perlakuan khusus terhadap sampel beras ketan sebelum difermentasi menjadi tapai. Sampel beras ketan dibersihkan kemudian ditimbang sebanyak 500 gram, selanjutnya dimasak dengan cara dikukus menggunakan dandang dengan air sebanyak 300 ml. Setelah masak kemudian didinginkan diatas nampan,

Teja Dwi Sutanto / Jurnal Gradien Vol. 2 No. 1 Januari 2006 :123-125

124

kemudian diberi sedikit larutan gula sebagai starter, selanjutnya dberi serbuk ragi atau yeast sebagai fermentatornya sebanyak 25 gram dan diaduk sampai rata. Langkah selanjutnya dimasukkan kedalam 3 macam wadah yaitu dari daun pisang, plastik dan gelas, ditutup rapat dan dibiarkan terjadi fermentasi selama 7 hari sehingga semua karbohidrat dalam beras ketan yang mula-mula diubah menjadi alkohol sudah berubah semuanya menjadi asam asetat [1]. Setelah fermentasi selama 7 hari selanjutnya kadar etanol dalam tapai yang dihasilkan ditentukan sebagai asam asetat dengan titrasi alkalimetri menggunakan larutan standard NaOH 0,1 N yang sudah distandardisasi lebih dahulu dengan asam oksalat [3].

3. Hasil Dan Pembahasan Variasi Wadah Fermentasi Terhadap Kadar Etanol dalam Tapai Ketan Putih Dari hasil penentuan kadar etanol dalam tapai ketan putih dengan fermentasi menggunakan 3 macam pembungkus yang berbeda ternyata terdapat perbedaan kadar etanol yang cukup signifikan. Hasil penentuan kadar etanol dalam setiap 500 gram beras ketan putih adalah seperti terlihat pada gambar berikut.

0.0751

0.0599

0.0339

0

0.02

0.04

0.06

0.08

Kad

ar E

tano

l (%

)

Daun Gelas Plastik

Wadah Fermentasi

Gambar 1. Variasi Wadah Fermentasi Terhadap Kadar Etanol

dalam Beras Ketan Putih Dari gambar tersebut terlihat bahwa kadar etanol yang paling tinggi diperoleh pada penggunaan daun pisang sebagai wadah fermentasi, kemudian diikuti penggunaan gelas dan yang yang paling kecil adalah pada penggunaan wadah dari plastik, hal ini disebabkan karena pembungkus atau wadah dari daun pisang akan memberikan suasana yang lebih cocok bagi mikrobia fermentator untuk berperan aktif dalam proses fermentasi karbohidrat menjadi etanol. Disamping itu dengan pembungkus daun

pisang yang relatif tidak begitu rapat dibanding pembungkus plastik dan gelas lebih memungkinkan bagi mikrobia azeto bacter yang merupakan mikrobia aerob untuk berperan maksimal dalam proses pengubahan etanol menjadi asam asetat. Variasi Wadah Fermentasi Terhadap Kadar Etanol dalam Tapai Ketan Hitam Hasil penentuan kadar etanol dalam tapai hasil fermentasi 500 gram beras ketan hitam dengan menggunakan tiga macam pembungkus yang berbeda yaitu daun pisang, plastik dan gelas dapat dilihat pada gambar berikut.

0.0407 0.0403 0.0388

0

0.02

0.04

0.06

0.08

Kad

ar E

tano

l (%

)

Daun Gelas Plastik

Wadah Fermentasi Gambar 2. Variasi Wadah Fermentasi Terhadap Kadar

Etanol dalam Beras Ketan Hitam Dari gambar tersebut terlihat bahwa kadar etanol yang paling tinggi sebagaimana pada beras ketan putih, diperoleh pada penggunaan pembungkus daun pisang kemudian diikuti dengan penggunaan gelas dan plastik. Namun demikian perbedaannya tidak terlalu besar sebagaimana pada tapai ketan putih Dari kedua gambar tersebut juga terlihat bahwa dengan wadah yang sama ternyata kadar etanol dalam tapai dari beras ketan putih lebih banyak jika dibandingkan dengan kadar etanol dalam tapai dari beras ketan hitam. Hal ini disebabkan karena struktur amilosa dari beras ketan putih lebih mudah dicerna oleh mikrobia fermentator dibanding ketan hitam, sehingga enzim amilase lebih mudah melakukan pemecahan karbohidrat menjadi etanol dan akibatnya dengan waktu fermentasi yang sama, kadar etanol dalam tapai dari beras ketan putih lebih besar dibanding dari ketan hitam.

4. Kesimpulan Berdasarkan hasil percobaan yang diperoleh dalam penelitian ini maka dapat dibuat kesimpulan sebagai berikut :

Teja Dwi Sutanto / Jurnal Gradien Vol. 2 No. 1 Januari 2006 :123-125

125

1. Kandungan etanol dalam tapai dari beras ketan putih lebih besar dibanding dengan tapai dari beras ketan hitam.

2. Wadah atau pembungkus dari daun pisang lebih baik jika dibandingkan dengan wadah atau pembungkus dari gelas dan plastik untuk fermentasi pembuatan tapai ketan.

3. Kandungan etanol dalam hasil fermentasi beras ketan putih dengan pembungkus daun, kaca dan plastik berturut-turut adalah 0,0751 %, 0,0599 % dan 0,0339 % sedangkan kandungan etanol dalam hasil fermentasi ketan hitam dengan pembungkus daun, kaca dan plastik berturut-turut adalah 0,0407 %, 0,0403 % dan 0,0388 %

Daftar Pustaka

[1] Brewster, R.Q. and Ewen, W.M., Organic Chemistry, 1961,

New Jersey, Prentice Hall Inc. [2] Fessenden, R.J. and Fessenden, J.S., Kimia Organik, Jilid

I, 1989, Jakarta, Erlangga [3] Mudjiran, Analisa Anorganik Kuantitatif, 1982,

Yogyakarta, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Gadjah Mada.

[4] Fox, B.A.ameron, A.G., Food Science, 1987, London, Hodder and Stoughton

[5] Vogel, A.I., Quantitative Inorganic Analysis Theory and Practice, 1958, London, Longmans Green and Co.

[6] Winarno, F.G., Kimia Pangan dan Gizi, 1988, Jakarta, PT. Gramedia


Efek Perendaman Larutan Alkali Terhadap Prilaku Film Kertas Dari Nata de coco Yang Dimodifikasi

Irfan Gustian, Teja Dwi Sutanto, Morina Adfa


Diterima 26 Nopember 2005; disetujui 15 Desember 2005

Abstrak - Telah dibuat film kertas dari nata de coco yang di modifikasi dengan memanfaatkan air kelapa sebagai bahan dasar yang difermentasi. Penelitian ini bertujuan untuk membuat selulosa bakteri dari air kelapa dengan bantuan bakteri Acetobacter xylinum, memodifikasi selulosa bakteri yang dihasilkan menjadi lembaran kertas, menguji sifat fisik serta melihat morphologi kertas yang dihasilkan. Gel nata de coco yang dihasilkan dilakukan dua variasi perlakuan. Variasi yang pertama direndam dalam larutan NaOH 2% pada temperatur 80 – 90 oC selama 2 jam. Variasi yang kedua hasil dari peralakuan pertama direndam lagi dalam larutan NaOH 5% pada temperatur 80 – 90 oC selama 2 jam. Masing-masing variasi di hot press pada tekanan 100 Kgf/cm2 pada suhu 120 oC selama 3 - 5 menit. Hasil pengujian di dapat kerapatan kertas variasi 1 berkisar antara 1,0570 – 1,3330 g/cm2 dan lebih besar dibandingkan dengan kertas variasi 2 yaitu berkisar antara 0,7813 – 1,1666 g/cm2. Sifat penggembungan kertas atau derajat swelling kertas oleh air, variasi 1 memiliki derajat swelling berkisar antara 61,1842 – 91,2162% lebih kecil dibandingkan dengan derajat swelling kertas nata de coco variasi 2 berkisar 77,1428 – 109,4595%. Hasil SEM terlihat permukaan cukup rapat dengan serat yang panjang dengan arah yang tidak beraturan. Kata Kunci: Nata de coco; Acetobacter xylinum; swelling.

1. Pendahuluan

Potensi air kelapa sebagai bahan yang memberikan nilai ekonomis sudah cukup lama diketahui orang. Potensi-potensi tersebut antara lain dapat difermentasi menjadi asam asetat sebagai bahan pembuat cuka, sebagai minuman kesehatan anti diare, minuman olahraga (new sport drink) dan sebagai bahan makanan yaitu nata de coco. Potensi-potensi ini berkaitan dengan kandungan yang ada dalam air kelapa. Kandungan utama air kelapa adalah karbohidrat, termasuk didalamnya gula sederhana, glukosa, fruktosa dan sukrosa maupun gula alkohol utamanya sorbitol (glucitol). Mengingat banyaknya potensi air kelapa sebagai bahan yang bernilai ekonomis, studi pemanfaatannya perlu terus dikembangkan. Salah satu dari beberapa potensi air kelapa yang banyak dikembangkan di Indonesia adalah nata de coco. Nata de coco adalah hasil proses fermentasi air kelapa menggunakan Acetobacter xylinum. Kandungan utama nata de coco adalah selulosa. Selulosa bakterial mempunyai beberapa keunggulan antara lain kemurnian

tinggi, derajat kristalinitas tinggi, mempunyai kerapatan antara 300 dan 900 kg/m3, kekuatan tarik tinggi, dan elastis.. Penelitian yang mengarah pada pengembangan selulosa bakterial sebagai material yang bernilai tambah sudah banyak dilakukan. Beberapa diantaranya adalah penggunaan selulosa bakterial sebagai bahan diafragma tranduser [3], bahan pencampur dalam industri kertas, karakterisasi sifat listrik dan magnetnya, sebagai support untuk sensor glukosa dan sebagai membran dialisis [5]. Pada tulisan ini telah mencoba memanfaatkan nata de coco sebagai material untuk pembuatan film kertas. Telah disadari bahwa selama ini kertas yang digunakan berasal dari cotton linters dan pulp wood hal ini menunjukkan adanya suatu ketergantungan yang sangat tinggi terhadap bahan hayati yang tidak dapat diperbaharui. Untuk mengurangi hal tersebut maka perlu dilakukannya suatu penelitian yang memanfaatkan bahan lain yang banyak terdapat di negara kita. Penelitian-penelitian ke arah mencari bahan baku alternatif lain untuk pembuatan kertas

Irfan Gustian / Jurnal Gradien Vol. 2 No. 1 Januari 2006 :126-129

127

yang berasal dari bahan alam hayati hingga saat ini sudah banyak yang melaporkannya Sedangkan untuk nata de coco masih sedikit peneliti yang melaporkan, antara lain [5], yang mempelajari kinerjanya sebagai membran dialisis dalam sistem larutan encer. Kemudian [4], mencoba mempelajari sifat fisis dan mekanis. Untuk itu maka perlu dilakukannya penelitian untuk mencari sumber bahan alternatif lainnya yang dapat dibuat sebagai kertas.

2. Metode Penelitian Pembuatan nata de coco Satu liter air kelapa yang telah disaring dididihkan. Setelah mendidih ditambahkan 6,7 gram gula pasir, 5 gram Ammonium Sulfat dan 1 mL Asam asetat Glasial. Dalam keadaan panas kemudian dipindahkan ke wadah plastik/baki berukuran 15 x 20 cm dan ditutup dengan kertas. Setelah dingin/ suhu kamar di inokulasi dengan starter bakteri Acetobacter xylinum (10% dari volume media fermentasi) dengan membuka penutup kertas dan ditutup kembali dan difermentasikan/diinkubasi pada suhu kamar selama 3-5 hari. Proses Pemurnian nata de coco Gel nata de coco hasil fermentasi selanjutnya di rendam dengan air mengalir selama 24 jam, kemudian direndam dalam air panas pada suhu 80-90oC selama 2 jam, seteleh itu dicuci dengan NaOH 2% pada suhu 80-90oC selama 2 jam dan direndam lagi dalam air panas pada 80-90oC selama 2 jam. Hasil pemurnian ini di buat 2 variasi, variasi I direndam dengan air mengalir selama 24 jam sebelum dibuat film membran. Variasi II, nata de aren hasil pemurnian di rendam dalam NaOH 5% pada suhu 80-90oC selama 2 jam setelah itu direndam dengan air mengalir selama 24 jam selanjutnya dibuat lembaran kertas. Pengukuran ketebalan gel dan kertas nata de coco Tebal gel nata de coco dan lembaran kertasnya diukur dengan menggunakan thickener pada berbagai bagian gel dan lembaran kertas dan dihitung ketebalan rata-ratanya. Pembuatan kertas untuk kromatografi dari nata de coco

Nata de coco yang telah dimurnikan dengan 2 variasi tersebut dan diukur ketebalannya, kemudian dibuat kertas untuk kromatografi dengan meletekkan nata de coco pada antara dua plat baja yang dilapisi dengan kasa stainles stell dan kain non woven selanjutnya ditekan dengan tekanan 100 kgf/cm2 pada suhu 120oC selama 3-5 menit. kertas yang dihasilkan di simpan dalam dalam plastik untuk proses selanjutnya. Uji Sifat Fisik Uji Kerapatan kertas untuk kromatografi Pengujian kerapatan kertas dilakukan dengan menimbang kertas kering nata de coco , kemudian hasilnya dibagi dengan volume kering. Penentuan volume dilakukan dengan perkalian luas alas x tebal.film. Uji Derajat Swelling kertas Lembaran kertas dari nata de coco dipotong dengan ukuran 2 x 2 cm dan diukur ketebalannya. Lembaran tersebut ditimbang berat awalnya dan direndam dalm aquades selama 7 hari ditimbang setiap selang 24 jam hingga didapat berat kesetimbangan. Setiap kali akan menimbang terlebih dahulu menyeka kedua permukaan film kertas dengan menggunakan tisue.

3. Hasil Dan Pembahasan Pembuatan nata de coco Pembuatan film nata de coco diawali dengan mencampurkan air kelapa dan gula kemudian ditambahkan starter (bakteri Acetobacter xylinum dalam medium cair) setelah melalui pendinginan pada suhu kamar. Setelah masa fermentasi selama 7 hari akan terbentuk gel pada permukaan media cairnya. Gel yang terbentuk ini disebut pellicle. Proses terbentuknya pellicle merupakan rangkaian aktifitas bakteri acetobacter xylinum dengan nutrien yang ada pada media cair. Karena Acetobacter xylinum adalah bakteri yang memproduksi selulosa, maka nutrien yang berperan adalah nutrien yang mengandung glukosa. Dalam penelitian ini nutrien yang mengandung glukosa adalah air


128

kelapa dan gula pasir. Pada gula pasir, glukosa terbentuk melalui reaksi hidrolisis sukrosa dengan air. Reaksi yang terjadi :

OHOH

OH

CH2OHO

O

CH2OH

CH2OHOH

OH

O

OH OHOH

OH

CH2OHO CH2OH

OH

OH

OH

OH2O

CH2OH

Gambar 1. Reaksi hidrolisis sukrosa Glukosa yang berperan dalam pembentukan selulosa adalah glukosa dalam bentuk β sehingga semua glukosa yang ada dalam bentuk α akan diubah dalam bentuk β melalui enzim isomerase yang berada pada bakteri Acetobacter xylinum. Perubahan pada bentuk β terjadi pada gugus OH pada atom C-1 ( ). Reaksinya adalah

OH OHOH

OH

CH 2OHO isomerase

OH

OH

OH

OH

CH 2OHO

Gambar 2. Reaksi Isomerase α-Glukosa menjadi β-Glukosa Tahap berikutnya glukosa berikatan dengan glukosa yang lain melalui ikatan 1,4 β-glikosida. Ikatan ini terjadi antara gugus OH pada atom C-1 ( ) dari satu glukosa β dengan gugus OH pada atom C-4 ( ) dari glukosa β yang lain.

OH

OH

OH

OH

CH2OHO

OH

OH

OH

OH

CH2OHO H2O

OH

O

OH

OH

CH2OHO

OH

OH

OH

CH2OHO

Gambar 3. Reaksi anatara glukosa dengan glukosa Tahap terakhir adalah tahap polimerisasi yaitu pembentukan selulosa . Polimerisasi ini terjadi melalui enzim polimerisasi yang ada pada bakteri acetobacter xylinum. Secara fisik pembentukan selulosa adalah terbentuknya pellicle. Reaksi polimerisasinya:

OH

O

OH

OH

CH2OHO

OH

OH

OH

CH2OHO

Polimerisasi

O

OH

OH

CH2OHO

O

OH

OH

CH2OHO

OH

OH

OH

CH2OHO

Gambar 4. Reaksi Polimerisasi

Gel nata yang dihasilkan ini kemudian di cuci dengan air mengalir selama 24 jam, kemudian dipanaskan dalam air 80 – 90 oC selama 2 jam. Gel nata yang dihasilkan mempunyai ketebalan 3,5 – 4 mm berwarna putih kekuning. Sedangkan tebal kertas yang dihasilkan dari dua variasi perlakuan dan ditekan panas pada tekanan 100 kgf/cm2 dengan suhu 120 oC selama 3 sampai 5 menit berkisar antara 0,03 – 0,035 mm. Sifat Fisik kertas dari nata de coco Uji Kerapatan Pada penentuan kerapatan kertas yang dihasilkan, kertas dipotong dengan ukuran 2x2 cm dan mengalikannya dengan tebal kertas untuk mendapatkan volumenya, kemudian ditimbang beratnya. Tabel 4.1 menunjukan kerapatan kertas nata de coco dari 2 variasi perlakuan, dari hasil pengujian dimana kerapatan kertas yang hanya dicuci dan direndam dengan NaOH 2% (b/v) kerapatannya berkisar antara 1,0570 – 1,3330 g/cm2 dan lebih besar dibandingkan dengan kertas hasil pencucian dengan NaOH 2% (b/v) dan direndam lagi dalam NaOH 5% (b/v) yaitu berkisar antara 0,7813 – 1,1666 g/cm2.

Tabel 1. Kerapatan kertasa Variasi

perendaman Tebal kertas

(mm) Kerapatan

(g/cm2) Variasi I 0,03 1.0570 –1,3330 Variasi II 0,03 0,7813 – 1,1666

Ukuran kertas 2x2 cm Disini terlihat bahwa kertas dari hasil perendaman NaOH 5% (b/v) lebih kecil dibandingkan dengan kertas yang hanya dicuci dan direndam dalam NaOH 2% (b/v). Perbedaan kerapatan ini mungkin disebabkan karena konsentrasi larutan alkali, pada perendaman dengan NaOH 5% (b/v) dapat menyebabkan pengerutan pada gel nata de coco ini terlihat dari bentuk fisiknya sehingga dapat mendegradasi fibril selulosa. Hal ini dapat menyebabkan putusnya hubungan antar fibril selulosa dan mengakibatkan tidak terbentuknya mikrofibril-mikrofibril selulosa, dengan demikian akan mengakibatkan berkurangnya kerapatan membran [1].

fruktosa

α-Glukosa β-Glukosa

Sukrosa glukosa fruktosa


129

Uji Derajat Swelling kertas Sama dengan pengujian kerapatan kertas, pada pengujian swelling kertas nata de coco dipotong dengan ukuran 2 x 2 cm. Berat awal kertas sebelum direndam kedalam air demineral ditentukan terlebih dahulu, sedangkan pertambahan berat kertas setelah direndam ditimbang setiap selang 24 jam selam 7 hari hingga didapat berat kesetimbangan. Tabel 2 menunjukkan derajat swelling kertas oleh air, dimana penyerapan air oleh kertas nata de coco yang hanya dicuci dengan NaOH 2% (b/v) memiliki derajat swelling berkisar antara 61,1842 – 91,2162% lebih kecil dibandingkan dengan derajat swelling kertas nata de coco yang dihasilkan yang direndam lagi dalam NaOH 5% (b/v) berkisar 77,1428 – 109,4595%. Hal ini terjadi disebabkan karena perbedaan perlakuan perendaman kepada kedua kertas. Pada perlakuan perendaman NaOH 5% (b/v) derajat swellingnya lebih besar ini dapat menyebabkan terjadinya pengerutan dan dapat mendegradasi fibril selulosa dan mikrifibril dari selulosa tidak terbentuk hubungan satu sama lainnya dan ketika dilakukannya perendaman dalam air. Kertas yang ini akan menyerap air lebih banyak dimana air akan masuk pada daerah mikrofibril yang tidak berhubungan ini, dibandingkan dengan kertas yang hanya dicuci dan direndam dalam NaOH 2% (b/v). Penyerapan air oleh kertas pada saat kesetimbangan sudah menunjukan derajat swelling yang konstan, ini berarti bahwa swelling air oleh kertas terjadi penngembungan terbatas. Dimana penggembungan ini terjadi, air masuk kedalam daerah-daerah mikrofibril yang susunannya tidak teratur (amorf), penggembungan jenis ini dikenal sebgai penggembungan interfibril [1].

Tabel 2. Derajat swelling kertas nata de coco Variasi perendaman Derajat swelling (%)

Variasi I 61,1842 –91,2162 Variasi II 77,1428 –109,4528 • Ukran kertas 2 x 2 cm • Penimbangan tiap selang 24 jam selama 7 hari

Penggembungan jenis ini disebabkan karena ikatan hidrogen pada daerah amorf lebih lemah dibandingkan dengan daerah kristalin. Derajat swelling ini bila dilihat

dari kerapatannya, kerapatan dari kertas nata de coco yang lebih ringan memungkinkan adanya daerah-daerah mikrofibril yang tidak berhubungan untuk dapat dimasuki oleh air sehingga terjadi swelling yang lebih besar. Hasil Foto SEM dapat dilihat pada Gambar di bawah ini.

.

Gambar. Hasil foto SEM bagian permukaan kertas

4. Kesimpulan Berdasarkan hasil dan pembahasan dapat disimpulkan sebagai berikut: 1. Kertas yang dihasilkan dari nata de coco yang di

modifikasi dengan memanfaatkan aktifitas mikroorganisme dapat dibuat dengan bahan dasar air kelapa.

2. Prilaku fisik yang dihasilkan akibat perendaman larutan alkali menunjukkan adanya perbedaan.

3. Dari foto SEM menunjukkan bahwa arah serat-serat selulosa dari kertas yang dihasilkan tidak satu orientasi/tidak beraturan dan panjang-panjang/ menyatu sehingga dapat mempengaruhi sifat kapileritas.

Daftar Pustaka

[1]. Adisesa, HT, Beberapa perubahan struktur dalam selulosa

pada pengeringan, 1993, Tesis , ITB. [2]. Brandrup, et al., Polymer handbook, ed 4th, 1999, Jhon

Wiley & sons, Inc, New York. [3]. Iguchi, M., et al., Review bacterial cellulose A Masterpiece

of Nature’s arts, 2000, J. Material Sciences, 35 [4]. Indarti. L, dan Yudianti. R, Pengaruh Alkali pada sifat

fisis dan mekanis dari lapisan tipis bioselulosa, 1995, Proseding simposium nasional himpunan polimer Indonesia.

[5]. Shibazaki. H, S.Kuga, F. Onabe and M.Usuda, Bacterial cellulose membrane as separation medium, 1993, Journal of applied Polymer sciences, 50, 965-969.


Efek Kenaikan pH Pada Mekanisme Ekstraksi Cair-Cair Terhadap Asam Asam Karboksilat

Agus Martono HP

Jurusan Kimia Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Bengkulu, Indonesia


Abstraks - Kesadaran manusia akan bersih lingkungan semakin hari kian meningkat, terbukti dengan terbentuknya beberapa organisasi yang bergerak dalam bidang penanggulangan pencemaran lingkungan. Kondisi ini muncul akibat dari kemajuan teknologi industri yang kadang melupakan faktor kesehatan lingkungan. Demikian halnya dengan kemajuan dibidang agroindustri yang sedikit banyak ikut andil dalam mencemari lingkungan khususnya perairan. Hal ini terjadi karena limbah cair dari agroindustri biasanya langsung dibuang ke sungai atau danau. Padahal dalam limbah cair tersebut masih mengandung beberapa asam-asam organik yang apabila dibuang begitu saja jelas akan mencemari lingkungan perairan, namun dengan teknologi ekstraksi cair-cair maka asam-asam organiknya dapat dipisahkan dari limbah cairnya bahkan dapat dimurnikan dan dijual kembali. Sehingga dapat mengimbangi beaya pengolahan limbahnya. Telah dilakukan penelitian dalam ekstraksi terhadap limbah cair agroindustri yang mengandung asam-asam organik dengan sistem ekstraksi cair-cair. Dalam penelitian ini telah ditemukan pengaruh pH inisial larutan terhadap mekanisme ekstrkasi terhadap asam-asam organik yang terdapat dalam limbah cair agroindustri. Adapun asam-asam tersebut adalah asam malat dan asam sitrat dengan konsentarsi sama dengan konsentrasi nyata dalam limbah cair agroindustri yaitu, 5 gram/l dan 10 gram/l larutan. Kisaran pH larutan inisial yang dipergunakan dalam penelitian ini antara 1 sampai 5. Sedangkan pelarut organik yang dipergunakan terdiri dari 70% vol Tributyl phosphate dan 30% vol Dodekan dan perbandingan koefisien antara pelarut organik dengan fasa cairnya mulai dari 1 sampai dengan 5. Hasil dari penelitian ini menunjukkan bahwa pada pH larutan inisial antara 1 – 2 dapat mengekstraksi 86%-90% asam sitart dan 73% - 85% asam malat dari dua macam konsetrasi tersebut diatas. Kata Kunci : Ekstraksi cair-cair; Asam-asam Karboksilat.

1. Pendahuluan Sesuai dengan konsep dasar Produksi Bersih salah satu langkah untuk mengurangi resiko terhadap kesehatan manusia dan pencemaran terhadap lingkungan adalah dengan cara meminimalisasi limbah yang dikeluarkan oleh industri [10]. Usaha untuk meminimalisasi limbah dapat dikelompokkan dalam dua langkah, yaitu mereduksi sumbernya dan pemanfaatan limbahnya. Disamping itu, keuntungan suatu industri antara lain diperoleh dengan adanya efisiensi pada proses produksi, baik dibidang pengolahan maupun pada sektor daur ulang bahan yang dipergunakan pada proses produksinya. Perkembangan teknologi dibidang agroindustri telah mengalami kemajuan yang sangat berarti, namun sisi negatif dari kemajuan tersebut telah menghasilkan limbah industri yang semakin sulit untuk dapat diatasi. Limbah cair agroindustri seperti industri gula tebu

misalnya, merupakan hasil buangan industri yang biasanya langsung dibuang ke sungai. Hal ini akan dapat mengganggu biota yang hidup di sungai, bahkan dapat mengganggu kesehatan manusia, apabila air sungai tersebut dikonsumsi langsung oleh penduduk di sekitarnya. Limbah cair dari industri gula tebu biasanya mengandung asam-asam karboksilat yang larut dalam air, sehingga perlu teknologi khusus untuk memisahkannya [9]. Hingga kini telah dikenal beberapa cara pemisahan asam-asam karboksilat yang terdapat didalam limbah cair dari industri gula tebu, seperti dengan cara pengendapan, destilasi, pengomplekskan, fermentasi dan ekstraksi [8]. Dari beberapa metoda tersebut, ekstraksi merupakan cara yang paling tepat, karena cara ini kecuali lebih hemat energi, sangat mudah dilakukan dan asam-asam dari hasil ekstraksi ini dapat dijual kembali karena memiliki kemurnian yang tinggi [4]. Ada berbagai macam ekstraksi, namun untuk pemisahan asam-asam karboksilat dari limbah cair, dipergunakan ekstraksi cair-cair dengan

Agus Martono HP / Jurnal Gradien Vol. 2 No. 1 Januari 2006 : 130-133

131

campuran pelarut organik baik sebagai ekstraktan maupun diluent-nya [5][6]. Untuk asam-asam organik yang volatile (mudah menguap) seperti asam asetat misalnya, alternatif penyelesaian yang dapat dipergunakan antara lain dengan cara distilasi dan azeotrop atau distilasi ekstraksi menggunakan pelarut dapat juga sistem adsorpsi [6][7]. Namun untuk asam-asam karboksilat yang sukar menguap, cara-cara diatas tidak efektif dikarenakan disamping prosesnya mahal, hasilnya tidak dapat optimal. Untuk pemisahan asam sitrat atau malat dari proses fermentasinya dapat dipergunakan dengan penambahan kalsium hidroksida agar terbentuk garam kalsiumnya, kemudian ditambahkan dengan asam sulfat pada pemurnian asam-asam karboksilatnya. Namun, proses ini akan mengkonsumsi bahan kimia dalam jumlah besar dan menimbulkan pencemaran baru. Ekstraksi dengan pelarut lebih efektif untuk asam-asam karboksilat yang sukar menguap dari limbah cairnya. Pelarut yang dapat dipergunakan antara lain, yang memiliki daya ekstraksi besar dan bersifat basa seperti amina tersier atau phosphine oksida [2].

2. Metode Penelitian

Tributylphosphat (TBP), Dodecane, asam malat, asam sitrat, NaOH, H2SO4 dan indikator dibeli dari Aldrich Chemical Company. Air dimineralisasi diperoleh dari Millipore Mille–Q. Analisis dilakukan dengan menggunakan High Performance Liquid Chromatography (HPLC) model Spectra Seri P 100, integrator Chromjet SP 4400 dan detektor UV spectra-Physics 100 UV-VIS, serta kolom analisa AMINEX HPX-87P pada suhu 35°C.

Berdasarkan hasil analisa limbah cair industri gula tebu [9], ternyata mengandung asam-asam karboksilat seperti; asam malat, asam sitrat dan lain-lain dengan kisaran konsentrasi 0,5 g/l sampai dengan 8 g/l, dalam penelitian ini digunakan dua jenis konsentrasi 5 g/l dari batas bawah dan 10 g/l sebagai batas atas konsentrasi untuk kedua jenis asam tersebut. Sebagai campuran pelarut organik dipergunakan TBP untuk ekstraktan dengan variasi persentasi volume 70 % vol TBP dan 30% vol Dodekan. 50 mililiter fasa organik (ekstraktan + diluent) dan 50 mililiter fasa cair yang berisi larutan asam malat atau sitrat

dimasukkan ke dalam corong pemisah, dikocok selama 3 jam, kemudian diamkan agar terjadi pemisahan antara fasa cair dan fasa organik selama kurang lebih 2 jam. Setelah kedua fasa terpisah sempurna, tempatkan fasa cair ke dalam gelas piala dan fasa organik ke dalam gelas piala yang lain. - Timbang 10 gram fasa cair sebanyak tiga buah, titrasi

masing dengan larutan NaOH 0,01 N, 0,1N atau 1 N sertakan Phenolfptalin (pp) sebagai indikatornya.

- Timbang 10 gram fasa organik sebanyak tiga buah, titrasi dengan larutan NaOH 0,01 N, 0,1 N atau 1 N. Untuk mempermudah jalannya titrasi tambahkan alkohol secukupnya dan sertakan Phenolfptalin (pp) sebagai indikatornya.

- Ulangi hal ini untuk masing-masing komposisi pelarut, terhadap fasa cairnya untuk setiap jenis asam tersebut diatas.

Verifikasi hasil titrasi ini dengan hasil analisa dari HPLC, agar mendapatkan besarnya persentasi asam yang terkestraksi dengan pasti. Analisa dan titrasi sebaiknya dilakukan pada hari yang sama dengan proses ekstraksi agar terhindar dari proses oksidasi dan biodegradasi sampel.

3. Hasil dan Pembahasan Kandungan gugus Phosphoryl yang merupakan ikatan antara Oksigen dan Phosphor dalam TBP, berfungsi sebagai basa Lewis, adalah merupakan ekstraktan yang lebih kuat daripada ekstraktan yang mengandung gugus karbonil, sehingga dapat mengikat asam-asam karboksilat dalam limbah cair. Dalam hal ini hanya asam-asam karboksilat yang tidak terdesosiasi saja yang dapat diikat oleh TBP. Namun karena karakteristiknya secara kimia-fisik maka untuk mempermudah proses ekstraksi diperlukan pelarut pendamping seperti DIPE, MIBK dan Dodecane sebagai diluent [1][11]. Sebagai diluent, DIPE dan MIBK adalah merupakan diluent yang baik dikarenakan kedua pelarut pendamping ini memiliki harga viskositas dan massa jenis rendah serta bersifat polar, sehingga akan dapat meningkatkan transfer asam-asam organik dari fasa cair ke fasa organik, terbukti dengan harga koefisien partisi yang tinggi. Namun karena kedua diluent ini larut dalam air, maka hal ini akan sangat


132

merugikan, karena proses produksinya menjadi tidak ekonomis. Sehingga diluent yang baik untuk proses ini adalah Dodecane, karena diluent ini merupakan hidrokarbon alifatik yang inert dan tidak larut dalam air [2]. Distribusi asam-asam karboksilat dari fasa cair ke fasa organik sangat dipengaruhi oleh kekuatan asamnya, untuk asam-asam poli karboksilat besarnya prosentase asam terekstraksi sangat tergantung dari harga derajad desosiasi pertamanya. Untuk pH inisial larutan antara 2 –2,4 dimana lebih kecil dari pKA asamnya, maka asam ini akan lebih mudah ditarik oleh pelarutnya, sehingga menghasilkan persentase ekstraksi lebih besar [5]. Penambahan NaOH pada fasa cair untuk menaikkan besarnya pH, akan mempengaruhi keberadaan asam-asam organiknya. Karena asam sitrat lebih kuat daripada asam malat, maka asam sitart akan lebih banyak bereaksi dengan NaOH dan mengakibatkan pengurangan jumlah asam sitrat dalam bentuk molekul. Sehingga, pelarut organik akan lebih banyak mengekstraksi asam malat yang tidak banyak mengalami ionisasi. Semakin tinggi harga pH larutan awal, akan menurunkan secara drastis kemampuan ekstraksinya. Fenomena ini juga disebabkan adanya netralisasi asam oleh NaOH tersebut, yang mengakibatkan sistem pengomplekskan antara pelarut organik dengan asam-asam karboksilat menjadi tidak sempurna.

Tabel 1 : Pengaruh Kenaikan pH inisial pada Persentase Ekstraksi Asam Kondisi operasional Asam sitrat : 5 g/l Acide malat : 5 g/l Perbandingan Fasa Pelarut : 1 Campuran Pelarut : (70% vol) T.B.P. +

(30%vol) Dodékan pH

Larutan Awal

% Ekstraksi Asam Sitrat

% Ekstraksi Asam Malat

% Kemurnian Asam Sitrat

1 46,15 25,30 64,80 2 53,47 39,00 57,70 3 34,53 28,27 54,40 4 29,13 31,95 47,62 5 23,52 26,64 46,10

Tabel 2 : Pengaruh Kenaikan pH inisial pada Persentase

Ekstraksi Asam Kondisi operasional Asam sitrat : 10 g/l Acide malat : 10 g/l

Perbandingan Fasa Pelarut : 1 Campuran Pelarut : (70% vol) T.B.P. +

(30% vol) Dodékan pH Larutan

Awal % Ekstraksi Asam Sitrat



1 42,87 35,17 54,60 2 52,93 46,60 50,46 3 26,77 24,52 51,69 4 14,77 15,68 52,21 5 13,15 14,07 47,25


Ekstraksi Asam Kondisi operasional Asam sitrat : 5 g/l Acide malat : 5 g/l Perbandingan Fasa Pelarut : 3 Campuran Pelarut : (70% vol) T.B.P. +

(30% vol) Dodékan pH

Larutan Awal




1 74,20 60,70 55,24 2 78,43 67,30 53,70 3 47,17 46,40 49,82 4 28,63 32,00 47,15 5 22,62 23,97 47,75



(30% vol) Dodékan pH

Larutan Awal




1 73,53 55,51 56,98 2 77,02 70,84 52,09 3 39,76 39,61 50,10 4 19,45 23,24 45,56 5 10,86 11,63 48,29



(30% vol) Dodékan


133

pH Larutan Awal




1 85,82 72,94 54,30 2 89,97 77,50 52,17 3 50,58 52,25 48,59 4 28,28 34,24 45,16 5 15,08 18,76 43,77



(30% vol) Dodékan pH Larutan

Awal % Ekstraksi Asam Sitrat



1 83,00 72,60 53,34 2 85,20 79,20 51,82 3 41,92 46,12 47,61 4 10,20 14,75 40,88 5 2,40 3,50 40,68

Dari hasil analisa (Tabel 1 – 6), nampak bahwa semakin tinggi pH larutan awal akan menurunkan harga/ persantase ekstraksinya. Untuk berbagai jenis perbandingan komposisi pelarut (ekstraktan + diluent) dengan fasa cairnya terhadap berbagai jenis asam yang diekstrak, akan dapat menaikkan persentase ekstraksi untuk kedua jenis asamnya. Hal ini logis, dikarenakan semakin besar konsentrasi ekstrkatan akan menaikkan daya ekstraksinya.

Daftar Pustaka [1] Achour J., Malmary G., Putranto A., Nezzal G.,

Mourgues J,F., Molinier J., Liquid-liquid Equilibrium of Lactic Acid between Water and Tris (6-methylheptyl) amine and Trybutylphosphate in Various Diluents., 1994, J. Chem. Eng. Data., 39., 711 – 713.,

[2] Baniel, A.M; Blumberg, R; Hajdu, K, Recovery of Acids from Aqueous Solutions, 1981, U.S. Patent, 4, 275, 234.

[3] Busche, R.M., The Business of Biomass., 1985, Biotechnol.Prog, 1, 165 – 180.

[4] Juang, R.S. ; Huang, R.H., -, 1997, Chem Eng J., 65, 47– 53.

[5] Kertes A.S., King C.J., - , 1986, Biotechnol. Bioeng., 28, 277.

[6] King C.J., Advances in Separation techniques : Recovery of Polar Organic from aqueous solution, 1993, Proceedings of 11th International Congress of Chemical Engineering. Chisa, Praha.

[7] Kuo, Y., Munson, C.L., Rixey, W.G., Garcia, A.A., Frierman, M., Use of Adsorbents for Recovery of Acetic Acid from aqueous Solutions., I – Factors Governing Capacity., 1987, Sep. Purif. Methods, 16, 31 –64.

[8] Lockwood, L.B., Production of Organic Acids by Fermentation. In Microbial Technology; Peppler, 1979, H.J., Perlman, D., Eds Academics, New York, pp 356 –387.

[9] Myrtil – Celestine D., Parfait A., 1988, Int. Sugar J., 90, 28–32.

[10] Pujiastuti, L., Produksi Bersih., 1999, DirJen DIKTI. [11] Putranto A., Malmary G., Albet J., Molinier J., Separation

of Organic Acids from Simulated Aqueous Effluents of Agroindustries through Solvent Extraction, 2000, Proceedings of 5th Indonesian Studetnt’s Scientific Meeting : 224 – 227.

[12] Yang S,T., White S, A., Tsiung Hsu S., -, 1991, Ind. Eng. Chem. Res., 30, 1335–1342.


Karakterisasi Matrik Leslie Ordo Tiga

Mudin Simanihuruk, Hartanto

Jurusan Matematika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Bengkulu, Indonesia


Abstrak - Misalkan L = [ lij ] adalah matrix bujursangkar berorder n sedemikian hingga l1j ≥ 0 untuk j = 1,2,...,n ; 0 < li(i-1) ≤ 1 untuk i = 2,3,...,n dan lij = 0 untuk i = 2,3,...,n ; dan j = 1,2,...,n tetapi j ≠ i - 1. Matrix L disebut matrix Leslie. Perhatikanlah bahwa setiap pemangkatan dari matrix Leslie L = 0 0 6

0 00 0

12

13

adalah I, L atau L2. Pada makalah ini,

akan diberikan karakterisasi matrix Leslie L berorder tiga sedemikian hingga setiap pemangkatan dari L adalah I, L atau L2 . Kata Kunci : matrix Leslie ; demografi.

1. Pendahuluan Salah satu model pertumbuhan populasi yang sering digunakan para ahli demografi adalah Model Leslie. Para perempuan didistribusikan kedalam kelompok berdasarkan usia. Apabila banyaknya kelompok sama

dengan n dan ( )kiX adalah banyaknya perempuan

pada Kelompok i pada pengamatan tk maka dapat ditunjukkan bahwa

1( )kX = a1

11( )kX − + a2 2

1( )kX − +

....+ an nkX ( )−1 , ai ≥ 0, i = 1,2,...,n dan

ikX +1

( ) =

bi ikX ( )−1 , 0 < bi ≤ 1, i = 1,2,..,n-1. Dua persamaan

yang terakhir dapat ditulis dalam bentuk matrik ( )kX = L ( )kX −1 , k = 1,2,... di mana

( )kX = 1

2

3

( )

( )

( )

( )

.

.

.

k

k

k

nk

XXX

X

dan L = a a a a ab

b

b

n n

n

1 2 3 1

1

2

1

0 0 0 00 0 0 0

0 0

0 0 0 0

. . .

. . .

. . .. . . . . .. . . . . . . .. . . . . . . .

. . .

−

−

.

Matrik L disebut matrik Leslie. Dari persamaan yang terakhir dapat ditunjukkan bahwa ( )kX = Lk (0)X .

Apabila Lk = I, maka ( )kX = (0)X . Jadi pada saat tk

jumlah populasi kembali sama dengan jumlah populasi awal. Pada makalah ini akan diberikan syarat perlu dan cukup bagi matrik Leslie L ordo tiga agar Lk

= I. Apabila Lk = I, maka setiap pemangkatan Lr akan menghasilkan L atau L2,..., atau Lk-1.

2. Pertumbuhan populasi model Leslie Misalkan L adalah usia maksimum yang dapat dicapai oleh perempuan pada suatu populasi. Apabila para perempuan itu dibagi bagi kedalam n kelompok berdasarkan usia, maka jarak interval masing-masing kelompok adalah L/n. Dengan demikian Kelompok 1, adalah mereka yang berusia [0,L/2), Kelompok 2 adalah mereka yang berusia [L/n,2L/n),..., Kelompok n adalah mereka yang berusia [(n-1)/L,L]. Setiap saat, komposisi jumlah perempuan dalam kelompok dipengaruhi oleh tiga faktor yaitu faktor kelahiran, kematian dan pertambahan usia. Misalkan ti adalah waktu pengamatan, i = 0,1,..., k,.. . Pada Model Leslie jarak pengamatan ti-1 ke ti sama

Mudin Simanihuruk / Jurnal Gradien Vol. 2 No. 1 Januari 2006: 134-138 135

dengan jarak interval kelompok. Jadi ti = iL/n, untuk i = 0,1,2,....

Misalkan ai adalah rata-rata banyaknya anak perempuan yang lahir dari setiap Kelompok i dan bi adalah perbandingan antara banyaknya perempuan yang bertahan hidup sehingga mampu masuk kedalam Kelompok i+1, dengan banyaknya perempuan dalam Kelompok i. Perhatikanlah bahwa ai ≥ 0, i=1,2,3,...,n dan 0 < bi ≤ 1. Perhatikan juga paling sedikit satu ai > 0, karena kalau tidak berarti proses kelahiran tidak terjadi dan bi > 0, karena kalau tidak, maka tidak ada perempuan yang bertahan masuk kedalam Kelompok i+1. Misalkan

ikX ( ) adalah banyaknya perempuan pada

Kelompok i pada pengamatan tk untuk i = 1,2,...,n. Pada saat pengamatan tk banyaknya perempuan pada Kelompok 1 sama dengan banyaknya anak perempuan yang lahir di Kelompok 1 dari waktu tk-1 ke tk + banyaknya anak perempuan yang lahir di Kelompok 2 dari waktu tk-1 ke tk + ... + banyaknya anak perempuan yang lahir di Kelompok n dari waktu tk-1 ke tk . Jadi

1( )kX = a1

11( )kX − + a2

21( )kX − + ....+ an

nkX ( )−1 (1)

Karena jarak interval setiap kelompok sama dengan jarak dua pengamatan yang berurutan, maka semua perempuan yang berada pada Kelompok i+1 pada saat pengamatan tk+1 berada pada Kelompok i pada saat pengamatan tk. Oleh karena itu banyaknya perempuan pada Kelompok i+1, i = 1,2,...,n-1, pada saat pengamatan tk sama dengan banyaknya perempuan yang masih hidup pada Kelompok i pada waktu tk-1 ke tk . Jadi

ikX +1

( ) = biikX ( )−1 , i = 1,2,..,n-1 (2)

Persamaan (1) dan (2) dapat ditulis dalam bentuk matrik

( )kX = L ( )kX −1 , k = 1,2,... (3)

di mana ( )kX = 1

2

3

( )

( )

( )

( )

.

.

.

k

k

k

nk

XXX

X

dan

L = a a a a ab

b

b

n n

n

1 2 3 1

1

2

1

0 0 0 00 0 0 0

0 0

0 0 0 0

. . .

. . .

. . .. . . . . .. . . . . . . .. . . . . . . .

. . .

−

−

.

Matrik L disebut matrik Leslie. Dari persamaan rekursif (3) dapat ditunjukkan bahwa

( )kX = Lk(0)X . Apabila Lk = I, maka ( )kX = (0)X .

Berarti pada saat tk jumlah populasi kembali sama denga jumlah populasi awal.

Matriks L = 0 0 60 0

0 0

12

13

adalah matrik Leslie yang

mempunyai sifat-sifat yang menarik. Bernadelli (1941) dalam [2] menunjukkan bahwa setiap pemangkatan dari L akan menghasilkan matrik identitas I, L atau L2. Dapat ditunjukkan bahwa nilai eigen dari matrik Leslie

L di atas adalah 1, − +12

12

3i , dan − −12

12

3i

dan L tidak mempunyai nilai eigen yang dominan di mana nilai eigen dominan adalah nilai eigen positip λ 1 sedemikian hingga λ 1

≥ λ k . Apakah setiap

matrik Leslie L ordo tiga yang nilai eigennya tidak dominan merupakan syarat cukup agar Lt = I untuk suatu bilangan bulat t ? Ternyata hal ini tidak benar

karena matrik Leslie L = 0 0 480 0

0 0

12

13

mempunyai nilai

eigen 2, -1+ i 3 , -1 - i 3 dan di antaranya tidak ada nilai eigen yang dominan, akan tetapi dengan mudah dapat ditunjukkan bahwa L3 ≠ I. Jadi nilai eigen yang dominan dari matrik Leslie L tidak dapat kita harapkan memberikan karakteristik bagi matrik Lt = I untuk suatu bilangan bulat t. Pada makalah ini kita tunjukkan syarat perlu dan cukup agar matrik Leslie L memiliki sifat Lt = I untuk suatu bilangan bulat positip t ≥ 2. Karakteritik lain dari matrik Leslie dapat dibaca pada Linear Algebra [1]. Notasi dan definisi lain yang digunakan dalam makalah ini mengikuti [1]. Yang dimaksud dengan


matrik diagonal adalah matrik bujursangkar yang entrinya diluar diagonal utama sama dengan nol. Matrik diagonal D ordo n dinyatakan dengan D = diag(d1, d2,...,dn). Apabila setiap di > 0 maka matrik diagonal D disebut matrik diagonal positip. Selanjutnya notasi (A)ij menunjukkan entri dari matrik A pada baris ke i dan kolom ke j.

3. Syarat perlu dan cukup Dua Lemma berikut akan menunjukkan bahwa entri dari matrik Leslie L ordo tiga adalah bilangan positip. Lemma 1: Misalkan A adalah matrik ordo tiga. Jika semua entri dari A pada kolom ketiga sama dengan nol, maka semua entri dari At pada kolom ketiga sama dengan nol.

Bukti: Misalkan At-1 = b b bb b bb b b

11 12 13

21 22 23

31 32 33

, t ≥ 2. Akan

ditunjukkan bahwa semua entri dari At pada kolom ketiga semuanya sama dengan nol. Perhatikanlah bahwa At = At-1A dan (At)13 = b11 a13 + b12a23 + b13a33, (At)23 = b21a13 + b22a23 + b23a33, (At)33 = b31 a13 + b32a23 + b33a33 . Karena a13 = a23 =a33 = 0 maka (At)13 = (At)23 = (At)33 = 0. Dengan demikian Lemma terbukti.

Lemma 2: Misalkan L = a b cx

y0 0

0 0

, a ≥ 0, b ≥ 0, c ≥

0, 0 < x ≤ 1 dan 0 < y ≤ 1. Jika Lt, t ≥ 2, adalah matrik diagonal positip maka c>0. Bukti: Misalkan c = 0. Aplikasikan Lemma 1 terhadap Lt diperoleh (Lt)13 = (Lt)23 = (Lt)33 = 0. Karena (Lt)33 = 0, maka Lt bukan matrik diagonal positip, kontradiksi dengan hipothesis dari Lemma. Oleh karena itu c > 0.


y0 0

0 0

, a ≥ 0, b ≥ 0, c > 0,

0 < x ≤ 1 dan 0 < y ≤ 1. Jika a > 0 atau b > 0 maka ada dua entri positip pada baris pertama dari Lt, t ≥ 2. Bukti: Lemma akan dibuktikan dengan induksi pada t.

Basis Induksi: Perhatikanlah bahwa baris pertama dari L2 yang dinyatakan dengan 2

1B sama dengan

21B = [a b c] a b c

xy0 0

0 0

= [a2 + bx ab + cy ac].

Perhatikanlah bahwa cy > 0. Jika a > 0 maka ac > 0. Jadi ada dua entri positip dari L2 pada baris pertama yaitu cy + ab dan ac. Demikian pula apabila b > 0, maka ada dua entri positip dari L2 pada baris pertama yaitu cy + ab > 0 dan a2 + bx > 0. Hipothesa Induksi: Misalkan ada dua entri positip pada baris pertama dari Lt-1. Langkah Induksi: Akan ditunjukkan ada dua entri positip pada baris pertama dari Lt.

Misalkan 11−tB = [p q r], adalah baris pertama Lt-1.

Berdasarkan induksi hipothesa diperoleh: (i). p > 0 dan q> 0, (ii) p > 0 dan r > 0, atasu (iii) q > 0 dan r > 0.

Perhatikanlah bahwa tB1 dari Lt adalah

tB1 = [p q r] a b cx

y0 0

0 0

= [ap + qx pb + yr pc].

Jika p > 0 dan q > 0 maka ap + qx > 0 dan pc > 0, jadi ada dua entri positip pada baris pertama dari Lt. Jika p > 0 dan r > 0, maka pb + yr > 0 dan pc > 0, jadi ada dua entri positip pada baris pertama Lt. Jika q > 0 dan r > 0, maka ap + qx > 0 dan pb +yr > 0, jadi ada dua entri positip pada baris pertama Lt. Dengan demikian Lemma terbukti.


y0 0

0 0

adalah matrik

Leslie di mana a ≥ 0, b ≥ 0, c ≥ 0, 0 < x ≤ 1 dan 0 < y ≤ 1. Jika Lt adalah matrik diagonal positip maka a = 0 dan b = 0. Bukti: Karena Lt adalah matrik diagonal positip maka Lt memenuhi kondisi Lemma 2. Oleh karena itu kita simpulkan c > 0. Selanjutnya misalkan a > 0 atau b >0. Kita akan menunjukkan bahwa pemisalan ini akan


menimbulkan kontradiksi seperti berikut. Karena a > 0 atau b > 0, maka dengan menggunakan Lemma 3 kita peroleh bahwa Lt mempunyai dua entri positip pada baris pertama. Berarti Lt bukan matrik diagonal yang dengan sendirinya bukan matrik diagonal positip, kontradiksi dengan hipothesis dari Lemma. Oleh karena itu pemisalan a > 0 atau b > 0 adalah salah. Jadi a = 0 dan b = 0. Dengan demikian Lemma terbukti. Selanjutnya akan kita buktikan bahwa syarat perlu dan cukup agar matrik Leslie L ordo tiga memiliki sifat L3k = I, untuk bilangan bulat k ≥1. Teorema berikut menunjukkan bahwa L3k = diag((cxy)k, (cxy)k, (cxy)k) bila L adalah matrik ordo tiga di mana entri bada baris pertama kolom pertama dan kolom kedua sama dengan nol.

Teorema 1: Misalkan L = a b cx

y0 0

0 0

, c > 0, 0 < x ≤ 1

dan 0 < y ≤ 1. Jika a = 0 dan b = 0 maka L3k = diag((cxy)k, (cxy)k, (cxy)k) untuk k ≥ 1. Bukti: Lemma akan dibuktikan dengan induksi pada k. Basis Induksi: Dengan mudah dapat ditunjukkan bahwa

L2 = 0 00 0

0 0

cycx

xy

dan L3 = L2.L =

diag(cxy,cxy,cxy). Jadi Teorma di atas benar untuk k= 1. Hipothesa Induksi: Misalkan teorema di atas benar untuk k-1 yaitu L3(k-1) = diag ((cxy)k-1, (cxy)k-1, (cxy)k-

1). Langkah induksi : Akan ditunjukkan bahwa L3k = diag ((cxy)k, (cxy)k, (cxy)k). Perhatikanlah bahwa L3(k-1)+1 = L3(k-1).

L= ( )

( )( )

cxycxy

cxy

k

k

k

−

−

−

1

1

1

0 00 00 0

.0 0

0 00 0

cx

y

= 0 0

0 00 0

1

1

1

( )( )

( )

cxy ccxy x

cxy y

k

k

k

−

−

−

,

L3(k-1) + 2 = L3(k-1)+1 .

L = 0 0

0 00 0

1

1

1

( )( )

( )

cxy ccxy x

cxy y

k

k

k

−

−

−

.0 0

0 00 0

cx

y

= 0 00 0

0

1

1

1 1

( )( )

( ) ( )

cxy cycxy cx

cxy xy cxy y

k

k

k k

−

−

− −

dan L3(k-1) + 3 = L3(k-1)+2 .

L = 0 00 0

0

1

1

1 1

( )( )

( ) ( )

cxy cycxy cx

cxy xy cxy y

k

k

k k

−

−

− −

.0 0

0 00 0

cx

y

= diag((cxy)k-1, (cxy)k-1, (cxy)k-1). Dengan demikian teorema di atas terbukti. Selanjutnya teorema berikut menjelaskan syarat perlu dan cukup bagi matrik Lesli L agar L3k adalah matrik diagonal positip.


y0 0

0 0

adalah matrik

Leslie di mana a ≥ 0, b ≥ 0, c ≥ 0, 0 < x ≤ 1 dan 0 < y ≤ 1. Matrik L3k adalah matrik diagonal positip jika dan hanya jika a = b = 0 dan c > 0. Bukti: Misalkan L3k adalah matrik diagonal positip. Dengan menggunakan Lemma 2 terhadap L3k diperoleh c > 0. Selanjutnya aplikasikan Lemma 4 terhadap L3k diperoleh a = b = 0. Sebaliknya misalkan a = b = 0 dan c > 0. Dengan menggunakan Teorema 1 kita peroleh L3k = diag((cxy)k, (cxy)k, (cxy)k). Karena c > 0 maka cxy > 0. Oleh karena itu L3k adalah matrik diagonal positip. dengan demikian teorema terbukti. Akhirnya kita sampai pada karakterisasi dari matrik L sehingga L3k = I sebagaimana dinyatakan oleh teorema berikut.


y0 0

0 0

adalah matrik

Leslie di mana a ≥ 0, b ≥ 0, c ≥ 0, 0 < x ≤ 1 dan 0 < y


≤ 1. Matrik L3k = I jika dan hanya a = b = 0 dan cxy = 1. Bukti: Misalkan L3k = I. Karena semua entri pada diagonal utama sama dengan 1, maka jelas L3k adalah matrik diagonal positip. Berdasarkan Lemma 4, kita peroleh a = b = 0. Selanjutnya berdasrkan Teorema 1 kita peroleh L3k = diag((cxy)k, (cxy)k, (cxy)k). Karena L3k = I, maka (cxy)k = 1. Akibatnya cxy = 1. Sebaliknya misalkan a = b = 0 dan cxy = 1. Karena cxy = 1 maka jelas c > 0. Perhatikanlah kondisi Teorema 2 yaitu a = b = 0 dan c > 0 terpenuhi. Oleh karena itu berdasarkan Teorema 2 kita peroleh bahwa L3k adalah matrik diagonal positip. Dengan menggunakan Teorema 1 kita peroleh L3k = diag((cxy)k, (cxy)k, (cxy)k). Karena cxy = 1 maka L3k = I. Dengan demikian teorema terbukti. Dari Teorema 3 kita lihat apabila a = b = 0 dan cxy = 1 maka L3k = I, L3k+ 1 = L dan L3k+2 = L2 dan L3k+3 = I untuk k ≥ 1.

4. Simpulan dan saran Penelitian ini telah berhasil mengkarakterisasi matriksLeslie L berordo tiga sebagaimana dinyatakan dalam Teorema 3. Karakterisasi tersebut adalah

sebagai berikut: Misalkan L = a b cx

y0 0

0 0

adalah

matrik Leslie di mana a ≥ 0, b ≥ 0, c ≥ 0, 0 < x ≤ 1 dan 0 < y ≤ 1. Matrik L3k = I jika dan hanya jika a = b = 0 dan cxy = 1. Tekhnik-tekhnik pembuktian yang digunakan pada penelitian ini belum dapat digeneralisasi untuk mengkarakterisasi matrik Leslie berordo n, untuk nilai n yang cukup besar. Karakterisasi matrik Leslie L berordo tiga di atas menunjukkan suatu dugaan (conjecture) berikut: Apabila a1, a2,..., an adalah entri pada baris pertama dari matrik Leslie L berordo n dan b1, b2,..., bn-1 di mana bi adalah entri pada baris ke i+1 dan pada kolom ke i, maka L3k = I jika dan hanya jika a1 = a2 = ...= an-1 = 0, b1b2...bn-1an = 1. Salah satu kesukaran yang ditemui untuk membuktikan Conjecture ini terletak pada kesulitan untuk mengeneralisasi Teorema 1 untuk

sebarang matrik Leslie L berordo n. Kami sarankan agar peneliti lain meneliti lebih lanjut tentang kebenaran conjecture di atas.

Daftar Pustaka [1] Anton,H. and Rorres,C., Elementary Linear Algebra

(7th ed.), 1994, John Wiley & Sons,Inc. New York. [2] Tarumingkeng,R., Dinamika Populasi, 1994, Pustaka

Sinar Harapan, Jakarta.


Bentuk Umum Perluasan Teorema Pythagoras

Mulia Astuti, Buyung Keraman, Ulfasari Rafflesia


Diterima 1 September 2005; disetujui 20 Desember 2005

Abstrak - Penelitian ini membahas perluasan teorema Pythagoras melalui pendekatan hubungan kesetaraan pada luas daerah. Secara matematis luas daerah diukur berdasarkan variabel-variabel yang terkait dan kesetaraan pada luas daerah dikembalikan pada kesetaraan fungsi-fungsi yang mengacu pada variabel tersebut. Perluasan teorema Pythagoras di 2R dengan pendekatan hubungan kesetaraan pada luas daerah dapat menjelaskan hubungan luas daerah yang berkaitan dengan panjang sisi-sisi segitiga siku-siku. Fokus dari penelitian ini adalah membahas perluasan teorema Pythagoras di nR . Kata Kunci : Hubungan Kesetaraan; Teorema Pythagoras.

1. Pendahuluan Ausri [3] telah menyelidiki perluasan teorema Pythagoras

di 2R dengan menggunakan hubungan kesetaraan pada luas daerah yang berkaitan dengan sisi segitiga siku-siku. Hasil yang diperoleh adalah penjumlahan luas daerah yang berkaitan dengan sisi siku-siku sama dengan luas daerah yang berkaitan dengan sisi miringnya. Makalah ini akan membahas bentuk umum perluasan

teorema Pythagoras yaitu di nR , dengan menggunakan hubungan kesetaraan pada luas daerah, pada persamaan dengan bentuk umum :

),...,,(... 2122

32

22

1 nn xxxgxxxx ++++= (1)

atau pada persamaan dengan bentuk umum 22

22

121 ...),...,,( nn xxxxxxf +++= (2)

Dalam hal ini ),...,,( n21 aaa , ia > 0, ni1 ≤≤

merupakan solusi persamaan (1) atau persamaan (2).

2. Hubungan Kesetaraan [4] Dalam mendefinisikan hubungan kesetaraan, diperlukan notasi-notasi berikut: 1. Notasi ( ) ( ) ( ) ( ) ( ) ( )x,x,xG,xF,xg,xf aaaaaa γφ

menyatakan fungsi kontinu untuk suatu konstanta R∈> aa ,0 dan [ ]a,0x∈

2. Notasi ( ) ( ) ( ) ( ) ( ) ( )x,x,xG,xF,xg,xf aaaaaa −−−−−− γφ

menyatakan fungsi kontinu untuk suatu konstanta R∈> aa ,0 dan [ ]0,ax −∈

3. Notasi ( ) ( )( )xg,xfR aa menyatakan daerah di bidang

kartesian yang dibatasi oleh: (i) Fungsi ( )xfa dan fungsi ( )xga dengan

( )xga ≤ ( )xfa [ ]a,0x∈∀

(ii) Garis vertikal 0x = dan ax = . 4. Notasi ( ) ( )( )xg,xfA aa menyatakan luas daerah

( ) ( )( )xg,xfR aa .

Definisi 2.1 (Hubungan kesetaraan pada selang) [3] Dua selang dikatakan setara jika kedua selang tersebut mempunyai panjang yang sama. Misalnya : 1. Selang 1xt0 ≤≤ setara dengan selang

21

2 xtxx

0 ≤≤ , untuk 0x1 ≠

2. Selang 40 ≤≤ t setara dengan selang 106 ≤≤ t . Definisi 2.2 (Hubungan kesetaraan pada fungsi) [3] (i) Fungsi ( )xfa dikatakan setara dengan ( )xFb ,

ditulis ( ) ( )xFxf ba π , jika:

( ) [ ]bxxbaf

abxF ab ,0∈∀

=

Mulia Astuti / Jurnal Gradien Vol. 2 No. 1 Januari 2006 : 139-143

140

(ii) Fungsi ( )xf a− dikatakan setara dengan

( )xF b− ditulis ( ) ( )xFxf ba −− π , jika:

( ) [ ]0,bxxbaf

abxF ab −∈∀

= −− .

Misalnya :

Fungsi ( ) axxaxfa ≤≤−= 022 setara dengan fungsi

( ) bxxbxFb ≤≤−= 022 .

Definisi 2.3 (Hubungan kesetaraan pada luas daerah) [3] (i) Daerah ( ) ( )( )xgxfR aa , setara dengan daerah

( ) ( )( )xGxFR bb , , ditulis

( ) ( )( ) ( ) ( )( )xGxFRxgxfR bbaa ,, π jika

( ) ( )xFxf ba π dan ( ) ( )xGxg ba π

(ii) Luas daerah ( ) ( )( )xgxfR aa , namakan

( )( ))(, xgxfA aa dihitung dengan:

( )( ) ( )( )dxxgxfxgxfAa

aaaa ∫ −=0

)()(, .

Definisi 2.4 [2] Misalkan n adalah bilangan bulat positif, n-tupel terurut didefinisikan sebagai urutan n bilangan riil

)a,,a,(a n21 … . Himpunan semua n-tupel terurut

dinamakan Ruang-n Euclidis dan dinotasikan dengan nR

3. Contoh Perluasan Teorema Pythagoras di 2R

Perluasan teorema Pythagoras di 2R menggunakan hubungan kesetaraan pada luas daerah yang berkaitan dengan sisi-sisi segitiga siku-siku. Berikut contoh

perluasan teorema Pythagoras di 2R . 1. Pada masing-masing sisi segitiga siku-siku terdapat

daerah beraturan yang berbentuk bujur sangkar. Hubungan luas bujur sangkar tersebut adalah a2 + b2 = c2.

2. Pada masing-masing sisi segitiga siku-siku terdapat daerah yang beraturan berbentuk setengah lingkaran dengan diameter sama dengan sisi segitiga siku-siku. Hubungan antara luas setengah lingkaran

tersebut adalah: 2228π

8π

8π cba =+

3. Pada masing-masing sisi segitiga siku-siku terdapat daerah yang beraturan berbentuk segitiga yang setara.

Misalkan segitiga itu membentuk sudut-sudut α,β dan θ. Hubungan luas ketiga segitiga setara di atas adalah:

αsin2βsinθsin2

αβθa +

αsin2βsinθsin2

αβθb =

αsin2βsinθsin2

αβθc

Jika pada masing-masing sisi segitiga siku-siku terdapat daerah yang tidak beraturan, maka teorema 3.1 memperlihatkan hubungan antara luas daerah yang terdapat pada masing-masing sisi segitiga siku-siku.

Teorema 3.1 [4]

Misalkan a dan b adalah panjang sisi siku-siku dan c adalah sisi miring segitiga siku-siku. Misalkan daerah

( ) ( )( ),, xgxfR aa ( ) ( )( )xGxFR bb , dan ( ) ( )( )xxR cc γ,φ

setara maka: ( ) ( )( ) ( ) ( )( ) ( ) ( )( )xxAxGxFAxgxfA ccbbaa γ,φ,, =+

Akibat 3.2 [4] Misalkan a ,b dan c masing-masing adalah panjang sisi segitiga sebarang, dengan sisi c berlawanan dengan sudut θ . Misalkan daerah ( ) ( )( )xgxfR aa , , ( ) ( )( )xGxFR bb ,

dan ( ) ( )( )xxR cc γ,φ setara. Maka:

( ) ( )( ) ( ) ( )( )

( ) ( )( ) ( ) ( )( )xxAc

abxGxFA

xgxfAxxA

ccbb

aacc

γ,φθcos2

,

,γ,φ

2−

+=

4. Bentuk Umum Perluasan Teorema Pythagoras

Bentuk umum perluasan teorema Pythagoras, yaitu di ruang-n Euclidis, dilakukan dengan memperluas pasangan

terurut 0,0),( >> baba di 2R yang merupakan solusi

persamaan lingkaran 222),( cyxyxf =+= , menjadi

n-tupel terurut niaaaa in ≤≤> 1,0,),...,,( 21 di nR

yang merupakan solusi persamaan (1) atau merupakan solusi persamaan (2). Sehingga bentuk umum perluasan teorema Pythagoras ini, diperoleh dengan memperluas Teorema 3.1 menjadi teorema 4.2 dan Akibat 3.2 menjadi teorema 4.1 yaitu sebagai berikut :


141

Teorema 4.1 Misalkan niaaaa in ≤≤> 1,0,),...,,( 21 adalah solusi

Persamaan (1). Misalkan daerah ),(),...,,(),,( 2211 nanaaaaa gfRgfRgfR setara,

maka berlaku:

),(),...,,(

),(...),(),(

1121

21

2211

aan

nanaaaaa

gfAa

aaag

gfAgfAgfA

+

++=

Bukti: Berdasarkan Definisi 2.3 (ii) diperoleh:

∫

∫

−

++−

=++

na

nana

a

aa

nanaaa

dxxgxf

dxxgxf

gfAgfA

0

2

022

22

.......))()((

...))()((

),(...),(

(3)

karena daerah ),(),...,,(),,( 2211 nanaaaaa gfRgfRgfR setara,

berdasarkan Definisi 2.2 dan Definisi 2.3 (i), maka Persamaan (3) ditulis sebagai:

∫

∫

−

++

−

=++

na

na

na

n

a

aa

nanaaa

dxxaagx

aaf

aa

dxxaagx

aaf

aa

gfAgfA

0

11

11

1

2

0 2

11

2

11

1

2

22

.......)()(

...)()(

),(...),(

(4)

Dengan melakukan perubahan variabel yaitu:

1,...,1,1

1 −==+

nixaa

ui

i maka Persamaan (4) dapat

pula ditulis sebagai:

( )

( )∫

∫

−−− −

++−

=++

1

01

11111

1

1

01

1

21111

1

2

22

.......)()(

...)()(

),(...),(

a

nn

nanan

a

aa

nanaaa

duaauguf

aa

duaauguf

aa

gfAgfA

( )

( )∫

∫

−

++−=

1

0112

1

2

1

0112

1

22

.......)()(

...)()(

a

aan

a

aa

dxxgxfaa

dxxgxfaa

( )dxxgxfa

aaaa

aan ∫ −

+++=

1

0112

1

223

22 )()(

... (5)

Karena ),...,,( 21 naaa adalah solusi Persamaan (1),

maka diperoleh:

).,...,,(...

),...,,(...

212

122

32

2

2122

32

22

1

nn

nn

aaagaaaa

aaagaaaa

−=+++

++++=

Sehingga Persamaan (5) dapat dituliskan sebagai:

∫ −

=++

−1

0112

1

),...,2,1(21

22

))()((

),(...),(a

aaa

naaaga

nanaaa

dxxgxf

gfAgfA

),( 1121

),...,2,1(21

aaa

naaaga gfA−=

),(),( 1121

),...,2,1(11 aa

anaaag

aa gfAgfA −=

Dengan demikian diperoleh persamaan:

),(),...,,(

),(...),(),(

1121

21

2211

aan

nanaaaaa

gfAa

aaag

gfAgfAgfA

+

++=

Teorema 4.2 Misalkan niaaaa in ≤≤> 1,0,),...,,( 21 adalah solusi

Persamaan (2). Misalkan daerah ),(),...,,(),,( 2211 nanaaaaa gfRgfRgfR setara,

maka :

).,(),(1

2),...,2,1(

iaia

n

j ianaaaf

jaja gfAgfA∑=

=

Jika 0),( ≠iaia gfA maka

∑=

=n

jjaja

iagiafAia

n gfAaaaf1

),(

221 ).,(),...,,(

Bukti: Berdasarkan Definisi 2.3 (ii) diperoleh:


142

∫

∫

∑

−

++−=

++==

na

nana

a

aanana

aa

n

jjaja

dxxgxf

dxxgxfgfA

gfAgfA

0

1

011

111

.......))()((

...))()((),(

...),(),(

(6)

karena daerah ),(),...,,(),,( 2211 nanaaaaa gfRgfRgfR setara,

berdasarkan Definisi 2.2 dan Definisi 2.3 (i), maka Persamaan (4.6) ditulis sebagai:

∫

∫∑

−++

−=

=

na

n

iia

n

iia

i

n

ai

iai

iai

n

jjaja

dxxaagx

aaf

aa

dxxaagx

aaf

aagfA

0

1

0 11

1

1

)()(...

)()(),(

(7)

Dengan melakukan perubahan variabel yaitu:

1,...,1, −== njxaa

uj

ij maka Persamaan (7) dapat

pula ditulis sebagai:

( )

( )∫

∫∑

−++

−==

ia

ni

nniania

i

n

ia

iiaia

i

n

jjaja

duaauguf

aa

duaauguf

aagfA

0

01

111

1

1

)()(...

)()(),(

( ) ( )∫∫ −++−=ia

iaiai

nia

iaiai

dxxgxfaadxxgxf

aa

02

2

02

21 )()(...)()(

( )dxxgxfa

aa ia

iaiai

n ∫ −++

=0

2

221 )()(

... (8)

Karena ),...,,( 21 naaa adalah solusi Persamaan (2),

maka: 22

22

121 ...),...,,( nn aaaaaaf +++= .

Dengan demikian Persamaan (4.8) dapat dituliskan sebagai:

).,(),(1

2),...,2,1(

iaia

n

j ianaaaf

jaja gfAgfA∑=

=

Jelas bahwa, bila 0),( ≠iaia gfA maka:

∑=

=n

jjaja

iagiafAia

n gfAaaaf1

),(

221 ).,(),...,,(

Akibat 4.3 Misalkan niaaaa in ≤≤> 1,0,),...,,( 21 adalah solusi

Persamaan (4.2). Misalkan daerah nkgfR kaka ≤≤1),( semuanya setara. Luas daerah

),( kaka gfR yaitu 0),( =kaka gfA jika dan hanya

jika 0),(1∑=

=n

jjaja gfA .

Bukti: )(⇒ Misalkan 0),( =kaka gfA . Maka menurut

Teorema 4.2, 0),(1∑=

=n

jaa jj

gfA .

)(⇐ Misalkan 0),(1∑=

=n

jjaja gfA dan andaikan

0),( ≠kaka gfA untuk suatu k dengan nk ≤≤1 .

Maka menurut Teorema 4.2 haruslah

0),...,,( 21 =naaaf hal ini bertentangan dengan

Persamaan (2). Oleh karena itu haruslah 0),( =kaka gfA .


Persamaan (2). Misalkan daerah njgfR jaja ≤≤1),(

semuanya setara. Jika 0),( ≠kaka gfA dan

0),( ≠iaia gfA untuk sebarang i,k dengan

nkni ≤≤≤≤ 1,1( maka ),(

2

),(

2

kagkafAka

iagiafAia =

Bukti: Misalkan 0),( ≠kaka gfA untuk nk ≤≤1 dan daerah

),( jaja gfR semuanya setara untuk nj ≤≤1 maka

menurut Teorema 4.2 berlaku:


143

∑=

=n

jjaja

kagkafAka

n gfAaaaf1

),(

221 ).,(),...,,(

(9) dengan cara yang sama berlaku pula:

∑=

=n

jjaja

iagiafAia

n gfAaaaf1

),(

221 ).,(),...,,( (10)

Dari Persamaan (9) dan (10) diperoleh:

∑=

=n

jjaja

iagiafAia gfA

1),(

2),(

∑=

n

jjaja

kagkafAka gfA

1),(

2),(

Dengan demikian diperoleh:

),(

2

),(

2

kagkafAka

iagiafAia =


Persamaan (4.2). Misalkan daerah ),( jaja gfR dan

daerah ),( jaja GFR setara untuk nj ≤≤1 . Jika

0),( ≠kaka gfA dan 0),( ≠kaka GFA maka

),(),(

),(),(

kaGkaFAkagkafA

iaGiaFAiagiafA

= untuk sebarang i,k dengan

nkni ≤≤≤≤ 1,1( .

Bukti: Misalkan 0),( ≠kaka gfA untuk nk ≤≤1 dan daerah

),( jaja gfR setara untuk nj ≤≤1 maka menurut

Akibat 4.4 diperoleh:

),(

2

),(

2

kagkafAka

iagiafAia = (11)

Misalkan 0),( ≠kaka GFA untuk nk ≤≤1 dan daerah

),( jaja GFR setara untuk nj ≤≤1 maka menurut

Akibat 4.4 diperoleh:

),(

2

),(

2

kaGkaFAka

iaGiaFAia = (12)

Dari Persamaan (11) dan Persamaan (12) diperoleh:

),(),(

22iaia

kagkafAka

i gfAa =

),(),(

2

iaiakaGkaFA

ka GFA=

Dengan demikian diperoleh:

),(),(

),(),(

kaGkaFAkagkafA

iaGiaFAiagiafA

= .

5. Kesimpulan

1. Telah diperoleh hasil perluasan teorema Pythagoras

di 2R dengan menggunakan hubungan kesetaraan pada luas daerah. Hasil yang diperoleh dari penggunaan hubungan kesetaraan tersebut adalah : Jika pada setiap sisi segitiga siku-siku terdapat daerah yang beraturan maupun yang tidak beraturan dan setara, maka terdapat hubungan diantara luas daerah tersebut dengan persamaan teorema pythagoras, yaitu:luas daerah yang terkait dengan sisi miring segitiga siku-siku, namakan sisi c merupakan penjumlahan dari luas daerah yang terkait dengan panjang sisi siku-siku, namakan sisi a dan sisi b.

2. Dari hasil perluasan teorema Pythagoras di 2R diperoleh bentuk perumusan umum perluasan

teorema Pythagoras, yaitu di nR seperti dijelaskan dalam teorema 4.1 dan teorema 4.2.

Daftar pustaka [1] Anton, H., Calculus with Analytic Geometry, 1988, John

Wiley. New York [2] Anton, H., Aljabar Linier Elementer, Edisi kelima, 1994,

Erlangga, Jakarta [3] Ausri, A., Pengembangan Sifat Pythagoras dengan

Menggunakan Hubungan Kesetaraan, 1997, Jumpa 6 : 83-90

[4] Clay, James R and Yuen Fong, Generalization of Pythagorean Theorem, 1995, Sea Bull. Math. 19:19-26

[5] Herstein, I.N., 1975. Topics In Algebra. 2thed. 1975, John Wiley. New York


Metode Diagonally Implicit Runge-Kutta Dalam Penyelesaian Model Matematika Sistem Cardiovascular

Fanani Haryo Widodo dan Awal Isgiyanto



Abstrak - Pokok pembahasan paper ini adalah tentang penentuan metode penyelesaian model matematika untuk sistem cardiovascular. Model tersebut disajikan sebagai ssistem non-linier persamaan diferensial biasa. Karateristik sistem tersebut adalah sangat “kaku” sebagaimana ditunjukkan oleh nilai eigennya yang riil negatif. Dalam pemilihan metode yang tepat untuk menyelesaikan model tersebut, tingkat akurasi dan efisiensi menjadi ukuran yang penting. Dalam hal ini, telah dijelajahi klas metode implisit Runge-Kutta. Fakta menunjukkan bahwa klas metode tersebut sangat “mahal” karena klas metode tersebut pada umumnya melibatkan persamaan aljabar yang ketat dan rumit. Dalam studi yang telah dilakukan, ditemukan bahwa metode Diagonally Implicit Runge-Kutta (DIRK) adalah tepat untuk menjawab model tersebut. Dalam hal ini digunakan Butcher array yang dihasilkan oleh penelitian Kevin Burrage. Yang menarik dari studi ini adalah bahwa matrik Jacobian memiliki struktur blok diagonal yang terklasifikasi dalam tiga blok sedemikian sehingga metode penyelesaian yang didapat memiliki potensi bagi penyelesaian dengan menggunakan algoritma paralel. Kata Kunci : DIRK; model matematika sistem cardiovascular; Butcher array Kevin Burrage; matrik berstruktur blok

diagonal.

1. Pendahuluan Sistem kontrol memainkan peranan yang sangat penting dalam perkembangan teknologi dan kemajuan peradaban modern. Hal ini banyak dijumpai di semua sektor industri seperti sistem transportasi, pengendalian kualitas produk suatu pabrik, sistem pembangkit energi dll. Demikian juga problem-problem di dalam pengendalian sistem sosial dan ekonomi dalam kehidupan kita dan pengendalian sistem sirkulasi darah dalam sistem cardiovascular manusia atau hewan dapat dikaji melalui pendekatan teori kontrol. Didukung oleh kemajuan dan meluasnya penggunaan teknologi komputer, peranan sistem kontrol telah secara dramatis mengalami kemajuan khususnya dalam analisis dan desain. Sistem cardiovascular pada tubuh manusia telah menjadi suatu studi yang intensif di dalam sistem kontrol khususnya dalam suatu kerangka yang terkait dengan problem Mayer wave, suatu osilasi dengan periode yang panjang dalam tekanan darah arteri. Serangan Mayer wave sering muncul ketika seseorang berada pada kondisi abnormal seperti kekurangan darah, akibat pendarahan

yang berlebihan, atau karena perubahan yang mendadak dalam suplai darah ke bagian-bagian tubuh [3]. Dalam suatu masyarakat yang berpenduduk padat dengan lingkungan yang tercemar, masalah Mayer wave mungkin cukup potensial sebagai penyakit masyarakat dikarenakan oleh seringnya ditemukan kondisi-kondisi abnormal seperti kemacetan lalu lintas, dan situasi penuh sesak oleh kerumunan manusia dengan udara yang tercemar. Pada kondisi tubuh yang kritis serangan Mayer wave dapat menyebabkan implikasi psiologis yang serius seperti pingsan [1]. Karenanya studi khusus tentang masalah ini akan memberikan kontribusi pada upaya penyelesaian penyakit masyarakat. Dari banyak eksperimen yang difokuskan pada Mayer wave banyak ditemukan fitur yang menjelaskan gejala sekitar Mayer wave. Meskipun begitu sebab timbulnya Mayer wave tetap merupakan problem tak terselesaikan [13]. Lebih jauh problem pengembangan metode penyelesaian untuk model matematika dari sistem cardiovascular memperoleh perhatian pada beberapa tahun terakhir ini. Abbiw-Jackson melakukan studi tentang sisem cardiovascular yang dimodelkan secara

Fanani / Jurnal Gradien Vol. 2 No. 1 Januari 2006 : 144-151

145

matematika dalam bentuk suatu sistem tiga persamaan diferensial biasa nonlinier [1]. Mereka menganalisa model tersebut menggunakan metode stabilitas Liapunov dan teori Hoft bifurcation dengan kesimpulan bahwa kenaikan di dalam gain dari baroreflex feedback loop memungkinkan timbulnya serangan osilasi, sedangkan perubahan dalam parameter-parameter model tidak mempengaruhi kestabilan dari titik keseimbangan. Penelitian ini mencoba mencari penyelesaian problem matematika dari sistem cardiovascular sebagaimana diberikan oleh [1] melalui pendekatan komputasi numerik. Basaruddin dalam penelitiannya mengajukan suatu metode untuk menjawab model demikian yang didasarkan pada kelas BDF [9]. Namun sebegitu jauh banyak penemuan menunjukkan bahwa metode BDF hanya tepat untuk sistem yang tidak begitu kaku, sementara model dari sistem cardiovascular mempunyai indikasi cukup stiff sebagaimana hasil yang diperoleh dari penelitian Abbiw-Jackson yang mengatakan bahwa perubahan-perubahan parameter tidak mempengaruhi kestabilan titik kesetimbangan [1]. Karenanya untuk menyelesaiakan model tersebut diperlukan daerah kestabilan yang luas. Sifat yang demikian dimiliki oleh sifat metode Diagonally Implicit Runge-Kutta. Pada sisi lain kelemahan pada metode Implicit Runge-Kutta adalah perhitungannya menggunakan persamaan aljabar yang rumit dan komplek. Dengan kata lain metode ini cukup mahal atau kurang efisien. Peneliti berusaha mengembangkan metode DIRK untuk menjelesaikan problem tersebut.

2. Sistem Cardiovascular dan Mayer Wave Sistem cardiovascular adalah sistem sirkulasi darah. Pada manusia dan mamalia tingkat tinggi sistem tersebut dibagi menjadi dua, yaitu sistem siskulasi besar (sistemik) yang mensuplai darah dari jantung keseluruh organ tubuh kemudian kembali ke jantung, dan sistem sirkulai kecil (pulmonary) yang membawa darah dari jantung ke paru-paru untuk pertukaran gas kemudian kembali ke jantung. Perjalanan darah pada sirkulasi sistemik, darah mengalir dari bilik kiri (atrium sinistrum) masuk ke serambi kiri (ventriculus sinistrum), kemudian melalui arteri utama (aorta) masuk ke jalinan kapiler organ tubuh. Sedangkan perjalanan darah pada sirkulasi kecil, darah mengalir dari bilik kanan (atrium dextrum) masuk ke serambi kanan

(ventriculus dextrum), kemudian melalui arteri paru-paru (pulmonalis) masuk ke jalinan kapiler paru-paru kemudian mengalir kembali melalui vena paru-paru masuk ke bilik kiri (atrium sinistrum). Dari sistem sirkulasi di atas terlihat bahwa ada empat volume darah yang memungkinkan dipilih sebagai variabel model sistem cardiovascular, yaitu volume darah pada arteri sistemik (systemic arteries) SAv , volume

darah di dalam vena sistemik (systemic veins) SVv ,

volume darah di dalam vena paru-paru (pulmonary veins)

PVv , dan volume darah di dalam arteri paru-paru

(pulmonary arteries) PAv . Akan tetapi penelitian yang

dilakukan oleh Abbiw-Jackson menunjukkan bahwa variabel volume darah di dalam arteri paru-paru merupakan volume terkecil. Kecuali itu jika dilihat dari fungsinya, darah yang berada di dalam arteri paru-paru hanya berfungsi sebagai transport pembuangan gas, sehingga ditetapkan untuk diabaikan. Jadi variabel pada penelitian ini SAv , SVv dan PVv .

3. Model Matematika Sistem Cardiovascular

Model matematika sistem cardiovascular meliputi sistem tiga persamaan diferensial biasa sebagai berikut [1]:

SCS

D

PV

PVL

SVS

SV

SAS

SASACR

VC

vFCCR

vCR

vdt

dv−++−=

( )DSVSV

R

SVSSAS

SASV VvCFC

CRCRv

dtdv

−

+−=

1

+−

−−−=

PV

L

PVPPAP

PVSVSAPVCFC

CRCRvvvV

dtdv 10

dimana CSA = Compliance in systemic arteries, CPA = Compliance in pulmonary arteries, CPV = Compliance in pulmonary veins, CSV = Systemic venous compliance, RS = Systemic capillary resistance, RP = Pulmonary resistance, CR = Compliance in right heart, CL = Compliance in left heart, F = Heart rate, VD = Systemic venous unstressed volume, volume of vessel at pSV = 0, pSV = Systemic vein pressure, VSA = volume of blood in systemic arteries, VSV = volume of blood in systemic veins, VPV = volume of blood in pulmonary veins. Dengan mensubstitusi persamaan F, RS, CSV, dan VD yang menyatakan model sistem Baroreflex control, yaitu:


146

• ( ) ( )( ) ( ) 2

121 1 F

vV

VFFBFF

nSA

nc

nc

n ++

=+−=

dimana F1 dan F2 konstan, dan F2 adalah nilai F untuk Bn = 1.

• systemic resistance Rs dimodelkan sebagai :

( ) ( )( ) ( ) 2

121 1 R

vVVRRBRR n

SAn

C

nC

nS ++

=+−=

• systemic venous compliance CSV dimodelkan sebagai:

( )( ) ( ) 2

121 C

vV

vCCBCC n

SAn

C

nSA

nSV ++

=+=

• systemic venous unstressed volume VD dimodelkan sebagai :

DVD =

( )( ) ( ) 2

121 C

vV

vDDBDV n

SAn

C

nSA

nD ++

=+=

maka didapat model nonlinier dari sistem cardiovascular sebagai berikut:

1. ( ) SA

nC

nSASA

SAn

CnSASA

CVRRVCRVVV

dtdV

212

1

++

+−=

+

( )PV

nC

nSAPV

PVLn

CnSA

CVVCCCVFFVF

+

+++ 212

( )( ) ( )( )( ) ( ) n

CnSA

nC

nSA

nC

nSASV

nC

nSASV

VCRRVVCRCCRRVRCCVDVDDVVVV

2221222121

2221

221

+++++++

−+−++

2. ( ) SA

nC

nSASA

SAn

CnSASV

CVRRVCRVVV

dtdV

212

1

++

+=

+

( ) ( ) ( )( ) ( ) 221212

221

CVVCCVRRVRVVVVVVDVDD

nC

nSA

nC

nSA

SVn

CnSA

nC

nSA

nC

nSA

++++

+−++++

( ) ( ) n

CnSA

SVRn

CnSA

VCVCCVCVFFVF

221

212

++

++−

( ) ( ) ( )( ) n

CnSA

nC

nSA

nC

nSA

nC

nSA

VCVCCVVVDVDDVFFVF

221

221212

+++

+++++

3. PAP

PVSVSAPVCR

VVVVdt

dV −−−= 0

( ) ( ) PV

PVn

CnSA

Ln

CnSA

PVPV

CVVCVFFVF

CR

+

+++− 2121

Sedangkan parameter yang relevan terhadap Mayer wave ada 16 buah sebagaimana ditunjukkan dalam Tabel 1 di bawah [1][12][15][17].

Tabel 1. Nilai parameter yang digunakan dalam penelitian. Simbol Nilai

CSA 0,01 lt/mmHG CPA 0,00667 lt/mmHG CPV 0,08 lt/mmHG RP 1,79 mmHG CR 0,035 lt/mmHG CL 0,014 lt/mmHG

F1 = 80 beats/min F F1 = 40 beats/min

R1 20 lt/min R2 7,5 lt/min D1 3,0 lt D2 0,5 lt C1 1,0 lt/mmHG C2 0,25 lt/mmHG VC 1,0 lt VD 2,0 lt

4. Masalah Nilai Awal

Masalah nilai awal disajikan sebagai :

( )( ) 00

0,,yty

ttytfy=

>=′

dimana y adalah suatu vektor di dalam Rm dan mm RRf →: adalah suatu fungsi nonlinier. Persamaan

ini juga dapat dinyatakan dalam bentuk : ( )

( ) 00

0,yty

ttyfy=

>=′

Menyelesaikan model ini adalah lebih mudah jika digunakan metode Runge-Kutta, karena tidak diperlukan peubah waktu secara khusus. Ini dilakukan dengan menampilkan peubah waktu sebagai satu komponen dalam sistem [9]. Masalah nilai awal tersebut dikatakan kaku jika terdapat perbedaan yang besar antara nilai-nilai eigen λ dari matrik Jacobian permasalahan tersebut. Masalah syarat awal

( ) ( ) ,, 00 yxyyfy ==′ y vektor berdimensi m, dengan

fx merupakan titik ujung dari integrasi. Ada masalah

syarat awal yang solusi eksaknya dapat ditentukan dengan


147

mudah, akan tetapi banyak juga masalah syarat awal yang solusi eksaknya tidak mungkin ditemukan [9]. Sedangkan Alhadi [2] menyatakan bahwa ( ) ( )( ),xyfxy =′

( ) mm RRfyxy →= :,00 dikatakan kaku bila

( ) 10 >>− Lxx f dimana L adalah Konstanta Lipschitz.

Atau jika nilai eigen dari matrik Jacobiannya negatif dan terdapat perbedaan yang besar antara nilai-nilai eigennya. Jika pada suatu kasus grafiknya turun, maka akan terjadi penurunan dengan laju sebesar iλ . Semakin besar iλ

berarti semakin besar penurunannya sebagai akibat perbedaan yang besar antara nilai iλ , sehingga terjadi

perubahan yang besar dan penurunan yang cepat pada suatu interval tertentu. Pada interval dengan penurunan yang cepat diperlukan strategi pengambilan ukuran langkah (step size) yang kecil, sedangkan pada selang yang solusinya berjalan lambat, maka digunakan ukuran langkah yang lebih besar. Tidak mudah untuk mengetahui nilai L awal pada sebagian besar kasus, sehingga sulit untuk menentukan kapan nilai ( )Lxx f 0− termasuk kategori besar [14]. Pendekatan lain

menyatakan bahwa kekakuan dapat dideduksi dari sifat sistem yang dimodelkan. Hainer dan Wanner menyatakan bahwa masalah nilai awal yang kaku memiliki solusi numerik dimana gerakan lambat yang mulus akan diganggu oleh solusi-solusi yang bergerak cepat didekatnya [11]. Untuk menyelesaikan problem yang kaku dapat digunakan metode Implisit Runge-Kutta (IRK), dengan berbagai variasinya. Dalam penelitian ini metode yang digunakan adalah metode Diagonally Implicit Runge-Kutta (DIRK).

5. Metode Runge-Kutta Umum Solusi menggunakan pendekatan numerik dapat dilakukan dengan metode Implisit Runge-Kutta (IRK). Ada beberapa metode Implisit Runge-Kutta yang dapat digunakan, metode IRK Labatto [9], metode IRK Radau [9] atau metode Diagonally Implicit Runge-Kutta (DIRK). Secara umum, metode IRK pada tahap-s untuk persamaan diferensial biasa ( )ytfy ,=′ , dapat ditulis dalam bentuk:

( ) siYhctfahyY jjn

s

jijni ≤≤++= −

=− ∑ 1,,1

11

( )iin

s

iinn Yhctfbhyy ,1

11 ++= −

=− ∑

Metode tersebut dapat dinyatakan juga dalam notasi : c1 a11 a12 … a1s c2 a21 a22 … a2s

… … … … … cs as1 as2 … ass b1 b2 … as

sehingga dapat ditulis ∑=

==s

jsiiji ac

1,...,1, si ,...,1= [8].

Sedangkan metode Diagonally Implicit Runge-Kutta (DIRK) 4-tahap sebagai ODE-solver adalah :

1. ( )∑=

+=4

1jjijni VfahvV 4,3,2,1; =i

2. ( )∑=

+ +=4

11

iiinn Vfbhvv

Kelas dari metode ini diklasifikasikan sesuai integrator yang dinyatakan sebagai Butchers array :

c A bT

Dimana c dan b adalah vektor yang panjangnya s dan A adalah suatu matrik berordo s x s. Basaruddin telah menunjukkan penggunaan Butchers array diatas dengan rumus Runge-Kutta [9].

( )

( )i

s

iinn

j

s

jijni

Yfbhyy

siYfahyY

∑

∑

=+

=

+=

=+=

11

1...1,

Ambil ( ) smTTs

Ti RYYY ∈= ,..., dan pilih

( ) ( ) ( )( )TTs

Ti YfYfYF ,...,= , maka jika metode Runge-

Kutta dituliskan dalam notasi tensor, dinyatakan sebagai : ( ) ( )( ) ( )YFIbhyy

YFIAhyeY

mT

nn

mn

⊗+=

⊗+⊗=

+1

Pada notasi ini tampak jelas bahwa sistem persamaan di atas memerlukan penyelesaian di setiap tahap. Dalam kasus f adalah suatu fungsi nonlinier, persamaan tersebut


148

menunjukkan munculnya sistem persamaan nonlinier. Dalam kasus ini, metode iterasi seperti metode Newton boleh digunakan untuk menentukan Yi. Jika pada akhir suatu iterasi Yi diganti oleh Yi + di iterasi Newton akan memerlukan vektor di yang baru agar memenuhi : ( ) Ψ=−⊗ dhIII ms

dimana J adalah blok matrik orde sm, dengan elemen (i,j) adalah aijf(Yj), i,j =1, … ,s dan

( ) ( )( ) siYfahyY

ddd

j

s

jijnii

TTs

TTTs

T

,...,1;

,...,,,...,

1

11

=++−=Ψ

ΨΨ=Ψ=

∑=

Bentuk persamaan ini juga mudah digunakan dengan aproksimasi matriks Jacobian ( )jYf menggunakan

( )nyf , yang telah dihitung pada tahap n.

6. Algoritma DIRK Untuk Model Matematika Sistem

Cardiovascular Metode Diagonally Implicit Runge-Kutta (DIRK) 4-tahap adalah :

( ) 43,2,1;..14

1andjVfahvV j

jijni =+= ∑

=

(1)

( ) 43,2,1;.24

1 andiVfbhvv ii

inn =+= ∑+ ( 2)

dengan Vi, vn dan f adalah vektor tiga demensi. Pengertian 4 tahap menunjuk pada proses untuk mendapatkan nilai vn melalui 4 tahap, yaitu V1, V2, V3 dan V4. Karenanya Persamaan (1) merupakan sistem persamaan yang terdiri dari 12 persamaan nonlinier. Untuk menentukan nilai V1...V4 yang memenuhi sistem persamaan tersebut dipakai iterasi Newton [18]. Untuk menerapkan metode Newton pada (1) maka Persamaan (1) tersebut ditulis sebagai:

( ) ( ) 43,2,1;,,,4

14321 andiVfahvVVVVVF j

jijnii =++−= ∑

=

(3) Algoritma Newton tersebut adalah sebagai berikut:

1. Choose a starting point Vi* 2. Calculate Fi(x*) 3. Stop if all components are sufficiently small 4. Calculate J(x*)

5. Solve for d the n x n system J(x*)d = - f(x*) ; x = x* + d

6. Set x* = x 7. Repeat from step 2

Perhitungan ki untuk persamaan U = AU + d(t) dapat diturunkan menjadi :

( ),)(1

hctdkahUAk ij

i

jiji +++= ∑

=

dengan i = 1,2,3,4. untuk selanjutnya didapat

( )( )( )( )

++++++++

==

4

3

2

1

hctdUAhctdUAhctdUAhctdUA

rkkN

yang merupakan vektor berukuran 4nx1 dengan N adalah

−−−−−−

−−−

=

hAaIhAahAahAaIhAahAa

hAaIhAahAaI

N

)4,4(0)2,4()1,4(0)3,3()2,3()1,3(00)2,2()1,2(000)1,1(

yang merupakan matrik berukuran 4nx4n, sedangkan I dan 0 adalah matrik identitas dan matrik nol yang berukuran n x n, dan k adalah vektor berukuran 4nx1. Algoritma DIRK 4-tahap tersebut adalah [16]: Langkah 1. Hitung rk = [rk1 rk2 rk3 rk4], dengan rkI =

AU + d(t+h*ci) adalah vektor berukuran nx1, i = 1,2,3 dan 4.

Langkah 2. Hitung Nij dengan Nii =I - aii*h*A, Nij = aij*h*A, i = 1,2,3 dan 4, j = 1, 2, i – 1.

Langkah 3. Selesaikan Nk = rk untuk mencari k = [k1 k2 k3 k4]T, dengan k1, k2, k3, k4 masing-masing berukuran nx1.

Langkah 4. Hitung solusi Um+1 = Um + ii

i bkh∑=

4

1,

dengan Um+1 dan Um masing-masing berukuran nx1.

7. Solusi Model

Untuk menyelesaikan model tersebut, ODE-solver yang digunakan adalah rumus Diagonally Implicit Runge-Kutta. Dalam perhitungannya ditunjukkan, bahwa modelnya melibatkan sejumlah parameter yang nilainya merupakan hasil kajian literature sebagaimana disajikan di atas.


149

Matriks Jacobian J(f) yang sesuai dengan tiga peubah ( =ASv . volume darah di dalam systemic arteries, =SVv

volume darah di dalam systemic veins dan =PVv volume

darah di dalam pulmonary veins) disajikan sbb:

Ambil vektor

=

PV

SV

SA

VVV

v . Maka ( )vfv =′ dan

=′

dtdV

dtdV

dtdV

v

PV

SV

SA

; ( )( )( )( )

=

vfvfvf

vf

3

2

1

( )

∂∂

∂∂

∂∂

∂∂

∂∂

∂∂

∂∂

∂∂

∂∂

=

PVSVSA

PVSVSA

PVSVSA

Vf

Vf

Vf

Vf

Vf

Vf

Vf

Vf

Vf

fJ

333

222

111

Metode implicit Runge-Kutta sangat layak untuk menyelesaikan masalah tersebut, karena metode ini akan mempunyai daerah stabilitas yang luas ('large stability region'). Alhadi [2] dalam penelitiannya menunjukkan bahwa integrator mempunyai sifat-sifat seperti yang disampaikan oleh Burrage, K. [10] 0 0 1 1/2 1/2 3/2 5/8 3/8 1/2 1 7/18 1/3 -2/9 1/2 7/18 1/3 -2/9 1/2 Penerapan Butcher array pada rumus DIRK menghasilkan model : 1. nvV =1

( ) ( )

++= 212 2

121 VfVfhvV n

( ) ( ) ( )

+++= 3213 2

183

85 VfVfVfhvV n

( ) ( ) ( ) ( )

+−++= 43214 2

192

31

187 VfVfVfVfhvV n

2. ( ) ( ) ( ) ( )

+−++=+ 43211 2

192

31

187 VfVfVfVfhvV nn

Substitusi f ke rumus RK menghasilkan

1).

nPV

SV

SA

PV

SV

SA

vvv

VVV

=

1

( ) ( )

++

=

21

2

,,21,,

21

PVSVSAPVSVSA

nPV

SV

SA

PV

SV

SAVVVfVVVfh

vvv

VVV

( ) ( ) ++

=

21

3

,,83,,

85

PVSVSAPVSVSA

nPV

SV

SA

PV

SV

SAVVVfVVVfh

vvv

VVV

( ) ],,21

3PVSVSA VVVf+

( ) ( )21

4

,,31,,

187

PVSVSAPVSVSA

nPV

SV

SA

PV

SV

SAVVVfVVVfh

vvv

VVV

+

+

=

( ) ( ) ],,21,,

92

43 PVSVSAPVSVSA VVVfVVVf +−

2)

( ) ( )21

1

,,31,,

187

PVSVSAPVSVSA

nPV

SV

SA

nPV

SV

SAVVVfVVVfh

vvv

VVV

+

+

=

+

( ) ( ) ],,21,,

92

43 PVSVSAPVSVSA VVVfVVVf +−

dimana : ( )

=

dtdV

dtdV

dtdV

VVVf

PV

SV

SA

PVSVSA ,,

Analisis pada penelitian ini menunjukkan bahwa model yang didapatkan adalah kaku, hal ini terlihat dari hasil nilai ciri (eigen values) dari matriks Jacobian :

( )

−−−

−=

740235.16775699967.832429676.6160767618368.72159440.303

9999931.7620761904695.5092704.773J f

pada sistem persamaan diferensial biasa, yaitu :

202184461.2108783187.51202184461.2108783187.51

841602.846

3

2

1

−−=+−=

−=

λλλ

Hasilnya menunjukkan bahwa sebaran nilai cirinya negatif dan nilai mutlaknya berbeda cukup lebar antara nilai ciri yang satu dengan yang lain. Ini menunjukkan bahwa


150

modelnya adalah kaku (stiff). Hasil lain adalah bahwa matriks Jacobian ordo 12 x 12 J(FF) menunjukkan suatu bentuk block diagonal seperti di bawah ini:

Dalam hal ini penggunaan parameter dan nilai awal sesuai dengan literatur [1][12][15][17]. Algoritma yang dibangun melalui DIRK-4 tersebut telah menghasilkan nilai konvergensi untuk vektor volume darah (VSA, VSV, VPV) sebagai berikut:

8. Kesimpulan

Dari kajian di atas dapat dirumuskan beberapa poin kesimpulan sebagai berikut: Untuk parameter-parameter yang diperoleh dari hasil survey literatur, model matematika dari sistem cardiovascular merupakan problem yang stiff seperti ditunjukkan oleh eigenvalue bernilai riil dari matrik Jacobiannya yang selalu negatif.

Meskipun titik konvergensi dari vektor volume darah bisa dicapai dengan algoritma yang dikembangkan namun akurasi konvergensi tersebut masih perlu diuji dengan mengembangkan estimasi kesalahan dalam algoritmanya. Matrik Jacobian yang diturunkan dari DIRK 4-tahap berstruktur blok diagonal. Karenanya problem tersebut dapat diselesaikan dengan metode paralel.

Daftar Pustaka

[1] Abbiw-Jackson, S.M and William F. Langford, Gain-induced Oscillations in Blood Pressure, 1998, Journal of Mathematical Biology, 37, 203-234

[2] Alhadi, B., Analisis Kinerja Integrator Singly Diagonally Implicit Runge-Kutta Method (SDIRK) dalam Menyelesaikan Sistem ODE yang Kaku, 2000, Jurnal Penelitian UI. (submitting).

[3] Anderson, B.A., R.A. Kenney, and E. Neil, The role of the Chemoceptors of the Carotid and Aortic Regions in the Productions of the Mayer Waves, 1950, Acta. Physiol. Scand., 20, 203-220

[4] Anonim, An Adaptive Method of Lines for Non-Linear System of Convection-Diffusion-Reaction Equations Arising from Contaminant Transport in Groundwater Pollution, 1998, Jakarta, Faculty of Computer Science - University of Indonesia.

[5] Anonim, Parallel Hepscotch Methods for Shallow Water Equations, 1999, Makara, Depok: No.6 Seri B. 2 Januari 1999. p. 60-68

[6] Anonim, Prekondisi Otomatis Metode Conjugate Gradient untuk Penyelesaian Persoalan Optimisasi, 2000, Makara, no. 7 seri B Mei, hal. 49-55.

[7] Anonim, Stiff Differential Equations Solved by Radau Methods, 1998, Switzerland : Dept. de Mathematiques, Universite de Geneve

[8] Ascher, U. M. and L. R. Petzold, Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations, 1998, Philadelphia : Siam.

[9] Basaruddin, T., Parallel Methods for Ordinary Differential Equations, Parallel Simulation Environment and Its Application Groundwater Flow, 2000a, (Monitoring and Self Evaluation Report Graduate Team Research Grant BATCH IV), Jakarta: Faculty of Computer Science-University of Indonesia.

[10] Burrage, K. et. al., A-Four Stage Index 2 Diagonally Implicit Runge-Kutta Method, 1999, Laporan riset dari CAPEC (Computer Aided Process Engineering Centre), Brisbane : University of Quensland.


151

[11] Hairer, E and G Wanner, Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems. Springer Seriea in Computational Mathematics, 1990, Springer Verlag.

[12] Hoppensteadt, F.C. and C.S. Peskin, Mathematics in Medicine and the Life Sciences, 1992, New York : Springer-Verlag.

[13] Penaz, J., Mayer Wave: History and Methodology, 1978, Automedica, 2, 135-141.

[14] Sandys, D. Y. dan T. Basaruddin, Evaluasi Kinerja Radau IIA Dalam Menyelesaikan Persamaan Diferensial Kaku, Parallel Simulation Environment and Its Application Groundwater Flow, 1999, (Monitoring and Self Evaluation Report Graduate Team Research Grant BATCH IV), Jakarta: Faculty of Computer Science-University of Indonesia.

[15] Schnittinger, I. et. al., Standardized intracardiac measurements of two-dimensional echocardiography, 1983, J Am Coll Cardiol. 5: 934.

[16] Suhartono, Efek Diskritisasi Spatial Methode Galerkin Semi Deskrit dan Metode Beda Hingga Terhadap Kinerja Metode Runge-Kutta Implisit Diagonal, 1999, Thesis, Depok : Fasilkom Universitas Indonesia.

g ISSN 0216-2393 GRADIEN · PDF fileKimia 4. Penentuan Adanya ... Abstrak - Telah dilakukan penelitian mengenai pemodelan fisis aplikasi metode geolistrik tahanan jenis 2 ... ciri

Documents