METODE ITERASI DUA LANGKAH BERDASARKAN ATURAN KUADRATUR BARU UNTUK MENYELESAIKAN PERSAMAAN NONLINEAR Meutia Raya Fitri1∗ 1
Mahasiswa Program Studi S1 Matematika Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Riau Kampus Binawidya Pekanbaru (28293), Indonesia ∗
[email protected]
ABSTRACT This article discusses the two-step iterative method derived following the pattern in obtaining the iterative method proposed by Ujevic [Appl. Math. Comput. 174 (2006): 1416-1426.]. Analytically using fixed point iteration, it is shown that the proposed method converges to the root of nonlinear equations. Then numerical simulations show that the proposed method is comparable to methods Ujevic’s method and better than Newton’s method. Keywords: Nonlinear equations, Newton method, Ujevic’s method ABSTRAK Artikel ini membahas tentang metode iterasi dua langkah yang diturunkan mengikuti pola penurunan metode iterasi yang dikemukakan oleh Ujevic [Appl. Math. Comput. 174(2006): 1416-1426.]. Secara analitik menggunakan iterasi titik tetap ditunjukkan bahwa metode yang dikemukakan konvergen ke akar persamaan nonlinear. Kemudian melalui simulasi numerik terlihat bahwa metode iterasi yang dikemukakan adalah sebanding dengan metode Ujevic dan lebih baik dari metode Newton. Kata kunci: Persamaan nonlinear, metode Newton, metode Ujevic 1. PENDAHULUAN Teknik mendapat metode untuk menyelesaikan persamaan nonlinear secara numerik berkembang sangat pesat. Bermula dari memodifikasi metode Newton yang bentuk iterasinya f (xk ) , k = 0, 1, 2, · · · , f ′ (xk ) ̸= 0, xk+1 = xk − ′ f (xk ) berlanjut dengan mengkombinasikan metode Newton dengan metode lainnya. Pada artikel ini diturunkan metode iterasi baru yang penurunnnya mengikuti penurunan metode Ujevic[6] yang merupakan review dari karya Zhu et al.[5]. Pembahasan dimulai dengan menurunkan metode iterasi dan kekonvergenannya ke akar Repository FMIPA
1
persamaan nonlinear dibagian dua. Dibagian tiga diberikan simulasi numerik dengan melakukan perbandingan dengan metode Ujevic dan metode Newton. 2. METODE ITERASI DUA LANGKAH BERDASARKAN ATURAN KUADRATUR BARU UNTUK MENYELESAIKAN PERSAMAAN NONLINEAR 2.1 Penurunan Metode Iterasi Dua Langkah Pandang fungsi K1 (x, t) yang didefinisikan dengan t − 2a + b , t ∈ [a, x], 3 K1 (x, t) = a + 2b t − , t ∈ [x, b]. 3 Selanjutnya bila dikalikan f ′ (t) pada kedua ruas persamaan (1), diperoleh ( ) 2a + b ′ f (t), t ∈ [a, x], t− 3 ) ′ ( K1 (x, t)f (t) = a + 2b ′ f (t), t ∈ [x, b], t− 3
(1)
(2)
dimana x ∈ [a, b]. Dengan mengintegralkan K1 (x, t)f ′ (t) persamaan (2) pada interval [a, b] dan menggunakan integral parsial [4], setelah penyederhanaan, diperoleh ) ) ∫ b ∫ x( ∫ b( 2a + b ′ a + 2b ′ ′ K1 (x, t)f (t)dt = f (t)dt + f (t)dt t− t− 3 3 a a x [( ) ]x ∫ x 2a + b f (t)dt = t− f (t) − 3 a a [( ) ]b ∫ b a + 2b + t− f (t) − f (t)dt 3 x x ( ) ∫ b ∫ b b−a ′ K1 (x, t)f (t)dt = f (a) + f (b) + f (x) − f (t)dt. (3) 3 a a a+b , kemudian substitusikan x ke persamaan (3) maka Jika dimisalkan x = 2 diperoleh ) ( ( )) ∫ b ∫ b ( a+b b−a a+b ′ f (a) + f (b) + f − f (t)dt. (4) K1 , t f (t)dt = 2 3 2 a a ( ) ∫ b a + b ′ Dari persamaan (4) terlihat bahwa a K1 , t f (t)dt dapat dipandang se2 bagai taksiran error dari ruas kanan persamaan (4). Repository FMIPA
2
Selanjutnya dengan menyusun ulang ruas kanan dari persamaan (4) didapat ) ∫ b ) (( )( ∫ b 2 b−a ′ K1 (x, t)f (t)dt = f (t)dt f (a) + f (b) − 3 2 a a ( ( ) ∫ b ) 1 a+b + (b − a)f − f (t)dt . 3 2 a Dengan menggunakan error aturan trapesium [2] dan aturan titik tengah [2], didapat ( ) ( ) ∫ b 2 (b − a)3 ′′ 1 (b − a)3 ′′ ′ K1 (x, t)f (t)dt = f + f . 3 12 3 24 a Jadi
∫ b ( ) ( ) 3 2 (b − a)3 1 (b − a) ′ ′′ K1 (x, t)f (t)dt ≤ |f | + |f ′′ | 3 12 3 24 a 1 2 ≤ (b − a)3 max|f ′′ | + (b − a)3 max|f ′′ | 72 ∫ b 36 5 K1 (x, t)f ′ (t)dt ≤ (b − a)3 ||f ′′ ||∞ . 72 a
Kemudian definisikan juga
{ t − a, K2 (x, t) = t − b,
t ∈ [a, x], t ∈ [x, b].
Jika f ′ (t) dikalikan ke kedua ruas persamaan (5) diperoleh { (t − a)f ′ (t), t ∈ [a, x], ′ K2 (x, t)f (t) = (t − b)f ′ (t), t ∈ [x, b], dimana x ∈ [a, b]. Bila diintegralkan persamaan (5) pada [a, b] diperoleh ∫ b ∫ x ∫ b K2 (x, t)dt = (t − a)dt + (t − b)dt a a x [ ]x [ ]b 1 2 1 2 = t − at + t − bt 2 2 x ) ( a ∫ b a+b K2 (x, t)dt = (b − a) x − . 2 a
(5)
(6)
(7)
Bila diintegralkan persamaan (6) secara parsial, diperoleh ) ) ∫ b ∫ x( ∫ b( ′ ′ K2 (x, t)f (t)dt = t − a f (t)dt + t − b f ′ (t)dt a a x ∫ b ∫ x b x f (t)dt f (t)dt + [(t − b)f (t)]x − = [(t − a)f (t)]a − x a ∫ b ∫ b ′ K2 (x, t)f (t)dt = (b − a)f (x) − f (t)dt. (8) a
Repository FMIPA
a
3
∫b 1 ∫b K2 (x, t)dt a f ′ (t)dt menggunakan a b−a persamaan (8), (7), dan Teorema Dasar Kalkulus [3], dapat disederhanakan menjadi ∫ b ∫ b ∫ b 1 ′ K2 (x, t)dt K2 (x, t)f (t)dt − f ′ (t)dt b−a a a ( a ( )) ∫ b ∫ b a+b 1 = (b − a)f (x) − f (t)dt − (b − a) x − f ′ (t)dt b − a 2 a a ( ( )) ∫ b 1 a+b = (b − a)f (x) − f (t)dt − (b − a) x − (f (b) − f (a)) b−a 2 a ∫ b ∫ b ∫ b 1 ′ K2 (x, t)dt f ′ (t)dt K2 (x, t)f (t)dt − b − a a a a ( ) ∫ b a+b = (b − a)f (x) − x − (f (b) − f (a)) − f (t)dt. (9) 2 a
Selanjutnya dihitung
∫b a
K2 (x, t)f ′ (t)dt −
Dari persamaan (9) diperoleh ∫ b ∫ b ∫ b 1 ′′ ′′ K2 (x, t)dt f (t)dt K2 (x, t)f (t)dt − b−a a a a ∫ b (∫ b )∫ b 1 ′′ ′′ K2 (x, s)ds f (t)dt = K2 (x, t)f (t)dt − b−a a a ∫a b ( ) ∫ b 1 K2 (x, s)ds f ′′ (t)dt = K2 (x, t) − b−a a a ∫ b( ) ∫ b 1 ′′ ≤ max|f (t)| K2 (x, s)ds dt K2 (x, t) − b−a a a ∫ b( ) ∫ b 1 = K2 (x, s)ds dt ||f ′′ ||∞ K2 (x, t) − b−a a a 2 (b − a) ||f ′′ ||∞ . ≤ 4 Bila dinyatakan ∫ R1 (x) =
b
K1 (x, t)f ′ (t)dt,
a
maka dari persamaan (3) diperoleh ( ) ∫ b b−a f (t)dt = f (a) + f (x) + f (b) − R1 (x). 3 a
(10)
Bila dinyatakan ∫ R2 = a
Repository FMIPA
b
1 K2 (x, t)f (t)dt − b−a ′
∫
∫
b
K2 (x, t)dt a
b
f ′ (t)dt,
a
4
maka dari persamaan (9) diperoleh ( )( ) ∫ b a+b f (t)dt = (b − a)f (x) − x − f (b) − f (a) − R2 (x). 2 a Berdasarkan persamaan (10), dengan mengabaikan nilai R1 diperoleh ) ( ∫ b b−a f (t)dt = f (a) + f (x) + f (b) , 3 a dengan mengabaikan nilai R2 dari persamaan (11), diperoleh ( )( ) ∫ b a+b f (t)dt = (b − a)f (x) − x − f (b) − f (a) , 2 a dan dengan menyamakan persamaan (12) dan (13) diperoleh ) ( ( )( ) b−a a+b f (a) + f (x) + f (b) = (b − a)f (x) − x − f (b) − f (a) . 3 2 Jika f (b) = 0, pada persamaan (14) diperoleh ( ) ( ) b−a a+b f (a) + f (x) = (b − a)f (x) + x − f (a). 3 2 Jika b = ¯b adalah solusi dari persamaan (15) maka diperoleh ( ¯b − a ¯b ) a + (f (a) + f (x)) = (¯b − a)f (x) + x − f (a). 3 2
(11)
(12)
(13)
(14)
(15)
(16)
Dengan melakukan manipulasi aljabar pada persamaan (16) dan menyelesaikan untuk ¯b diperoleh ¯b = a + 6(x − a)
f (a) . 5f (a) − 4f (x)
Misalkan x=a−α
f (a) , f ′ (a)
0 < α ≤ 1,
dan dengan mengambil ¯b = xk+1 , a = xk dan x = zk persamaan (16) diperoleh xk+1 = xk + 6(zk − xk )
f (xk ) , 5f (xk ) − 4f (zk )
(17)
dengan zk = xk − α
Repository FMIPA
f (xk ) , f ′ (xk )
0 < α ≤ 1.
5
Persamaan (17) inilah yang disebut dengan bentuk umum dari metode iterasi dua langkah untuk persamaan nonlinear. 2.2 Kekonvergenan Metode Iterasi Dua Langkah Misalkan f ∈ C 2 (c, d), f (b) = 0, b ∈ (c, d), f ′ (x) ̸= 0, f ′′ (x) ̸= 0, x ∈ (c, d) didefinisikan sebagai ψ(a) = a + 6(x − a)
f (a) , 5f (a) − 4f (x)
(18)
dengan x=a−α
f (a) , f ′ (a)
0 < α ≤ 1.
(19)
Teorema 1 [5] Misalkan f ∈ C 2 (c, d), f (b) = 0, b ∈ (c, d), f ′ (x) ̸= 0, f ′′ (x) ̸= 0, x ∈ (c, d) dan a0 ∈ U (b) dengan U (b) adalah lingkungan disekitar b sedemikian hingga |ψ ′ (a)| < 1 untuk setiap a ∈ U (b). Kemudian barisan xk+1 = ψ(xk ), dengan x0 = a0 dan ψ definisi dari persamaan (18). Bukti: Berdasarkan persamaan (19) diperoleh ( ) f (a) lim x = lim a − α ′ a→b a→b f (a) f (a) = lim a − lim α ′ a→b a→b f (a) f (a) =b−α lima→b f ′ (a) lim x = b. a→b
′
Perhatikan turunan dari ψ (a) di lingkungan di sekitar b. Jika persamaan (18) diturunkan terhadap a maka diperoleh f (a) ψ ′ (a) = 1 + 6(x′a − 1) 5f (a) − 4f (x) ( ) f ′ (a) f (a)(5f ′ (a) − 4f ′ (x)x′a ) + 6(x − a) − , 5f (a) − 4f (x) (5f (a) − 4f (x))2 dx . da Kemudian persamaan (19) diturunkan terhadap a maka diperoleh ( ′ 2 ) dx f (a) − f (a)f ′′ (a) =1−α . da f ′ (a)2
(20)
dengan x′a =
(21)
Selanjutnya dengan mengambil limit pada persamaan (21), diperoleh ( ) f ′ (a)2 − f (a)f ′′ (a) dx = lim 1 − α lim a→b a→b da f ′ (a)2 ( ) f ′ (a)2 dx f (a)f ′′ (a) lim = lim 1 − lim α lim ′ 2 − lim , a→b a→b a→b f (a) a→b da a→b f ′ (a)2 Repository FMIPA
6
karena lima→b f (a) = 0, maka dx =1−α·1−0 a→b da dx = 1 − α. lim a→b da Berdasarkan aturan L’Hospital [1] lim
f (x) f ′ (x) dx = lim ′ a→b f (a) da a→b f (a) f (x) = 1 − α, lim a→b f (a) lim
dan f (x) 1 = lim a→b 3f (a) − 2f (x) a→b 3 − 2 f (x) lim
f (a)
= =
1 3 − lima→b 2 ff (x) (a) 1
3 − 2 lima→b ff (x) (a) 1 = 3 − 2(1 − α) f (x) 1 lim = . a→b 3f (a) − 2f (x) 1 + 2α Kemudian dengan mengambil limit pada persamaan (20) untuk a → b adalah ) ) ( ( f (a) dx ′ −1 lim ψ (a) = 1 + lim 6 a→b a→b da 5f (a) − 4f (x) ( ( )) f (a)(5f ′ (a) − 4f ′ (x) dx f ′ (a) da + lim 6(x − a) − a→b 5f (a) − 4f (x) (5f (a) − 4f (x))2 ′ f (a) = 1 + lim 6(x − 1) lim a→b a→b 5f (a) − 4f (x) ( ) f (a)(5f ′ (a) − 4f ′ (x) dx f ′ (a) da + lim 6(x − a) lim − a→b a→b 5f (a) − 4f (x) (5f (a) − 4f (x))2 ( ) f ′ (x) dx 5−4 ′ 6α α f (a) da + − 6α lim ( =1−6 )2 a→b 1 + 4α 1 + 4α 4f (x) 5− f (a) 1 = 1 − 6α +0 1 + 4α 1 lim ψ ′ (a) = 1 − 6α + 0. a→b 1 + 4α Repository FMIPA
7
1 Jika α ∈ [0, 1] maka |1 − 6α 1+4α | < 1. Dengan demikian dapat disimpulkan bahwa terdapat lingkungan disekitar b sedemikian sehingga
|ψ ′ (a)| < β < 1 untuk setiap a ∈ U (b), dengan menggunakan Teorema Nilai Rata-rata [3] ψ(x) − ψ(y) = ψ ′ (ξ)(x − y), ξ ∈ (x, y), x, y ∈ U (b), atau ψ(x) − ψ(y) = β|(x − y)| untuk 0 < β < 1, x, y ∈ U (b). Oleh karena itu disimpulkan bahwa ψ adalah fungsi kontraktif. Jadi terdapat titik tetap ˆb yaitu ψ(ˆb) = ˆb, sehingga f (ˆb) f (ˆb) ψ(ˆb) = ˆb − 6α · = ˆb, f ′ (ˆb) 5f (ˆb) − 4f (ˆ x) dengan xˆ = ˆb − α
f (ˆb) , f ′ (ˆb)
f (ˆb) = 0,
dengan demikian ˆb = b. 3. KOMPUTASI NUMERIK Pada bagian ini dilakukan perbandingan numerik untuk membandingkan jumlah iterasi dan cots komputasi dari metode Newton, metode Ujevic, dan metode Iterasi Dua Langkah dalam menghampiri nilai akar suatu persamaan. Dalam melakukan perbandingan antara metode Newton, metode Ujevic, dan metode Iterasi Dua Langkah, persamaan nonlinear yang digunakan adalah 1 1. f (x) = − 1 x 1 2. f (x) = − sin(x) + 1 x 3. f (x) = exp(1 − x) − 1 4. f (x) = exp(x2 + 7x − 30) − 1 Hasil dari simulasi numerik terhadap ke empat fungsi di atas, ditunjukkan pada Tabel 1. Pada Tabel 1 dapat dilihat, k menunjukkan jumlah iterasi, x0 menunjukkan tebakan awal, metode Newton disingkat dengan MN, metode Ujevic disingkat dengan MU, dan metode iterasi dua langkah disingkat dengan MIDL, kolom pertama merupakan contoh 1 sampai 4 buah contoh yang berbeda. Pada Tabel 1 juga ditunjukkan akar dan tingkat error dari setiap metode yang dibandingkan. Repository FMIPA
8
Tabel 1: Perbandingan metode Newton, metode Ujevic, dan metode iterasi dua langkah f (x) f1
f2
f3
f4
x0
M etode MN 2.7 NU MIDL MN −1.3 NU MIDL MN 2.9 NU MIDL MN 2.8 NU MIDL
k Akar |f (xk )| |xk − xk−1 | 5 1.0000e + 00 2.3674e + 07 7 1.0000000000000000 4.1467e − 31 7.4357e − 16 7 1.0000000000000000 7.0075e − 29 1.0252e − 14 26 −0.6294464840733333 4.3354e − 36 1.0036e − 18 6 −0.6294464840733333 1.5631e − 20 6.9585e − 11 6 −0.6294464840733333 7.0660e − 31 4.9624e − 16 10 1.0000000000000000 2.4132e − 24 2.1969e − 12 6 1.0000000000000000 2.9764e − 36 2.8173e − 18 6 1.0000000000000000 5.2942e − 25 1.2603e − 12 17 3.0000000000000000 8.2506e − 33 9.8234e − 18 8 3.0000000000000000 1.0360e − 26 1.2710e − 14 8 3.0000000000000000 1.2758e − 28 1.4960e − 15
Dari Tabel 1 maka dapat dilihat bahwa Metode Newton memerlukan iterasi yang lebih banyak dibandingkan dengan metode Ujevic dan metode iterasi dua langkah untuk hampir semua fungsi dari persamaan nonlinear yang diberikan. Hal ini dapat dilihat pada contoh 1 untuk tebakan awal 2.7 metode Newton tidak menemukan akar sampai iterasi 100, sedangkan metode Ujevic dan metode iterasi dua langkah hanya memerlukan 7 iterasi. Pada fungsi kedua untuk tebakan awal −1.3 metode Newton memerlukan 26 iterasi, sedangkan metode Ujevic dan metode iterasi dua langkah hanya memerlukan 6 iterasi. Pada fungsi ketiga untuk tebakan awal 2.9 metode Newton memerlukan 10 iterasi, sedangkan metode Ujevic dan metode iterasi dua langkah hanya memerlukan 6 iterasi. Pada contoh 4 untuk tebakan awal 2.8 metode Newton memerlukan 17 iterasi, sedangkan metode Ujevic dan metode iterasi dua langkah hanya memerlukan 8 iterasi. Ucapan Terimakasih Penulis mengucapkan terimakasih kepada Dr. Imran M., M.Sc. dan Zulkarnain, M.Si. yang telah memberikan arahan dan bimbingan dalam penulisan skripsi yang menjadi acuan artikel ini. DAFTAR PUSTAKA [1] Bartle, R. G. & D. R. Shebert. 2000. Introduction to Real Analysis, 3rd Ed. John Wiley dan Sons, New York. [2] Burden, R. & J. D. Faires. 2003. Numerical Methods, 3th Ed. Brooks Cole, Belmont. [3] Martono, K. 1999. Kalkulus. Penerbit Erlangga. Jakarta. Repository FMIPA
9
[4] Stewart, J. 2012. Calculus, 7th Ed. Brooks Cole, Belmont. [5] Zhu, J., Y. Yao., & P. Guo. 2014. A New Two-Step Iterative Method Based on a New Quadrature Rule For Solving Nonlinear Equations. Journal of Interdisciplinary Mathematics. 17:1–10. [6] Ujevic, N. 2006. A Method for Solving Nonlinear Equations. Appl. Math. Comput, 174: 1416–1426.
Repository FMIPA
10