III PEMBAHASAN 3.1 Analisis Metode Perhatikan persamaan integral Volterra berikut. x
y( x) = f ( x) + λ ∫ K ( x, t , y(t ))dt ,
(8)
atau x ⎡Kn +(x−xn)Jn +(t −tn)Qn ⎤ y(x) = f (x) +λ∫ ⎢ ⎥ dt. +(y− yn)Zn ⎦ a ⎣
a
dengan y(x) merupakan fungsi yang akan ditentukan, a suatu konstanta, f(x) fungsi sembarang yang diketahui dan terdefinisi pada R, dan K(x,t,y(t)) suatu fungsi yang bentuknya tak linear dan bergantung pada variabel x, t, dan fungsi y(t). Berdasarkan bentuk fungsi K(x,t,y(t)), maka berikut ini akan dianalisis dua kasus.
3.1.1 Kasus pertama Fungsi K(x,t,y(t)) terdefinisi pada [a, b] . Selang [ a, b] dibagi menjadi n selang bagian, yaitu [ xi , xi +1 ) , i=0,1,2,3,...,n-1. Uraian deret Taylor fungsi K(x,t,y(t)) di sekitar titik
( xn , tn , yn ) dengan y(xn ) = yn adalah sebagai
Jn =
∂K ( xn , tn , yn ) , ∂x
Qn =
∂K ( xn , t n , y n ) , ∂t
Zn =
∂K ( xn , tn , y n ) . ∂y
K(x, t, y(t)) = K(xn , tn , yn ) ∂K(xn , tn , yn ) ∂x ∂K(xn , tn , yn ) + (t − tn ) ∂t ∂K(xn , tn , yn ) + ( y − yn ) . ∂y + (x − xn )
Jika persamaan (9) disubstitusikan persamaan (8) diperoleh : K(xn , tn , yn ) ⎡ ⎤ ⎢ ⎥ ∂ K ( x , t , y ) n n n ⎥ ⎢+(x − xn ) ⎥ ∂x x⎢ y(x) = f (x) + λ∫ ⎢ ∂K(xn , tn , yn ) ⎥ dt ⎢+(t −tn ) ⎥ a ∂t ⎢ ⎥ ⎢ ∂K(xn , tn , yn ) ⎥ ⎢+( y − yn ) ⎥ ∂y ⎣⎢ ⎦⎥
(12)
Karena dalam integrand pada persamaan (11), t merupakan variabel bebas, y variabel tak bebas, dan x sebuah parameter, maka persamaan (11) dapat ditulis y( x) = f (x) + λ Zn
x
∫ y(t)dt a
x
+ λ [ Kn + ( x − xn )Jn − Zn yn ] ∫ dt
(13)
a
+ λ Qn
berikut:
(11)
dengan K n = K ( x n , t n , y n ) ,
x
∫ ( t − t ) dt. n
a
(9)
Jika integral pada persamaan (13) disederhanakan, maka diperoleh bentuk linear persamaan integral Volterra sebagai berikut. ⎡ K n + ( x − xn ) J n ⎤ y ( x) = f ( x) + λ ⎢ ⎥ ( x − a) ⎣ − Z n yn ⎦ +
ke
λ 2
(14)
Qn ⎡⎣( x − tn ) 2 − (a − tn ) 2 ⎤⎦
+ λ Zn
x
∫ y(t ) dt
; xn ≤ x ≤ xn +1
a
(10)
Kemudian, jika kedua ruas pada persamaan (14) diturunkan terhadap x, maka diperoleh y '( x ) = f '( x ) + λ K n + λ (2 x − x n − a ) J n (15) − λ Z y + λQ (x − x ) n
n
+ λ Z n y( x)
n
n
atau ⎡Kn +(2x− xn −a) Jn ⎤ y'(x) −λ Zn y(x) = f '(x) + λ ⎢ ⎥. ⎣−yn Zn +Qn (x − xn )⎦
(16)
Persamaan (16) merupakan persamaan differensial biasa orde satu dengan penyelesaian analitik dalam bentuk.
Untuk memudahkan dalam perhitungan, misalkan x0 = a , tn = xn dan selang [a, b] dibagi n dengan panjang selang bagian ∆x = xn +1 − xn . Jadi persamaan (19) dapat ditulis : xn+1
yn+1 = yn + fn+1 + λZn exp(λZn xn+1) ∫ f (t)exp(−λZnt)dt
x
y( x) = y( xn ) + f ( x) + λ Z n exp(λ Z n x) ∫ f (t ) exp(−λ Z n t )dt
xn
⎧ ⎡ 2 ⎤ ⎫ ⎪Kn + ⎢∆x + (xn+1 − x0 ) + ⎥ Jn ⎪ λ Zn ⎦ ⎪ 1⎪ ⎣ − ⎨ ⎬ Zn ⎪ 1 ⎪ + (∆x + )Qn ⎪ ⎪ λZn ⎩ ⎭ ⎧ ⎡ ⎫ 2 ⎤ ⎪ ⎢Kn + (xn − x0 + λZ )Jn ⎥ ⎪ ⎪1 ⎢ n ⎥ − f (xn )⎪⎬ + exp(λZn∆x) ⎨ ⎥ 1 ⎪Zn ⎢ ⎪ + Qn ⎢ ⎥ Zn ⎦ ⎩⎪ ⎣ ⎭⎪
xn
⎧ ⎡ 2 ⎤ ⎫ ⎪ K n + ⎢( x − xn ) + ( x − a) + ⎥J ⎪ λ Zn ⎦ n ⎪ 1 ⎪ ⎣ − ⎨ ⎬ Zn ⎪ 1 ⎪ )Qn + ( x − tn + ⎪ ⎪ λ Zn ⎩ ⎭ + exp(λ Z n ( x − xn )) ⎧ ⎪ ⎪1 ⎨ ⎪ Zn ⎩⎪
⎫ 2 ⎡ ⎤ ⎪ ⎢ K n + ( xn − a + λ Z ) J n ⎥ n ⎢ ⎥ − f ( xn ) ⎪⎬ ⎢ ⎥ 1 ⎪ ⎢ + ( xn − tn + )Qn ⎥ Zn ⎣ ⎦ ⎭⎪
(17) Dengan demikian untuk x=xn+1 diperoleh : y(xn+1 ) = y(xn ) + f (xn+1 ) xn+1
+ λZn exp(λZn xn+1 ) ∫ f (t)exp(−λZnt)dt
(20) Jika diasumsikan selang [a, b] dibagi n sama panjang dengan panjang selang bagian ∆xn = h , maka xn = x0 + nh , sehingga persamaan (20) menjadi seperti berikut. xn+1
yn +1 = yn + f n +1 + λ Zn exp(λ Zn xn +1 ) ∫ f (t ) exp(−λ Zn t )dt
xn
⎧ ⎡ 2 ⎤ ⎫ ⎪Kn + ⎢ ( xn+1 − xn ) + ( xn+1 − a) + ⎥J ⎪ λZn ⎦ n ⎪ 1⎪ ⎣ − ⎨ ⎬ Zn ⎪ 1 ⎪ + ( xn+1 − tn + )Qn ⎪ ⎪ λZn ⎩ ⎭ + exp(λZn ( xn+1 − xn ))
xn
⎡ 2 ⎤ 1 ⎪⎧ ⎪⎫ )Qn ⎬ ⎨Kn + ⎢(n + 2)h + ⎥ J n + (h + λ Z λ Z ⎪⎩ n n ⎣ ⎦ ⎭⎪ ⎧ ⎡ ⎫ 2 ⎤ ⎪ ⎢ Kn + (nh + λ Z ) J n ⎥ ⎪ ⎪1 ⎢ ⎪ n ⎥ + exp(λ Zn h) ⎨ − f ( xn ) ⎬ ⎥ 1 ⎪ Zn ⎢ ⎪ + Q ) n ⎥ ⎪⎩ ⎢⎣ ⎪⎭ Zn ⎦ 1 − Zn
⎧ ⎡ 2 ⎫ ⎤ ⎪ ⎢Kn + (xn − a + λZ )Jn ⎥ ⎪ ⎪1 ⎢ ⎪ n ⎥ − f (xn )⎬ ⎨ ⎪ Zn ⎢ + ( x − t + 1 )Q ⎥ ⎪ n n n ⎥ ⎪⎩ ⎢⎣ ⎪⎭ Zn ⎦
(18)
(21) Aplikasi ke dalam perhitungan numerik membutuhkan nilai f(x), λ , K n , J n , Qn , Z n dan h. Hasilnya berupa fungsi y=y(x).
atau xn+1
yn+1 = yn + fn+1 + λZn exp(λZn xn+1) ∫ f (t)exp(−λZnt)dt xn
3.1.2 Kasus kedua Misalkan K(x,t,y(t)) memiliki titik singular
xn , maka formulasi yang diturunkan di atas
⎧ ⎡ 2 ⎤ ⎫ ⎪Kn + ⎢(xn+1 − xn ) + (xn+1 − a) + ⎥ Jn ⎪ λ Zn ⎦ ⎪ 1⎪ ⎣ − ⎨ ⎬ Zn ⎪ 1 ⎪ + (xn+1 −tn + )Qn ⎪ ⎪ λZn ⎩ ⎭
di
+ exp(λZn (xn+1 − xn ))
y(xn ) = yn , sebagai berikut:
1 ⎤ 2 ⎪⎧ 1 ⎡ ⎪⎫ )J + (x −t + )Q ⎥ − f (xn )⎬ ⎨ ⎢Kn + (xn − a + λZn n n n Zn n ⎦ ⎪⎩Zn ⎣ ⎪⎭
dengan
yn = y ( xn ) ,
dan
(19) f n = f ( xn ) .
Persamaan (19) merupakan pendekatan penyelesaian persamaan integral Volterra (8).
tidak berlaku. Dalam kasus ini, misalkan fungsi K(x,t,y(t)) didekati dengan uraian deret Taylor di sekitar titik ( xn +1 , tn +1 , yn ) dengan K ( x, t , y (t )) = K ( xn +1 , tn +1 , yn ) ∂K ( xn +1 , tn +1 , yn ) ∂x ∂K ( xn +1 , tn +1 , yn ) + (t − tn +1 ) (22) ∂t ∂K ( xn +1 , tn +1 , yn ) + ( y − yn ) ∂y + ( x − xn +1 )
Jika persamaan (22) disubstitusi ke dalam persamaan (8), kemudian dilakukan penyederhanaan seperti yang dilakukan pada penurunan persamaan (14), maka diperoleh persamaan berikut.
y ( xn +1 ) = yn + f n +1 + λ Z n +1 exp(λ Z n +1 xn +1 ) xn+1
∫
⎧ ⎫ ⎪ ⎪ K n +1 ⎪ ⎪ 1 ⎪⎪ 2 ⎪⎪ − + ( x − x + ) J ⎨ n +1 0 n ⎬ Z n +1 ⎪ λ Z n +1 +1 ⎪ ⎪ ⎪ 1 Qn +1 + ⎪ ⎪ λ Z n +1 ⎩⎪ ⎭⎪ + exp(λ Z n +1∆xn )
y(x) = f (x) + λ[ Kn+1 + (x − xn+1)Jn+1 − Zn+1 yn ] (x − a)
λ
+ Qn+1 ⎡⎣(x −tn+1)2 − (a −tn+1)2 ⎤⎦ 2
(23)
x
+ λZn+1 ∫ y(t)dt ; xn ≤ x ≤ xn+1
⎧ ⎫ ⎡ ⎤ ⎪ ⎪ ⎢ ⎥ K n +1 ⎪ ⎪ ⎢ ⎥ ⎪⎪ 1 ⎢ ⎪⎪ ⎥ 2 + − ∆ − + − x x x J f x ( ) ( ) ⎢ ⎥ ⎨ n n 0 n +1 n ⎬ λ Z n +1 ⎥ ⎪ Z n +1 ⎢ ⎪ ⎢ ⎥ ⎪ ⎪ 1 + ( − ∆xn )Qn +1 ⎢ ⎥ ⎪ ⎪ Z n ⎣⎢ ⎦⎥ ⎩⎪ ⎭⎪
a
dengan K n +1 = K ( xn +1 , tn +1 , yn ) ,
J n +1 =
∂K ( xn +1 , t n +1 , yn ) , ∂x
Qn +1 =
∂K ( xn +1 , tn +1 , yn ) , ∂t
Z n +1 =
(24)
∂K ( xn +1 , tn +1 , yn ) , ∂y
Persamaan (23) merupakan persamaan differensial biasa orde satu. Jika kedua ruas persamaan (23) diturunkan, maka diperoleh penyelesaian analitik yang serupa dengan persamaan (16) seperti berikut. y( x) = y( xn ) + f ( x) + λ Z n +1 exp(λ Z n +1 x) x
∫ f (t ) exp(−λ Z
t )dt
n +1
xn
⎧ ⎫ ⎡(2 x − xn +1 − a) ⎤ ⎪ ⎢ ⎥J ⎪ K + 2 n +1 ⎢ + ⎥ n +1 ⎪⎪ 1 ⎪⎪ ⎢⎣ λ Z n +1 ⎥⎦ − ⎨ ⎬ Zn +1 ⎪ ⎪ 1 ⎪ + ( x − tn +1 + )Qn +1 ⎪ λ Zn +1 ⎩⎪ ⎭⎪ + exp(λ Z n +1 ( x − xn )) ⎧ ⎫ ⎡ ⎤ ⎪ ⎪ ⎢ ⎥ K n + 1 ⎪ ⎪ ⎢ ⎥ ⎪⎪ 1 ⎢ ⎪⎪ ⎥ 2 ) J n +1 ⎥ − f ( xn ) ⎬ ; ⎢+(2 xn − xn +1 − a + ⎨ λ Z n +1 ⎥ ⎪ Z n +1 ⎢ ⎪ ⎢ ⎥ ⎪ ⎪ 1 )Qn +1 ⎥ ⎢ + ( xn − tn +1 + ⎪ ⎪ Z ⎢ ⎥ n +1 ⎣ ⎦ ⎩⎪ ⎭⎪
(25) Persamaan (25) merupakan pendekatan penyelesaian persamaan integral Volterra (15) dengan yn = y ( xn ) dan f n = f ( xn ) . Untuk memudahkan dalam perhitungan, misalkan x0 = a
,
f (t ) exp(−λ Z n +1t )dt
xn
tn +1 = xn +1 dan selang [a, b]
dibagi n, dengan panjang selang bagian ∆x = xn +1 − xn , diperoleh:
(26) Jika diasumsikan selang [a, b] dibagi n sama panjang dengan panjang selang bagian ∆xn = h , maka xn = x0 + nh , sehingga persamaan (26) menjadi : yn +1 = y n + f n +1 + λ Z n +1 exp(λ Z n +1 xn +1 ) xn +1
∫
f (t ) exp( − λ Z n +1t ) dt
xn
⎧ ⎫ ⎛ ( n + 1) h ⎞ ⎪ ⎜ ⎟J ⎪ K + 2 n +1 ⎜ ⎟ n +1 ⎪⎪ 1 ⎪⎪ ⎜ + λZ ⎟ − ⎨ ⎬ n +1 ⎠ ⎝ Z n +1 ⎪ ⎪ 1 ⎪ ⎪ Qn + 1 + λ Z n +1 ⎩⎪ ⎭⎪ ⎧ ⎪ ⎪ ⎪ ⎪ 1 + exp(λ Z n +1 h ) ⎨ ⎪ Z n +1 ⎪ ⎪ ⎪ ⎩
⎫ ⎡ ⎤ ⎪ ⎢ ⎥ K n +1 ⎪ ⎢ ⎥ ⎪ ⎢ ⎛ ( n − 1) h ⎞ ⎥ ⎪ ⎢ ⎜ ⎥ ⎟ ⎢ + ⎜ + 2 ⎟ J n +1 ⎥ − f ( x n ) ⎬ ⎪ ⎢ ⎜⎝ λ Z n +1 ⎟⎠ ⎥ ⎪ ⎢ ⎥ ⎪ 1 ⎢ ⎥ ⎪ ⎢ + Z Qn + 1 ⎥ n +1 ⎣ ⎦ ⎭
(27) Aplikasi ke dalam perhitungan numerik membutuhkan nilai f(x), λ , K n +1 , J n +1 , Qn +1 , Zn +1 dan h. Hasilnya berupa fungsi y=y(x). Untuk memahami kedua kasus yang diberikan di atas, maka berikut ini dibahas beberapa contoh aplikasi. 3.2 Aplikasi dan Hasil Numerik Dalam bagian ini akan dibahas tiga contoh kasus. Contoh yang kedua dan keempat menjelaskan kasus pertama dari metode ini. Sedangkan contoh kasus yang ketiga menjelaskan kasus kedua dari metode ini.
Contoh 2. Tinjau persamaan integral Volterra berikut pada selang [0,1]. x
1 y( x) = exp( x) − (exp(2x) −1) + ∫ y2 (t )dt 2 0
(28)
Penyelesaian eksak dari persamaan integral (28) adalah (29) y ( x) = exp( x) Berikut ini akan ditentukan penyelesaian hampiran persamaan integral (28) dengan menggunakan metode yang diuraikan di atas. Oleh karena itu, misalkan
y n +1 =
K ( x , t , y ) = y 2 (t ) 1 f ( x ) = exp( x) − (exp(2 x) − 1) 2 λ =1
Dengan demikian persamaan (15) memberikan bentuk linear persamaan integral Volterra dari persamaan (28) sebagai berikut. x
1 y( x) = exp( x) − (exp(2x) −1) − yn2 x + 2 yn ∫ y(t )dt ; 2 0
xn ≤ x ≤ xn +1 Kemudian persamaan persamaan berikut.
(27)
yn y 1 + exp( x n +1 ) − (exp(2 x n +1 ) − 1) + exp(2 y n h )( n − exp( x n ) 2 2 2
1 + exp(2 x n ) − 1)) + 2 y n exp(2 y n x n +1 ) 2
x n +1
∫
xn
(30)
(31) memberikan
(32)
1 (exp( t ) − (exp(2 t ) − 1)) exp( − 2 y n t ) dt 2
Persamaan (32) akan digunakan untuk menentukan hampiran penyelesaian persamaan (28). Dengan menggunakan software MATLAB akan ditentukan penyelesaian y(x) untuk nilai h yang berbeda. Untuk nilai h=0.1, h=0.01 dan h=0.01 dan beberapa nilai x pada [0,1] , diperoleh nilai y(x) seperti diberikan dalam Tabel 1. Tabel 1. Penyelesaian numerik contoh 2 x h=0.1 h=0.01 h=0.001 0.0 1.000000 1.000000 1.000000 0.1 1.104800 1.105200 1.105200 0.2 1.220500 1.221400 1.221400 0.3 1.348100 1.349800 1.349900 0.4 1.488800 1.491800 1.491800 0.5 1.643700 1.648700 1.648700 0.6 1.814000 1.822000 1.822100 0.7 2.000500 2.013600 2.013800 0.8 2.203900 2.225300 2.225500 0.9 2.423300 2.459200 2.459600 1.0 2.655600 2.717600 2.718300 Berikut ini diberikan grafik penyelesaian eksak dan hampiran penyelesaian persamaan (28) dengan nilai h yang berbeda, yakni untuk h=0.1, h=0.01, dan h=0.001.
Gambar 2. Grafik penyelesaian eksak dan hampiran contoh 2 dengan h=0.1
Gambar 3. Grafik penyelesaian eksak dan hampiran contoh 2 dengan h=0.01
Tabel 2 menunjukkan bahwa semakin kecil h, penyelesaian hampiran semakin mendekati penyelesaian eksak, dengan kesalahan mendekati 0. Contoh 3. Tinjau persamaan integral Volterra Hammerstein berikut pada selang [0,1]. x
y( x) = −
x4 5 2 3 1 + x + + ∫ y 2 (t )dt 10 6 8 0 2x
(33)
Titik singular dari persamaan integral (33) adalah x0 = 0 . Penyelesaian eksak dari persamaan integral (33) adalah 1 y ( x) = x 2 + 2 Gambar 4. Grafik penyelesaian eksak dan hampiran contoh 2 dengan h=0.001 Gambar 2, 3, dan 4 memperlihatkan bahwa hampiran penyelesaian eksak akan cukup baik untuk h yang semakin kecil. Kesalahan yang ditimbulkan untuk pemilihan h yang berbedabeda dapat dilihat pada Tabel 2. Tabel 2. Tingkat kesalahan metode linearisasi contoh 2 x h=0.1 h=0.01 h=0.001 0.0 0 0 0 0.1 0.0004 0 0 0.2 0.0009 0 0 0.3 0.0018 0.0001 0 0.4 0.003 0 0 0.5 0.005 0 0 0.6 0.0081 0.0001 0 0.7 0.0133 0.0002 0 0.8 0.0216 0.0002 0 0.9 0.0363 0.0004 0 1.0 0.0627 0.0007 0
(34)
Berikut ini akan ditentukan penyelesaian hampiran dari persamaan integral (33) dengan menggunakan metode yang diuraikan di atas. Oleh karena itu, misalkan 1 2 y (t ) 2x 4 x 5 3 f ( x) = − + x 2 + , 10 6 8 λ =1
K ( x, t , y ) =
(35)
Dengan demikian persamaan (23) memberikan bentuk persamaan integral Volterra linear dari persamaan (33) sebagai berikut. y ( x) = − +
y2 x − x n +1 x4 5 2 3 + x + − nn +1 [1 + ]x 10 6 8 2x x n +1
yn x n +1
x
∫ y (t )dt
xn ≤ x ≤ xn +1
;
0
Kemudian persamaan persamaan berikut.
(27)
(36) memberikan
n +1 yn x4 y t4 5 ( n + 1) h 5 11 y 3 (1 + ) − n +1 + xn2+1 + + nn+1 exp( y n ) ∫ (( − + t 2 + ) exp( − n t )) dt xn +1 xn +1 2 10 6 8 x 10 6 8 xn
x
y ( x) =
+ exp(
yn y (n − 1)h xn4 5 2 11 ) + − xn − ) h)( n (1 − 2 10 6 8 xn +1 xn +1
(37) Persamaan (37) akan digunakan untuk menentukan penyelesaian numerik persamaan (33). Dengan menggunakan software MATLAB akan ditentukan penyelesaian y(x) untuk nilai h=0.1, dan beberapa nilai x pada [0,1]. .
Perhitungan MATLAB untuk penyelesaian eksak, penyelesaian numerik serta tingkat kesalahannya dengan nilai h=0.1 diberikan dalam Tabel 3.
Tabel 3. Penyelesaian numerik, eksak dan tingkat kesalahan contoh 3 dengan h=0.1 Solusi Solusi Tingkat x analitik numerik kesalahan 0.0 1 1 0 0.1 0.51 1.1528 0.6428 0.2 0.54 1.2068 0.6668 0.3 0.59 1.2633 0.6733 0.4 0.66 1.3325 0.6725 0.5 0.75 1.4156 0.6656 0.6 0.86 1.5124 0.6524 0.7 0.99 1.6217 0.6317 0.8 1.14 1.7421 0.6021 0.9 1.31 1.8717 0.5617 1.0 1.5 2.0085 0.5085 Dalam Tabel 4 berikut diberikan tingkat kesalahan untuk beberapa nilai h yang berbeda, yakni untuk h=0.1, h=0.01, dan h=0.001. Tabel 4. Tingkat kesalahan metode linearisasi contoh 3. tingkat kesalahan X h=0.1 h=0.01 h=0.001 0 0 0 0 0.1 0.6428 0.672 0.6719 0.2 0.6668 0.6684 0.6668 0.3 0.6733 0.6605 0.6579 0.4 0.6725 0.6481 0.6446 0.5 0.6656 0.6305 0.626 0.6 0.6524 0.6066 0.6011 0.7 0.6317 0.575 0.5685 0.8 0.6021 0.5342 0.5267 0.9 0.5617 0.4824 0.4738 1 0.5085 0.4176 0.4079 Berdasarkan Tabel 4 diperoleh bahwa semakin kecil nilai h, maka penyelesaian hampiran dengan metode ini semakin mendekati penyelesaian eksaknya. Contoh 4. Tinjau persamaan integral Volterra berikut pada selang [0,1]. x 1 y( x) = exp(− x) + exp(−2 x) + ∫ y2 (t )dt (38) 2 0
Penyelesaian eksak dari persamaan integral (38) adalah (39) y ( x) = exp(− x) Berikut ini akan ditentukan penyelesaian hampiran dari persamaan integral (38) dengan menggunakan metode yang diuraikan di atas. Oleh karena itu, misalkan K ( x, t , y ) = y 2 (t ) 1 f ( x ) = exp( − x ) + exp( −2 x) (40) 2 λ = −1 Dengan demikian persamaan (20) memberikan persamaan berikut. yn 1 + exp(−x(n +1)) + exp(−2x(n +1)) 2 2
yn+1 =
xn+1
+ 2 yn exp(2 yn xn+1 ) ∫ [ exp(−t ) + exp(−2t )] dt xn
1 ⎡y ⎤ + exp(2 yn h) ⎢ n − exp(−xn ) − exp(−2xn )⎥ 2 ⎣2 ⎦
(41) Persamaan (41) akan digunakan untuk menentukan hampiran penyelesaian persamaan (38). Dengan menggunakan software MATLAB akan ditentukan penyelesaian y(x) untuk nilai h yang berbeda. Untuk nilai h=0.1, h=0.01 dan h=0.01 dan beberapa nilai x pada [0,1] diperoleh nilai y(x) seperti diberikan dalam Tabel 5. Tabel 5. Penyelesaian numerik dari contoh 4 dengan h=0.1, h=0.01, dan h=0.001. x
h=0.1
h=0.01
h=0.001
0.0
1.000000
1.000000
1.000000
0.1
0.904500
0.904800
0.904800
0.2
0.818100
0.818700
0.818700
0.3
0.739800
0.740800
0.740800
0.4
0.669000
0.670300
0.670300
0.5
0.604900
0.606500
0.606500
0.6
0.546900
0.548800
0.548800
0.7
0.494300
0.496600
0.496600
0.8
0.446800
0.449300
0.449300
0.9
0.403700
0.406600
0.406600
1.0
0.364700
0.367900
0.367900
Berikut ini diberikan grafik penyelesaian eksak dan hampiran persamaan (38) yakni untuk h=0.1 dan h=0.01. 1
• numerik -- analitik
0.9
Dari Gambar 5 dan 6 terlihat bahwa penyelesaian hampiran mendekati penyelesaian eksak dengan sangat baik. Kesalahan yang ditimbulkan untuk pemilihan h yang berbeda-beda dapat dilihat pada Tabel 6. Tabel 6. Tingkat kesalahan metode linearisasi contoh 4
0.8
0.7
x 0.6
0.5
0.4
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Gambar 5. Grafik penyelesaian eksak dan hampiran contoh 4 dengan h=0.1. 1
• --
0.9
numerik analitik
0.8
0.7
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Tingkat Kesalahan h=0.1 h=0.01 h=0.001 0 0 0 0.0003 0 0 0.0006 0 0 0.001 0 0 0.0013 0 0 0.0016 0 0 0.0019 0 0 0.0023 0 0 0.0025 0 0 0.0029 0 0 0.0032 0 0
Tabel 6 menunjukkan bahwa semakin kecil h, penyelesaian hampiran semakin mendekati penyelesaian eksak, dengan kesalahan mendekati 0.
0.6
0.5
0.4
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Gambar 6. Grafik penyelesaian eksak dan hampiran contoh 4 dengan h=0.01.