12 0 599 KB
INTEGRAL NUMERIK ATURAN SIMPSON π/π Diajukan untuk memenuhi salah satu tugas mata kuliah Metode Numerik Dosen Pengampu: Darta, S.Pd., M.Pd. / Thesa Kandaga, S.Si., M.Pd.
Oleh : Jayanti Purwati
(175050004)
Nida Aulia Pangesti Limbangan
(175050025)
Risya Fadilah
(175050031)
PROGAM STUDI PENDIDIKAN MATEMATIKA FAKULTAS KEGURUAN DAN ILMU PENDIDIKAN UNIVERSITAS PASUNDAN BANDUNG 2020
Dengan menyebut nama Allah yang Maha Pengasih lagi Maha Penyayang
KATA PENGANTAR Puji dan Syukur kita panjatkan kehadirat Allah SWT karena berkat rahmat dan hidayah-Nya akhirnya penyusun dapat menyelesaikan makalah Integrasi Numerik; Aturan Simpson 3/8 ini dengan baik dan lancar. Penyelesaian penyusunan makalah Integrasi Numerik; Aturan Simpson 3/8 ini ditujukan untuk memenuhi salah satu tugas yang diberikan oleh Dosen Mata Kuliah Metode Numerik, Program Studi Pendidikan Matematika Fakultas Keguruan dan Ilmu Pendidikan Universitas Pasundan, Bapak Darta, S.Pd., M.Pd. dan Bapak Thesa Kandaga, S.Si., M.Pd. Penyusun menyadari bahwa makalah ini tidak akan berhasil tanpa bantuan berbagai pihak baik secara langsung maupun tidak langsung. Maka dari itu, penyusun mengucapkan banyak terima kasih kepada semua pihak yang telah membantu dalam penyusunan makalah ini. Penyusun menyadari bahwa kerterbatasan kemampuan dan pengetahuan telah menjauhkan makalah ini dari kata sempurna. Untuk itu penyusun sangat mengharapkan kritik dan saran yang bersifat membangun dari pembaca demi perbaikan makalah ini kedepannya. Akhirnya besar harapan penyusun agar makalah ini dapat bermanfaat bagi semua pihak yang bergerak dan terlibat dari dunia pendidikan pada umumnya.
Bandung,
Mei 2020
Penyusun
ii
DAFTAR ISI
KATA PENGANTAR ............................................................................................ ii DAFTAR ISI .......................................................................................................... iii BAB I PENDAHULUAN ....................................................................................... 1 1. 1 Latar Belakang ........................................................................................ 1 1. 2 Rumusan Masalah ................................................................................... 1 1. 3 Tujuan ..................................................................................................... 2 BAB II KAJIAN PUSTAKA .................................................................................. 3 2. 1 Metode Numerik ..................................................................................... 3 2. 2 Angka Signifikan / Bena ........................................................................ 3 2. 3 Deret Taylor ............................................................................................ 6 2. 4 Deret Mc. Laurin .................................................................................... 7 2. 5 Galat (Error) ........................................................................................... 7 2. 6 Interpolasi Polinom NewtonβGregory Maju ........................................ 10 2. 7 Interpolasi Polinom NewtonβGregory Mundur .................................... 12 2. 8 Integrasi Numerik ................................................................................. 14 2. 9 Kaidah Trapesium................................................................................. 15 BAB III PEMBAHASAN ..................................................................................... 18 3. 1 Kaidah Metode Simpson 3/8................................................................. 18 3. 2 Contoh Soal Metode Integrasi Simpson 3/8 ......................................... 22 3. 3 Manfaat Metode Simpson 3/8 .............................................................. 27 BAB IV PENUTUP .............................................................................................. 28 4. 1 Kesimpulan ........................................................................................... 28 DAFTAR PUSTAKA ........................................................................................... 29
iii
BAB I PENDAHULUAN
1. 1 Latar Belakang Integrasi numerik adalah suatu metode yang digunakan untuk mendapatkan nilai β nilai hampiran dari beberapa integral tentu yang memerlukan penyelesaian numerik sebagai hampirannya. Suatu integral dapat diselesaikan dengan 2 cara, yaitu secara analitik dan secara numerik. Banyak permasalahan yang dihadapi dapat dimodelkan ke dalam persamaan integral. Namun persamaan integral ini terkadang memiliki bentuk yang terlalu rumit sehingga sulit untuk diselesaikan dengan menggunakan kaidah β kaidah kalkulus secara analitik. Untuk menangani persamaan yang rumit tersebut dapat menggunakan metode numeriks. Dengan menggunakan metode numerik, solusi eksak dari persoalan yang dihadapi tidak akan diperoleh. Metode numerik hanya bisa memberikan solusi yang mendekati atau menghampiri solusi sejati sehingga solusi numerik dinamakan juga solusi hampiran (approximation solution). Pendekatan solusi ini tentu saja tidak tepat sama dengan solusi sejati, sehingga ada selisih antara keduanya yang disebut dengan galat (error). Semakin kecil galat yang diperoleh berarti semakin dekat solusi hampiran yang diperoleh dengan solusi sejatinya. Tidak seperti metode analitik yang solusinya biasanya disajikan dalam bentuk fungsi matematik yang selanjutnya dapat di evaluasi untuk menghasilkan nilai dalam bentuk angka, sesuai dengan namanya solusi dari metode numerik selalu berbentuk angka.
1. 2 Rumusan Masalah Berdasarkan latar belakang yang telah diuraikan diatas, maka rumusan masalahnya adalah sebagai berikut. 1. Apa itu Metode Simpson 3/8? 2. Bagaimana algoritma dari Metode Simpson 3/8?
1
3. Bagaimana penyelesaian masalah dengan Metode Simpson 3/8? 4. Bagaimana penerapan dari Metode Simpson dalam kehidupan sehari β hari? 1. 3 Tujuan 1. Memahami pengertian dari Metode Simpson 3/8. 2. Mengetahui algoritma dari Metode Simpson 3/8. 3. Memahami langkah β langkah penyelesaian Metode Simpson 3/8. 4. Memahami penerapan dari Metode Simpson dalam kehidupan sehari β hari.
2
BAB II KAJIAN PUSTAKA
2. 1 Metode Numerik Metode numerik adalah satu β satunya metode alternative yang ada dalam upaya menyelesaikan persoalan β persoalan matematis. Metode yang lain dikenal dengan sebutan metode analitik. Ada dua alasan mengapa pilihan dijatuhkan kepada metode numerik. Alasan pertama metode ini memberikan keefesiensinan dan keefektifan di dalam menyelesaikan persoalan β persoalan matematis dikarenakan berkembangnya perangkat keras dan lunak komputer akhir β akhir ini. Alasan yang lain adalah metode numerik memungkinkan untuk mengkaji parametrik dari persoalan dengan medan yang bersifat sembarang. Alasan yang terakhir
ini
lebih
bermakna
ketidakmampuan
metode
analitik
untuk
menyelesaikan persoalan β persoalan matematis aplikasi yang kompleks. Dalam banyak literatur analisa numerik diungkapkan bahwa di dalam metode numerik keputusan menerima atau menolak suatu jawaban aproksimasi berdasarkan kepada toleransi kedekatan yang disepakati. Toleransi yang dibuat menyangkut kesepakatan kesalahan atau galat yang ditimbulkan oleh rumus atau formula yang digunakan. Tentu semakin kecil kesalahan atau galat yang ditimbulkan oleh penggunaan suatu rumus atau formula maka semakin baik hasil aproksimasi yang dihasilkan.
2. 2 Angka Signifikan / Bena Angka signifikan (Bena atau angka penting) adalah bilangan yang diperoleh dari hasil pengukuran yang terdiri dari angka β angka penting yang sudah pasti (terbaca pada alat ukur) dan satu angka terakhir yang ditafsir atau diragukan. Sedangkan angka eksak atau pasti adalah angka yang sudah pasti (tidak diragukan nilainya), yang diperoleh dari kegiatan membilang (menghitung).
3
1. Ketentuan Penulisan Angka Penting a. Semua angka yang bukan nol adalah angka penting. Contoh: 1) 14569 = 5 angka penting. 2) 2546 = 4 angka penting. 3) 6,89 = 3 angka penting. b. Semua angka nol yang berada di antara angka bukan nol termasuk angka penting. Contoh: 1) 2,0067 = 5 angka penting. 2) 7000,2003 = 9 angka penting. 3) 0,005006 = 4 angka penting. c. Semua angka nol yang terletak di belakang angka bukan nol termasuk angka penting. Contoh: 1) 2500, = 4 angka penting. 2) 70000, = 5 angka penting. d. Angka nol yang terletak di belakang angka bukan nol yang terakhir dan di belakang tanda decimal adalah angka penting. Contoh: 23,50000 = 7 angka penting. e. Angka nol yang terletak di belakang angka bukan nol yang terakhir dan tidak dengan tanda desimal adalah dalah angka tidak penting. Contoh: 1) 350000 = 2 angka penting. 2) 141441000 = 6 angka penting. f. Angka nol yang terletak di depan angka bukan nol yang pertama adalah angka tidak penting. Contoh: 0,0000352 = 3 angka penting.
4
2. Kaidah Pembulatan a. Jika angka pertama setelah angka yang hendak dipertahankan adalah 4 angka atau lebih kecil, maka angka itu dan seluruh angka di sebelah kanannya ditiadakan. Contoh: 1) 75,494 = 75,49 (angka 4 yang terakhir ditiadakan). 2) 1,00839 = 1,008 (kedua angka dari paling terakhir ditiadakan). b. Jika angka pertama setelah angka yang akan dipertahankan adalah 5 atau lebih besar, maka angka tersebut dan seluruh angka di bagian kanannya ditiadakan. Angka terakhir yang dipertahankan bertambah satu. 3. Kaidah Penjumlahan dan Pengurangan Apabila melakukan operasi penjumlaha atau pengurangan, maka hasilnya hanya boleh mengandung satu angka taksiran. Angka taksiran adalah angka terakhir dari suatu angka penting. Contoh: Hasil dari 273,219 + 15,5 + 8,43 ? Jawab: Jumlahkan ketiga angka diatas seperti biasanya, diperoleh: 273,219 + 15,5 + 8,43 = 297,149 Selanjutnya bulatkan hasilnya hingga hanya terdapat satu angka taksiran. Dari hasil penjumlahan angka 4 dan 9 ditiadakan, hasilnya menjadi 297,1. 4. Kaidah Perkalian dan Pembagian a. Pada operasi perkalian atau pembagian, hasil yang diperoleh hanya boleh memiliki jumlah angka penting sebanyak bilangan yang angka pentingnya paling sedikit. Contoh: Hitunglah opeasi perkalian dari 0,6283 Γ 2,2 Jawab: Lakukan prosedur perkalian dengan cara biasa, diperoleh: 0,6283 Γ 2,2 = 1,38226 Kemudian bulatkan hasilnya hingga memiliki angka penting sebanyak salah satu bilangan yang memiliki angka penting paling sedikit. 0,6283 = 4 angka penting. 2,2 = 2 angka penting.
5
Jadi yang diambil adalah 2 angka penting. Sehingga 1,38226 hasilnya dibulatkan menjadi 1,4 (dua angka penting). b. Hasil perkalian atau pembagian antara bilangan penting dengan bilangan eksak atau pasti hanya boleh memiliki angka penting sebanya jumlah angka penting pada bilangan penting. Contoh: Hitunglah opeasi perkalian dari 25 Γ 8,95 Jawab: 25 Γ 8,95 = 223,75 Hasilnya dibulatkan menjadi 224 (tiga angka penting) agar sama dengan banyak angka penting pada bilangan penting 8,95.
2. 3 Deret Taylor Deret Taylor adalah deret pangkat π(π§) = π(π§0 ) + ββ π=1
π (π) (π§0 ) π!
(π§ β π§0 )π
yang analitik pada daerah π· = {π§: |π§ β π§0 | < π} (lingkaran dengan pusat di π§0 dan jari β jari π). Maka untuk setiap titik π§ pada lingkaran itu, π(π§) dapat dinyatakan sebagai: π π(π§) = ββ π=0 ππ (π§ β π§0 ) β¦β¦ (|π§ β π§0 | < π)
dengan ππ =
π (π) (π§0 ) π!
β¦β¦ (π = 0, 1, 2, β¦ )
atau dituliskan, π(π§) = π(π§0 ) +
πβ²(π§π ) πβ²β²(π§π ) (π§ β π§0 ) + (π§ β π§0 )2 + β― (|π§ β π§0 | < π) 1! 2! β
= π(π§0 ) + β π=1
π (π) (π§0 ) (π§ β π§0 )π π! 1
Contoh : Tentukan deret Taylor π(π§) = π§ di π§ = 1 Jawab
1
1
π π : π(π§) = π§ = 1+(π§β1) = ββ π=0(β1) (π§ β 1) β¦β¦ (|π§ β 1| < 1)
6
2. 4 Deret Mc. Laurin Deret Mc. Laurin merupakan Deret Taylor pada saat π§0 = 0 berbentuk: β
π (π) (0) π π(π§) = β π§ β¦ (|π§| < π) π! π=0
Teorema: π§π
1. π π§ = ββ π=0 π! = 1 + π§ +
π§2
β¦ , |π§| < β
2!
(β1)π
π§3
2π+1 2. sin π§ = ββ =π§β π=0 (2π+1)! π§ (β1)π
2π 3. cos π§ = ββ =1β π=0 (2π)! π§ π§ 2π+1
4. sinh π§ = ββ π=0 (2π+1)! = π§ + π§ 2π
5. cosh π§ = ββ π=0 (2π)! = 1 + 6. 7.
1 1βπ§ 1 1+π§
π§3
π§2 2!
π§2
3!
2!
+
+
+
3!
4!
π§4
π§5
π§4
+
5!
4!
π§5 5!
β¦ , |π§| < β
β¦ , |π§| < β
β¦ , |π§| < β
β¦ , |π§| < β
π 2 3 = ββ π=0 π§ = 1 + π§ + π§ + π§ + β― , |π§| < 1 π π 2 3 = ββ π=0(β1) π§ = 1 β π§ + π§ β π§ + β― , |π§| < 1 1
Contoh : Perderetan π(π§) = 1+π§ dalam deret Mc. Laurin. Jawab
1
: Fungsi π(π§) = 1+π§ tidak analitik di π§ = β1, sehingga daerah keanalitikan |π§| < 1 β
β
π=0
π=0
1 1 π(π§) = = = β(βπ§)π = β(β1)π π§ π β¦ |π§| < 1 1 + π§ 1 β (βπ§)
2. 5 Galat (Error) Data numerik adalah suatu aproksimasi (taksiran) yang sesuai sampai dengan dua, tiga, atau lebih tempat desimal. Kadang metode yang digunakan pun adalah suatu aproksimasi. Oleh sebab itu galat dalam hasil perhitungan mungkin disebabkan oleh galat data, atau galat di dalam pemakaian suatu metode, atau kedua β duanya. Dalam bagian ini akan dibicarakan ide dasar tentang galat.
7
1. Tipe Galat a. Galat Inheren (Inheren Error) Galat inheren merupakan galat bawaan akibat penggunaan suatu metode numerik. Akibat perhitungan numerik sebagian besar adalah tidak eksak, dapat menyebabkan data yang diperoleh adalah data aproksimasi. Selain itu, keterbatasan dari alat komputasi seperti tabel matematika, kalkulator atau komputer digital juga membuat perhitungan numerik tidak eksak. Karena keterbatasan tersebut, bilangan β bilangan yang diperoleh adalah hasil pembulatan. Di dalam perhitungan, galat inheren dapat diperkecil melalui pengguaan data yang besar, pemeriksaan galat yang jelas dalam data, dan penggunaan alat komputasi dengan ketelitian yang tinggi. b. Galat pemotongan (Truncation Error) Galat ini disebabkan oleh adanya penghilangan sebarisan suku dari suatu deret atau ekspansi untuk tujuan peringkasan pekerjaan perhitungan. Galat pemotongan adalah galat yang tak dapat dihindarkan. 2. Jenis Galat a. Galat Mutlak (Absolute Error) Galat mutlak adalah selisih numerik antara besar nilai sebenarnya dengan nilai aproksimasinya. Jadi, bila π₯ besar nilai yang sebenarnya, dengan π₯1 nilai pendekatannya (aproksimaksinya), maka galat mutlak πΈπΌ didefinisikan dengan Galat = Nilai sebenarnya β Pendekatan πΈπΌ = π₯ β π₯1 = πΏπ₯ b. Galat Relatif Galat relatif πΈπ
didefinisikan dengan πΊππππ‘
Galat relatif = πππππ π ππππππππ¦π πΈπ
=
πΈπΌ πΏπ₯ = π₯ π₯
8
Kemudian persentase galat dihitung dari galat relatif yang diberikan dalam bentuk, Persentase πΈπ
=
πΊππππ‘ π ππππππππ¦π πππππ π ππππππππ¦π
Γ 100%
ππ
= 100πΈπ
c. Galat Global Misal π’ = π(π₯1 , π₯2 , β¦ , π₯π ) adalah fungsi dengan variabel banyak π₯1 = (1, 2, β¦ , π), dam misalkan galat dari tiap π₯1 adalah βπ₯1 . Galat βπ’ dari π’ diberikan dalam bentuk, π’ + βπ’ = π(π₯1 + βπ₯1 , π₯2 + βπ₯2 , β¦ , π₯π + βπ₯π ) Perluasan ruas kanan dari galat global tersebut oleh deret Taylor menghasilkan π
π’ + βπ’ = π(π₯1 , π₯2 , β¦ , π₯π ) + β π=1
ππ βπ₯ ππ₯π π
+ semua suku yang memuat (βπ₯π )2 + semua suku yang lain Anggap bahwa galat dalam π₯π adalah kecil dan
βπ₯π π₯π
< 1. Kemudian semua
suku setelah suku kedua pada ruas kanan persamaan diatas diabaikan. Persamaan menjadi, π
βπ’ β β π=1
ππ ππ ππ ππ βπ₯π = βπ₯π + βπ₯2 + β― + βπ₯ ππ₯π ππ₯π ππ₯2 ππ₯π π
Bila diperhatikan formula (1.12) bentuknya sama dengan diferensial total dari π’. Formula untuk galat relatif adalah sebagai berikut, πΈπ
=
βπ’ ππ’ βπ₯π ππ’ βπ₯2 ππ’ βπ₯π = β + β + β―+ β π’ ππ₯π π’ ππ₯2 π’ ππ₯π π’
d. Galat Dalam Aproksimaksi Deret Galat yang ada dalam aproksimaksi suatu deret dapat di evaluasi oleh sisa sesudah suku β suku ke π. Pandang deret Taylor untuk π(π₯) pada π₯ = π yang diberikan dalam bentuk, π(π₯) = π(π) + (π₯ β π)π β²(π) +
(π₯ β π)2 β²β²(π) (π₯ β π)πβ1 πβ1 (π) + π
π (π₯) π + π (π β 1)! 2!
9
Suku terakhir dalam deret diatas dikenal dengan sebutan suku sisa deret Taylor yang didefinisikan sebagai berikut, π
π (π₯) =
(π₯ β π)π π π (π), π < πΌ < π₯ π!
Untuk suatu barisan yang konvergen, suku β suku sisa akan mendekati nol untuk π β β. Jadi, bila kita mengaproksimasi π(π₯) oleh π suku pertama dari deret tersebut maka galat maksimum yang dibuat dalam aproksimasi tersebut diberikan oleh suku sisa.
2. 6 Interpolasi Polinom NewtonβGregory Maju Polinom NewtonβGregory Maju diturunkan dari tabel selisih maju. Penurunan rumus polinom NewtonβGregory Maju dikembangkan berdasarkan pada tabel selisih maju. 1. Penurunan Rumus Polinom NewtonβGregory Maju π[π₯1 , π₯2 ] =
π(π₯1 ) β π(π₯0 ) π₯1 β π₯0
βπ(π₯0 ) β βπ0 = 1! β π[π₯2 , π₯1 ] β π[π₯1 , π₯0 ] π[π₯1 , π₯2 , π₯0 ] = π₯2 β π₯0 =
π(π₯2 ) β π(π₯2 ) π(π₯1 ) β π(π₯0 ) β π₯2 β π₯1 π₯1 = π₯2 β π₯0 βπ1 β βπ0 β = 2β =
β2 π0 β2 π0
=
β2 π0 2! β2
Bentuk Umum: π[π₯π , β¦ , π₯1 , π₯0 ] =
βπ π(π₯0 ) π!βπ
βπ π
= π!βπ0
10
dengan demikian polinom Newton untuk data berjara sama dapat ditulis sebagai, ππ (π₯) = π(π₯) + (π₯ β π₯0 )π[π₯1 , π₯2 ] + (π₯ β π₯0 )(π₯ β π₯1 )π(π₯2 , π₯1 , π₯0 ) + β― + (π₯ β π₯0 )(π₯ β π₯1 ) β¦ (π₯ β π₯πβ1 )π[π₯π , π₯πβ1 , β¦ , π₯1 , π₯0 ] βπ0 β2 π0 + (π₯ β π₯0 )(π₯ β π₯1 ) +β― 1! β 2! β2 βπ π0 (π₯ )(π₯ ) (π₯ ) + β π₯0 β π₯1 β¦ β π₯πβ1 π! βπ
= π0 + (π₯ β π₯0 )
Persamaan ini dinamakan polinom NewtonβGregory maju. Persmaan diatas juga dapat ditulis sebagai relasi rekursif: ππ (π₯) = ππβ1 (π₯) + (π₯ β π₯0 )(π₯ β π₯1 ) β¦ (π₯ β π₯πβ1 )
βπ π0 π! βπ
Jika titik β titik berjarak sama dinyatakan sebagai, π₯π = π₯0 + πβ, π = 0, 1, 2, β¦ , π dan nilai π₯ yang diinterpolasikan adalah π₯ = π₯0 + π β, π β π
maka persamaan polinom NewtonβGregory Maju dapat juga ditulis dalam parameter π sebagai, ππ (π₯) = π0 +
π β π (π β 1)β2 2 π (π β 1)(π β 2) β¦ (π β π + 1)βπ π βπ0 + β π + β― + β π0 0 1! β 2! β2 π! βπ
yang menghasilkan ππ (π₯) = π0 +
π π (π β 1) 2 π (π β 1)(π β 2) β¦ (π β π + 1) π βπ + β π0 + β― + β π0 1! 0 2! π!
2. Algoritma Polinom Interpolasi Maju a. Definisikan fungsi π(π₯) b. Tentukan selang π(π₯) c. Tentukan jarak antar selang atau β d. Tentukan derajat π e. Buatlah tabel selisih maju
11
f. Tentukan π π =
π₯ β π₯0 β
g. Cari ππ (π₯) = π0 +
π π (π β 1) 2 π (π β 1)(π β 2) 3 βπ0 + β π0 + β π0 + β― 1! 2! 2! π (π β 1)(π β 2) β¦ (π β π + 1) π + β π0 π!
2. 7 Interpolasi Polinom NewtonβGregory Mundur Polinom NewtonβGregory Mundur (NewtonβGregory Backward) dibentuk dari tabel selisih mundur. Polinom ini sering digunakan pada perhitungan nilai turunan (derivatif) secara numerik. Titik-titik yang digunakan berjarak sama, yaitu π₯0 , π₯β1 , π₯β2 , β¦ , π₯βπ yang dalam hal ini, π₯π = π₯0 + πβ dengan π = 0, β1, β2, β¦ , βπ dan nilai π₯ yang diinterpolasikan adalah π₯ = π₯0 β π β dengan π β β 1. Penurunan Rumus Interpolasi Newton Gregory Mundur π(π₯0 , π₯1 ) =
π(π₯0 ) β π(π₯1 ) π₯0 β π₯1
βπ0 β βπ0 = 1! β =
π(π₯0 ) β π(π₯β1 ) π(π₯β1 ) β π(π₯β2 ) β π₯0 β π₯β1 π₯β1 β π₯β2 π(π₯0 , π₯1 , π₯2 ) = π₯0 β π₯β2 βπ0 βπ1 β β = β 2β =
β2 π0 2! β2
Bentuk umum: π[π₯0 , π₯β1 , π₯β2 , β― , π₯βπ ] =
βn π(π₯0) π!βπ
βn π
= π!βπ0 , π = 0,1,2, β¦
12
Selanjutnya, π(π₯) β ππ (π₯) ππ (π₯) = π(π₯0 ) + (π₯ β π₯0 )π[π₯0 , π₯1 ] + (π₯ β π₯0 )(π₯0 β π₯1 )π[π₯0 , π₯1 , π₯2 ] + β¦ + (π₯ β π₯0 )(π₯0 β π₯1 ) β¦ (π₯ β π₯π+1 )π[π₯π , π₯π+1 , β¦ , π₯1 , π₯0 ] β2 π
βπ
= π(π₯0 ) + (π₯ β π₯0 ) 1!β0 + (π₯ β π₯0 )(π₯0 β π₯1 ) 2!β20 + β¦ + (π₯ β π₯0 )(π₯0 β π₯1 ) β¦ (π₯ β π₯π+1 ) = π(π₯0 ) + π
βπ0 1!
+ π (π + 1)
β2 π0 2!
βn π0 π!βπ
+ β¦ + π (π + 1) β¦ (π + π β 1)
βn π0 π!
Polinom Newton-Gregory mundur yang menginterpolasi (π + 1) titik data adalah sebagai berikut: π
π(π₯) β ππ (π₯) = β ( π=0
π +πβ1 π ) β π0 π
π βπ0 π (π + 2 β 1)β2 π0 π (π + 2 β 1)(π + 3 β 1)β3 π0 ππ (π₯) = π0 + + + +β― 1! 2! 3! π (π + 2 β 1)(π + 3 β 1) β¦ (π + π β 1)βπ π0 + π! = π0 + π βπ0 +
π (π + 1)β2 π0 π (π + 1)(π + 2)β3 π0 + +β― 2! 3! π (π + 1)(π + 2) β¦ (π + π β 1)βπ π0 + π!
2. Algoritma Interpolasi Newton Gregory Mundur a. Tentukan fungsi π(π₯) b. Tentukan selang (π₯π ) c. Tentukan jarak antar selang atau β d. Tentukan derajat π e. Buatlah tabel selisih mundur f. Tentukan π π =
π₯ β π₯0 β
13
g. Cari ππ (π₯) = π0 +
π π (π + 2 β 1) 2 π (π + 2 β 1)(π + 3 β 1) 3 βπ0 + β π0 + β π0 1! 2! 3! π (π + 2 β 1)(π + 3 β 1) β¦ (π β π + 1) π + β―+ β π0 π!
2. 8 Integrasi Numerik Integrasi numerik adalah proses mencari hampiran luas bidang yang dibatasi oleh π(π₯) dan sumbu π₯ pada selang tertutup [π, π]. Jika π(π₯) dihampiri dengan polinomial ππ (π₯), maka integrasi numerik ditulis dalam bentuk, π
πΌ = β« π(π₯)ππ₯ π π
β β« ππ (π₯)ππ₯ π
Proses pencarian nilai hampiran πΌ dilakukan jika: 1. Fungsi π(π₯) disebut integran, mempunyai bentuk yang sulit untuk dilakukan proses integrasi. 2. Nilai π₯ dan π(π₯) hanya dalam bentuk tabel diskrit.
Gambar 2.1
Gambar 2.2
Luas bidang yang dibatasi π(π₯)
Hampiran luas bidang yang dibatasi ππ (π₯)
Proses menentukan nilai hampiran integrasi numerik dilakukan dengan beberapa cara atau metode, yaitu kaidah trapesium, kaidah titik tengah, kaidah Simpson, serta Kuadratur Gauss.
14
2. 9 Kaidah Trapesium Kaidah trapesium merupakan kaidah integrasi numerik yang didasarkan pada penjumlahan segmen-segmen berbentuk trapesium.
Gambar 2.3 Luas satu trapesium/pias π₯1
β« π(π₯)ππ₯ = π₯0
β [π(π₯0 ) + π(π₯0 )] 2
Persamaan diatas dikenal dengan nama kaidah trapesium.
Gambar 2.4 Luas beberapa n buah trapesium/pias
15
π
π₯1
π₯2
π₯π
β« π(π₯)ππ₯ β β« π(π₯)ππ₯ + β« π(π₯) + β― + β« π(π₯) π
π₯0
π₯1
π₯πβ1
β β β [π(π₯0 ) + π(π₯1 )] + [π(π₯1 ) + π(π₯2 )]+. . . + [π(π₯πβ1 ) + π(π₯π )] 2 2 2 β β [π(π₯0 ) + 2π(π₯1 ) + 2π(π₯2 )+. . . +2π(π₯πβ1 ) + π(π₯π )] 2 β
πβ1
β β (π0 + 2 β ππ + ππ ) 2 π=1
dengan ππ = π(π₯π ), π = 0, 1, 2, β¦ , π Gala
1
: πΈ β β 12 β3 π β²β²(π‘) , 0 < π‘ < β β3
Galat total : πΈπ‘ππ‘ β β 12 (π β²β² 0 + π β²β²1 + π β²β² 2 + β― + π β²β² πβ1 ) β3
β βπ 12 π β²β²(π‘) , π < π‘ < π π
β
πβ1
β« π(π₯) = (π0 + 2 β ππ + ππ ) + ()(β2 ) π
2
π=1
Algoritma Kaidah Trapesium 1. Mendefinisikan fungsi yang akan diintegrasikan π¦ = π(π₯) 2. Menentukan batas bawah (π) dan batas atas (π) integrasi 3. Menentukan jumlah segmen atau pias π 4. Menghitung lebar segmen yaitu β =
πβπ π
5. Buatlah tabel kaidah trapesium 6. Menentukan nilai integrasi menggunakan kaidah trapesium πβ1
β πΏ = [π(π₯0 ) + 2 β π(π₯π ) + π (π₯π )] 2 π=1
7. Menentukan nilai integrasi sejatinya 8. Menentukan galat kaidah trapesium β3
πΈ = β 12 (π β π)πβ²β²(π‘), π < π‘ < π
16
9. Menentukan nilai sejati (terletak diantara batas galat minimum dan maksimum) Nilai integrasi menggunakan kaidah trapesium β batas galat maksimum, dan nilai integrasi menggunakan kaidah trapesium β batas galat minimum. π
10. Menentukan galat hasil integrasi β«π π(π₯)ππ₯
17
BAB III PEMBAHASAN
Integral numerik merupakan suatu cara untuk menghitung aproksimasi luas daerah dibawah fungsi yang dimaksud pada selang yang diberikan. Aturan simpson adalah metode untuk integrasi numerik, pendekatan integral tertentu. Jika pada metode Simpson 1/3 digunakan pendekatan polinomial berderajat 2 untuk mencari luas dibawah kurva, pada metode Simpson 3/8 digunakan pendekatan polinomial berderajat 3 untuk memperoleh hasil yang lebih baik. Metode simpson 1/3 biasanya lebih disukai karena mencapai ketelitian order tiga dan hanya memerlukan tiga titik, dibandingkan metode simpson 3/8 yang membutuhkan empat titik. Dalam pemakaian banyak pias, metode simpson 1/3 hanya berlaku untuk jumlah pias genap. Apabila dikehendaki jumlah pias ganjil, maka dapat digunakan metode trapesium. Tetapi metode ini tidak begitu baik karena adanya kesalahan yang cukup besar. Untuk itu kedua metode dapat digabung, yaitu sejumlah genap pias digunakan metode simpson 1/3 sedang tiga pias sisanya digunakan metode simpson 3/8.
3. 1 Kaidah Metode Simpson 3/8 Seperti halnya pada kaidah Simpson 1/3, hampiran nilai integrasi yang lebih teliti dapat ditingkatkan terus dengan mengunakan polinom interpolasi berderajat lebih tinggi pula. Misalkan sekarang fungsi π(π₯) kita hampiri dengan polinom interpolasi derajat 3. Luas daerah yang dihitung sebagai hampiran nilai integrasi adalah daerah di bawah kurva polinom derajat 3 tersebut parabola (gambar). Untuk membentuk polinom interpolasi derajat 3, dibutuhkan 4 buah titik data, misalkan titik-titk tersebut (0, π(0)), (β, π(β)), (2β, π(2β)), πππ (3β, π(3β)).
18
Gambar 3.1 Kaidah Simpson 3/8
Polinom interpolasi NewtonβGregory derajat 3 yang melalui keempat buah titik itu adalah, π₯ π₯(π₯ β β) 2 π₯(π₯ β β)(π₯ β 2β) 3 π3 (π₯) = π(π₯0 ) + βπ(π₯0 ) + β π(π₯0 ) + β π(π₯0 ) 2 β 2! β 3! β3 π₯ π₯(π₯ β β) 2 π₯(π₯ β β)(π₯ β 2β) 3 = π0 + βπ0 + β π + β π0 β¦ (1) 0 β 2! β2 3! β3 Integrasi π3 (π₯) di dalam selang [0,3β] adalah, 3β
πΌ=β«
3β
π(π₯)ππ₯ β β« π3 (π₯)ππ₯
0 3β
ββ« 0 3β
ββ« 0
0
π₯ π₯(π₯ β β) 2 π₯(π₯ β β)(π₯ β 2β) 3 [π0 + βπ0 + β π0 + β π0 ] ππ₯ 2 β 2! β 3! β3 π₯ π₯ 2 β π₯β 2 π₯ 3 β 3π₯ 2 β + 2π₯β2 3 [π0 + βπ0 + β π + β π0 ] ππ₯ 0 β 2β2 6β3
β [π₯π0 +
π₯2 π₯3 π₯2β π₯4 3π₯ 3 β π₯ 2 β2 3 3β βπ0 + ( 2 β 2 ) β2 π0 + ( β + ) β π0 ] 2β 6β 4β 24β3 18β3 6β3 0
β 3βπ0 +
9β2 27β3 9β3 81β4 81β4 9β4 3 βπ0 + ( 2 β 2 ) β2 π0 + ( β + ) β π0 β 0 2β 6β 4β 24β3 18β3 6β3
19
9β 27β 9β 2 81β 81β 9β 3 βπ0 + ( β ) β π0 + ( β + ) β π0 2 6 4 24 18 6 9β 27β 2 27β 3 β 3βπ0 + βπ0 + β π0 + β π0 2 12 72 β 3βπ0 +
Mengingat, βπ0 = π1 β π0 β2 π0 = βπ1 β βπ0 = (π2 β π1 ) β (π1 β π0 ) = π2 β 2π1 + π0 dan β3 π0 = β2 π1 β β2 π0 = (βπ2 β βπ1 ) β (π2 β 2π1 + π0 ) = ((π3 β π2 ) β (π2 β π1 )) β (π2 β 2π1 + π0 ) = (π3 β 2π2 + π1 ) β (π2 β 2π1 + π0 ) = π3 β 3π2 + 3π1 β π0 maka selanjutnya, β 3βπ0 +
9β 27β 27β (π1 β π0 ) + (π2 β 2π1 + π0 ) + (π β 3π2 + 3π1 β π0 ) 2 12 72 3
9β 9β 27β 54β 27β 27β 81β 81β 27β π1 β π0 + π2 β π1 + π0 + π3 β π2 + π1 β π 2 2 12 12 12 72 72 72 72 0 9β 27β 27β 9β 54β 81β 27β 81β 27β ) π0 + ( β ) π1 + ( ) π2 + β (3β β + β + β π 2 12 12 2 12 72 12 72 72 3 β 3βπ0 +
β
27β 72
π0 +
81β 72
π1 +
81β 72
π2 +
27β 72
π3
3β 9β 9β 3β π0 + π1 + π2 + π 8 8 8 8 3 3β (π + π1 + π2 + π3 ) β¦ (2) β 8 0 β
20
Sedangkan kaidah Simpson 3/8 gabungan adalah, π
β« π(π₯)ππ₯ β π
3β (π + 3π1 + 3π2 + 2π3 + 3π4 + 3π5 + 2π6 + 3π7 + 3π8 + 2π9 8 0 + β― + 2ππβ3 + 3ππβ2 + 3ππβ1 + ππ ) πβ1
πβ3
π=1 πβ 3,6,9
π=3,6,9
3β (π0 + 3 β ππ + 2 β ππ + ππ ) β¦ (3) β 8 Persamaan (3) ini mudah dihafalkan dengan mengingat pola suku β sukunya, 1, 3, 3, 2,
3, 3, 2,
3, 3, 3, β¦ , 2, 3, 3, 1
Namun penggunaan kaidah Simpson 3/8 mensyaratkan jumlah upselang (π) harus kelipatan 3. Kaidah simpson 3/8 memiliki orde galat yang sama dengan orde galat kaidah simpson 1/3 namun dalam parktek, kaidah simpson 1/3 lebih disukai daripada kaidah simpson 3/8, karena dengan tiga titik (simpson 1/3) sudah diperoleh orde ketelitian yang sama dengan 4 titik (simpson 3/8). Tetapi untuk π kelipatan tiga, kita hanya dapat menggunakan kaidah simpson 3/8, dan bukan simpson 1/3. Algoritma Metode Integrasi Simpson 3/8 1. Masukkan nilai π (jumlah upselang) dan (π₯, π¦) dengan π₯ sebagai titik dan π¦ = π(π₯) 2. Hitung β =
(π₯0 βπ₯πβ1 ) π
dimana π₯0 merupakan titik awal, dan π₯πβ1 merupakan
titik akhir 3. Hitung πΌ = π¦0 β π¦πβ1 4. Jika π kelipatan 3, maka πβ3
β π=3.6,9,β¦
πβ3
= 2 β ππ π=3,6,9β¦
21
5. Dan jika bukan kelipatan 3, maka πβ1
πβ1
β
= 3 β ππ
πβ 3.6,9,β¦
πβ 3,6,9β¦
6. Hitung I π
πβ1
πβ3
πβ 3,6,9,..
π=3,6,9,β¦
3β (π + 3 β ππ + 2 β ππ +ππ ) + 0(β4 ) πΌ = β« π(π₯) = 8 0 π
3. 2 Contoh Soal Metode Integrasi Simpson 3/8 1. Hitunglah 3
β«
πππ (π₯)
2
βπ₯
3
dengan menggunakan kaidah 8 dan up selang yang digunakan adalah π = 7. Penyelesaian : 3
πππ (π₯)
β«
βπ₯
2
π = 2, π = 3, π = 7 β=
πβπ 3β2 = = 0,142857 π 7
Tabel titik di dalam selang [2,3] dengan h = 0,142857 π
π₯π
ππ
0
2
-0,29426
1
2,142857
-0,36982
2
2,285714
-0,43361
3
2,428571
-0,48537
4
2,571429
-0,52496
5
2,714286
-0,5524
6
2,857143
-0,56783
7
3
-0,57157
22
Nilai integrasi π(π₯) didalam selang [2,3] 3
πΌ=β«
πππ (π₯)
2
βπ₯
ππ₯
3β (π + 3π1 + 3π2 +2π3 + 3π4 + 3π5 + 2π6 + 3π7 ) 8 0 3 πΌ = . 0,142857(β2,29426 + 3. β0,36982 + 3. β0,43361 + 2. β0,48537 8 πΌ=
+ 3. β0,52496 + 3. β0,5524 + 2. β0,56783 + 1. β0,57157) πΌ = 0,053571. (β0,975777) πΌ = β0,52274
2. Hitunglah integral dari 1 2
β« πβπ₯ ππ₯ 0 3
dengan menggunakan kaidah 8 dan up selang yang digunakan adalah n = 12 Penyelesaian : 1 2
β« πβπ₯ ππ₯ 0
π = 12 β=
πβπ 1β0 = = 0,08333 π 12
Tabel titik β titik di dalam selang [0,1]dengan h = 0,08333 r
xr
fr
0
0
1
1
0,08333
0,99308
2
0,16666
097261
3
0,24999
0,93942
4
0,33332
0,89485
5
0,41665
0,84064
6
0,49998
0,77882
7
0,58331
0,71159
23
8
0,66664
0,6412
9
0,74997
0,56981
10
0,8333
0,49938
11
0,91663
0,43162
12
1
0,36788
Nilai integrasinya π(π₯) di dalam selang [0,1] adalah 1 2
πΌ = β« π βπ₯ ππ₯ 0 1 2 β« πβπ₯ ππ₯
0
πβ1
πβ3
π=1 π β 3,6,9
π=3,6,9
3β = π0 + 3 β ππ + 2 β ππ + ππ 8 [
=
3(0,8333) 8
]
[1 + 3(0,99308) + 3(0,97261) + 3(0,89485) +
3(0,84064) + 3(0,71159) + 3(0,6412) + 3(0,49938) + 3(0,43162) + 2(0,93942) + 2(0,7782) + 2(0,56981 + 0,36788 ] = (0,03124)(23,74889) = 0,74191
3. Hitunglah integral dari 4
πΌ = β« π π₯ ππ₯ 0
dengan menggunakan 5 pias dan βπ₯ = 0,8. Penyelesaian : Metode Simpson 3/8 dengan satu pias a. Integral dihitung dengan persamaan berikut. π(π₯0 ) + 3 π(π₯1 ) + 3 π(π₯2 ) + π( π₯1) πΌ = (π β π) [ ] 8
24
πΌ = (4 β 0) [
π 0 + 3 π 1,3333 + 3 π 2,667 + π 4 ] 8
= 55,07798
b. Besar kesalahan adalah 5,598150 β 55,07798 π= Γ 100% = β2,761 % 5,598150 Metode Simpson 3/8 dengan lima pias π(0) = π 0 = 1 π(0,8) = π 0,8 = 2,22554 π(1,6) = π 1,6 = 4,9530 π(2,4) = π 2,4 = 11,02318 π(3,2) = π 3,2 = 24,53253 = π 4 = 54,59815
π(4) β’
Integral untuk 2 pias pertama dihitung dengan metode simpson 1/3 dengan persamaan berikut : π1 =
πβπ [π(π) + 4π(π) + π(π)] 6
1,6
πΌ = 6 (1 + (4 x 2,22554) + 4,95303) = 3,96138 β’
Tiga pias terakhir digunakan aturan simpson 3/8 πΌ = (π β π) [ πΌ = 2,4 [
π(π₯0 )+3 π(π₯1 )+3 π(π₯2 )+π( π₯1) 8
]
4,95303+(3 π₯ 11,02318)+(3 π₯ 24,53253)+54,59815) 8
]
= 49,86549 β’
Integral total adalah jumlah dari kedua hasil diatas : πΌ = 3,96138 + 49,86549 = 53,82687
β’
Kesalahan terhadap nilai eksak : ππ‘ =
53,598150β53,826873 5,598150
Γ 100% = β0,427%
25
Cara lain 1) Fungsi integrasinya adalah π(π₯) = π π₯ 2) Batas bawah (π) = 0 3) Batas atas (π) = 4 4) Jumlah pias adalah π =
πβπ 4β0 = = 5 β 0,8
5) β = 0,8 6) Tabel aturan simpson I
xi
f(xi)
0
0
1
1
0,8
2,22554
2
1,6
4,9530
3
2,4
11,02318
4
3,2
24,53253
5
4
54,59815
Aturan Simpson 3/8 4
πβ1
πβ3
π=1 π β 3,6,9
π=3,6,9
3β β« π ππ₯ = [π + 3 β ππ + 2 β ππ + ππ ] 8 0 π₯
0
=
3(0,8) [1 + 3(2,22554) + 3(4,9530) + 3(24,53253) 8 + 3(54,59815) + 2(11,02318) + 54,59815 ]
= (0,3)(1 + 6,67662 + 14,85909 + 73,59759 + 163,79445 + 22,04634 + 54,59815) = 100,971
26
7) Nilai integrasi sejatinya 4
β« ππ₯ ππ₯ = ππ₯ ]40 0 4
= π π₯ ]0 = 54,59815 β 1 = 53,59815
3. 3 Manfaat Metode Simpson 3/8 1. Aturan simpson adalah seperangkat aturan yang digunakan dalam stabilitas kapal dan arsitektur angkatan laut, untuk menghitung luas dan volume angka yang tidak teratur. 2. Aturan simpson digunakan untuk menghitung volume sekoci dan oleh surveyor untuk menghitung volume lumpur tangka minyak kapal. 3. Aturan ke-3 simpson digunakan untuk menemukAan volume antara dua koordinat. 4. Aturan simpson digunakan petugas kapal untuk memeriksa bahwaarea dibawah kurva GZ kapal memenuhi kriteria stabilitas IMO. 5. Aturan 3/8 adalah tentang dua kali lebih akurat dari metode standar,tetapi menggunakan satu nilai fungsi lagi.
27
BAB IV PENUTUP
4. 1 Kesimpulan Metode Integrasi Simpson 3/8 adalah hampiran nilai integrasi yang lebih teliti dapat ditingkatkan terus dengan mengunakan polinom 3 pada empat titik data diskrit yang mempunyai jarak yang sama. Algoritma Metode Integrasi Simpson 3/8 1. Masukkan nilai π (jumlah upselang) dan (π₯, π¦) dengan π₯ sebagai titik dan π¦ = π(π₯) 2. Hitung β =
(π₯0 βπ₯πβ1 ) π
dimana π₯0 merupakan titik awal, dan π₯πβ1 merupakan
titik akhir 3. Hitung πΌ = π¦0 β π¦πβ1 4. Jika π kelipatan 3, maka πβ3
β
πβ3
= 2 β ππ
π=3.6,9,β¦
π=3,6,9β¦
5. Dan jika bukan kelipatan 3, maka πβ1
β
πβ1
= 3 β ππ
πβ 3.6,9,β¦
πβ 3,6,9β¦
6. Hitung I π
πβ1
πβ3
πβ 3,6,9,..
π=3,6,9,β¦
3β (π + 3 β ππ + 2 β ππ +ππ ) + 0(β4 ) πΌ = β« π(π₯) = 8 0 π
28
DAFTAR PUSTAKA https://en.wikipedia.org/wiki/Simpson%27s_rule https://www.academia.edu/35588482/METODE_SIMPSON_1_3_DAN_3_8 https://www.mathworks.com/matlabcentral/fileexchange/72562-simpson-s-3-8rule-composite https://aimprof08.wordpress.com/2012/09/07/aturan-simpson-3-per-8-simpsonrule/
29
1