The_Weeknd-Starboy-2016-C4 | Total Commander 9.21 Final [PL] [Cracked Sayala! & ZWT] [azjatycki] | Smallfoot (2018) 1080p Telecharger

BAb- 06 Integrasi Numerik

July 10, 2018 | Author: Novita Rusda Yani | Category: N/A
Share Embed Donate


Short Description

matematika...

Description

Bab 6 Integrasi Numerik Pelajarilah jagad raya ini. Jangan kecewa karena dunia tidak mengenal anda, tetapi kecewalah karena anda tidak mengenal dunia. (Kong Fu Tse - filusuf China)

Di dalam kalkulus, integral adalah satu dari dua pokok bahasan yang mendasar disamping turunan (derivative (derivative). ). Dalam kuliah kalkulus integral, anda telah diajarkan cara memperoleh solusi analitik (dan eksak) dari integral Tak-tentu maupun integral Tentu. Integral Tak-tentu dinyatakan sebagai  x) + C   ∫  f ( x)dx = F ( x)

(P.6.1)

Solusinya, F  Solusinya,  F ( x)  x) ,   ,  adalah fungsi menerus sedemikian sehingga F' ( x)  x) = f ( x)  x) , dan  ,  dan C  adalah sebuah konstanta. Integral Tentu menangani perhitungan integral di antara  batas-batas yang telah ditentukan, yang dinyatakan sebagai b

 I =

∫  f ( x)dx  

(P.6.2)

a

Menurut teorema dasar kalkulus integral, persamaan (P.6.2) dihitung sebagai b

∫ 

 x)  f ( x) dx = F ( x)

a

b a

=  F (b) - F  -  F (a)

Secara geometri, integrasi integrasi Tentu sama dengan luas daerah yang dibatasi oleh kurv a  y = f ( x)  x) ,   ,  garis x = a  dan garis x = b  (Gambar 6.1). Daerah yang dimaksud ditunjukkan oleh bagian yang diarsir.

Bab 6 Integrasi Numerik

269

 y

 y = f ( x)  x )

b

a

 x

Gambar 6.1 Tafsiran geometri integral integral Tentu

Fungsi-fungsi yang dapat diintegrasikan dapat dikelompokkan sebagai 1. Fungsi menerus yang sederhana, seperti polinomial, eksponensial, atau fungsi trigonometri. trigo nometri. Misalnya, Misalnya, 2

∫ (6 x

3

− x 2 + cos( x ) − e x )dx

0

Fungsi sederhana seperti ini mudah dihitung integralnya secara eksak dengan menggunakan metode analitik. Metode-metode analitik untuk menghitung integral fungsi yang demikian sudah tersedia, yaitu

∫ axn dx = axn+1/(n /(n+1) + C  ∫ eax dx = eax/a+ C  ∫  sin(ax  sin(ax+ +b) dx = -1/a -1/a  cos(ax cos(ax+ +b) + C  ∫  cos(ax  cos(ax+ +b) dx = dx = 1/a 1/ a  sin(ax sin( ax+ +b) + C  ∫ dx/ dx/ x =  x = ln x + C  ∫  ln xdx = dx = x  x ln x - x  -  x +  + C  2. Fungsi menerus yang rumit, misalnya 2

∫  0

270

(

2 + cos 1 + x

3

2

1 + 0. 5 sin sin  x

)e

0.5 x

dx

Metode Numerik

 y

 y = f ( x)  x )

b

a

 x

Gambar 6.1 Tafsiran geometri integral integral Tentu

Fungsi-fungsi yang dapat diintegrasikan dapat dikelompokkan sebagai 1. Fungsi menerus yang sederhana, seperti polinomial, eksponensial, atau fungsi trigonometri. trigo nometri. Misalnya, Misalnya, 2

∫ (6 x

3

− x 2 + cos( x ) − e x )dx

0

Fungsi sederhana seperti ini mudah dihitung integralnya secara eksak dengan menggunakan metode analitik. Metode-metode analitik untuk menghitung integral fungsi yang demikian sudah tersedia, yaitu

∫ axn dx = axn+1/(n /(n+1) + C  ∫ eax dx = eax/a+ C  ∫  sin(ax  sin(ax+ +b) dx = -1/a -1/a  cos(ax cos(ax+ +b) + C  ∫  cos(ax  cos(ax+ +b) dx = dx = 1/a 1/ a  sin(ax sin( ax+ +b) + C  ∫ dx/ dx/ x =  x = ln x + C  ∫  ln xdx = dx = x  x ln x - x  -  x +  + C  2. Fungsi menerus yang rumit, misalnya 2

∫  0

270

(

2 + cos 1 + x

3

2

1 + 0. 5 sin sin  x

)e

0.5 x

dx

Metode Numerik

Fungsi yang rumit seperti ini jelas sulit, bahkan tidak mungkin, diselesaikan dengan metode-metode integrasi yang sederhana. Karena itu, solusinya hanya dapat dapat dihitung dengan metode numerik. 3. Fungsi yang ditabulasikan, yang dalam hal ini nilai  x dan  x dan f   f ( x)  x) diberikan dalam sejumlah titik diskrit. Fungsi seperti ini sering dijumpai pada data hasil eksperimen di laboratorium atau berupa data pengamatan di lapangan. Pada kasus terakhir ini, umumnya fungsi  f ( x)  x) tidak diketahui secara eksplisit. Yang dapat diukur hanyalah besaran fisisnya saja. Misalnya,  x

f ( x   x )

0.00 0.25 0.50 0.75 1.00

6.0 7.5 8.0 9.0 8.5

Integrasi fungsi seperti ini jelas harus didikerjakan secara numerik.

6.1 Terapan Integral dalam Bidang Sains dan Rekayasa Integral mempunyai banyak terapan dalam bidang sains dan rekayasa. Dalam  praktek  prakt ek reka r eka yasa, seringkali seringkal i fungsi yang diintegrasikan diintegras ikan (integrand  (integrand ) adalah fungsi empirik yang diberikan dalam bentuk tabel, atau integrand -nya -nya tidak dalam  bentuk fungsi elementer (seperti sinh  x,  x, fungsi Gamma Γ (α), dsb), atau fungsi eksplisit f  eksplisit  f  yang  yang terlalu rumit untuk diintegralkan diintegralkan [KRE88 [KRE88]. ]. Oleh sebab sebab itu, metode metode numerik dapat digunakan untuk menghampiri integrasi. Di bawah ini diberikan beberapa contoh persoalan dalam bidang sains dan rekayasa. 1.

Dalam bidang fisika, integral digunakan untuk menghitung persamaan kecepatan. Misalkan kecepatan sebuah partikel merupakan fungsi waktu menerus yang diketahui terhadap waktu, v(t ). ). Jarak total d   yang ditempuh oleh partikel ini selama waktu t  diberikan   diberikan oleh: t 

d  =  =

∫ v(t )dt  0

Bab 6 Integrasi Numerik

271

2.

Dalam bidang teknik elektro/kelistrikan, telah diketahui bahwa harga rata-rata suatu arus listrik yang berosilasi sepanjang satu periode boleh nol. Disamping kenyataan bahwa hasil netto adalah nol, arus terse but mampu menimbulkan kerja dan menghasilkan panas. Karena itu para rekayasawan listrik sering mencirikan arus yang demikian dengan persamaan T 

 I  RMS  =

∫ i

2

(t )dt 

0



yang dalam hal ini  I  RMS   adalah arus  RMS  (root-mean-square), T   adalah  periode, dan i(t ) adalah arus pada rangkaian, misalnya i(t ) = 5e-2t  sin 2πt  =0

3.

untuk 0 ≤ t ≤ T /2 untuk T /2 ≤ t  ≤ T 

Contoh fungsi dalam bentuk tabel adalah pengukuran fluks panas matahari yang diberikan oleh tabel berikut: Waktu, jam

Fluks panas q , kalori/cm/jam

0

0.1

1

1.62

2

5.32

3

6.29

4

7.8

5

8.81

6

8.00

7

8.57

8

8.03

9

7.04

10

6.27

11

5.56

12

3.54

13

1.0

14

0.2

Data yang ditabulasikan pada tabel ini memberikan pengukuran fluks panas q setiap jam pada permukaan sebuah kolektor sinar matahari. Diminta

272

Metode Numerik

memperkiraan panas total yang diserap oleh panel kolektor seluas 150.000 cm2  selama waktu 14 jam. Panel mempunyai kemangkusan penyerapan (absorption), eab, sebesar 45%. Panas total yang diserap diberikan oleh  persamaan t 

 H  = eab

∫ qAdt  0

Demikianlah beberapa contoh terapan integral dalam bidang sains dan rekayasa. Umumnya fungsi yang diintegralkan bentuknya rumit sehingga sukar diselesaikan secara analitik. Karena itu, perhitungan integral secara numerik lebih banyak dipraktekkan oleh para insinyur.

6.2 Persoalan Integrasi Numerik Persoalan integrasi numerik ialah menghitung secara numerik integral Tentu b

 I =

∫  f ( x)dx a

yang dalam hal ini a dan b batas-batas integrasi, f   adalah fungsi yang dapat diberikan secara eksplisit dalam bentuk persamaan ataupun secara empirik dalam bentuk tabel nilai. Terdapat tiga pendekatan dalam menurunkan rumus integrasi numerik. Pendekatan  pertama adalah berdasarkan tafsiran geometri integral Tentu. Daerah integrasi dibagi atas sejumlah pias ( strip) yang berbentuk segiempat. Luas daerah integrasi dihampiri dengan luas seluruh pias. Rumus, dalam bab ini disebut kaidah, integrasi numerik yang diturunkan dengan pendekatan ini digolongkan ke dalam metode pias. Pendekatan kedua adalah berdasarkan polinom interpolasi. Di sini fungsi integrand  f ( x) dihampiri dengan polinom interpolasi  pn( x). Selanjutnya, integrasi dilakukan terhadap  pn( x) karena polinom lebih mudah diintegralkan ketimbang mengintegralkan  f ( x). Rumus integrasi numerik yang diturunkan dengan  pendekatan ini digolongkan ke dalam metode Newton-Cotes , yaitu metode yang umum untuk menurunkan rumus integarsi numerik.. Pendekatan ketiga sama sekali tidak menggunakan titik-titik diskrit sebagaimana  pada kedua pendekatan di atas. Nilai integral diperoleh dengan mengevaluasi nilai fungsi pada sejumlah titik tertentu di dalam selang [-1, 1], mengalikannya

Bab 6 Integrasi Numerik

273

dengan suatu konstanta, kemudian menjumlahkan keseluruhan perhitungan. Pendekatan ketiga ini dinamakan Kuadratur Gauss , yang akan dibahas pada  bagian akhir bab ini.

6.3 Metode Pias Pada umumnya, metode perhitungan integral secara numerik bekerja dengan sejumlah titik diskrit. Karena data yang ditabulasikan sudah berbentuk demikian, maka secara alami ia sesuai dengan kebanyakan metode integrasi numerik. Untuk fungsi menerus, titik-titik diskrit itu diperoleh dengan menggunakan persamaan fungsi yang diberikan untuk menghasilkan tabel nilai. Dihubungkan dengan tafsiran geometri inttegral Tentu, titik-titik pada tabel sama dengan membagi selang integrasi [a, b] menjadi n buah pias ( strip) atau segmen (Gambar 6.2). Lebar tiap pias adalah h =

b − a n

 

(P.6.3)

Titik absis pias dinyatakan sebagai  xr  = a + rh,

r  = 0, 1, 2, ..., n 

(P.6.4)

dan nilai fungsi pada titik absis pias adalah  f r  =  f ( xr )

(P.6.5)

Luas daerah integrasi [a, b] dihampiri sebagai luas n buah pias. Metode integrasi numerik yang berbasis pias ini disebut metode pias. Ada juga buku yang menyebutnya metode kuadratur, karena pias berbentuk segiempat.

r

x r 

f r 

0 1 2 3 4 ... n-2 n-1 n

 x 0  x 1  x 2  x 3  x 4 ...  x n-2  x n-1 xn

 f 0  f 1  f 2  f 3  f 4 ...  f n-2  f n-1 f n

 y

 f n

 f n-1  f 1

 y = f ( x)

 f 2

 f 0

h

h

a = x 0  x 1  x 2

h  x

 xn-1  xn=b

Gambar 6.2 Metode pias

274

Metode Numerik

Kaidah integrasi numerik yang dapat diturunkan dengan metode pias adalah: 1. Kaidah segiempat (rectangle rule) 2. Kaidah trapesium (trapezoidal rule) 3. Kaidah titik tengah (midpoint rule) Dua kaidah pertama pada hakekatnya sama, hanya cara penurunan rumusnya yang berbeda, sedangkan kaidah yang ketiga, kaidah titik tengah, merupakan  bentuk kompromi untuk memperoleh nilai hampiran yang lebih baik.

6.3.1 Kaidah Segiempat Pandang sebuah pias berbentuk empat persegi panjang dari  x = x0 sampai  x = x1  berikut (Gambar 6.3).

 y

 y = f ( x)

h

 x 0

 x1

 x

Gambar 6.3 Kaidah segiempat

Luas satu pias adalah (tinggi pias =  f ( x0) )  x1

∫  f ( x)dx ≈ h f ( x ) 0

(P.6.6)

 x0

atau (tinggi pias = f ( x1) )  x1

∫  f ( x)dx ≈ h f ( x ) 1

(P.6.7)

 x0

Bab 6 Integrasi Numerik

275

Jadi,  x1

∫  f ( x)dx ≈ hf ( x ) 0

 x0  x1

∫  f ( x)dx ≈ hf ( x )

+

1

 x0  x1

2

∫  f ( x)dx ≈ h [ f ( x ) + f ( x )] 0

1

 x0

Bagi setiap ruas persamaan hasil penjumlahan di atas dengan 2, untuk menghasilkan  x1

∫ 

 f ( x) dx ≈

 x0

h 2

 [ f ( x0) + f ( x1)]

(P.6.8)

Persamaan (P.6.8) ini dinamakan kaidah segiempat . Kaidah segiempat untuk satu pias dapat kita perluas untuk menghitung b

 I =

∫  f ( x)dx a

yang dalam hal ini, I   sama dengan luas daerah integrasi dalam selang [a, b]. Luas daerah tersebut diperoleh dengan membagi selang [a, b] menjadi n buah pias segiempat dengan lebar h, yaitu pias dengan absis [ x0 ,  x1], [ x1 ,  x2], [ x2 ,  x3], ... , dan pias [ xn-1 , xn]. Jumlah luas seluruh pias segiempat itu adalah hampiran luas I (Gambar 6.4). Kaidah integrasi yang diperoleh adalah kaidah segiempat gabungan (composite rectangle's rule): b

∫  f ( x)dx ≈ hf ( x ) + hf  ( x ) + hf ( x ) + ... + hf  ( x 0

1

2

n-1)

a b

∫  f ( x)dx ≈ hf  ( x ) + hf  ( x ) + hf  ( x ) + ... + hf  ( x ) 1

2

3

n

+

a

b

2

∫  f ( x)dx ≈

hf ( x0) + 2hf ( x1) + 2hf ( x2) + ... + 2hf ( xn-1) + hf ( xn)

a

276

Metode Numerik

Bagi setiap ruas persamaan hasil penjumlahan di atas dengan 2, untuk menghasilkan b

∫ 

 f ( x) dx ≈

a

h h  f  ( x0) + hf ( x1) + hf ( x2) + ... + hf ( xn-1) +  f  ( xn) 2 2

Jadi, kaidah segiempat gabungan adalah b

∫  f ( x)dx ≈ a

h 2

 ( f 0 + 2 f 1 + 2 f 2+ ... + 2 f n-1 + f n) =

dengan  f r  =  f ( xr ) ,

h 2

n −1

( f 0 + 2

∑ f  + f ) i

n

(P.6.9)

i =1

r  = 0, 1, 2, ..., n .

 y

 y = f ( x)

...

a = x0  x 1  x 2  x3

...  x n-2

Gambar 6.4 Kaidah segiempat

 x n-1  x n = b

 x

gabungan

6.3.2 Kaidah Trapesium Pandang sebuah pias berbentuk trapesium dari  x = x0  sampai x = x1  berikut (Gambar 6.5): Luas satu trapesium adalah  x1

∫ 

 x0

 f ( x) dx ≈

h 2

Bab 6 Integrasi Numerik

[ f ( x0) + f ( x0)]

(P.6.10)

277

Persamaan (P.6.10) ini dikenal dengan nama kaidah trapesium. Catatlah bahwa kaidah trapesium sama dengan kaidah segiempat.

 y

h

 x 0

 x1

 x

Gambar 6.5 Kaidah trapesium

Bila selang [a, b] dibagi atas n buah pias trapesium, kaidah integrasi yang diperoleh adalah kaidah trapesium gabungan  (composite trapezoidal's rule): b

∫  f ( x)dx

 x1

 x2

 x n

 x0

 x1

 xn −1

∫ 

≈  f ( x) dx  +

a



h



h



h

2

2

2

∫  f ( x)dx  + ... + ∫  f ( x)dx

 [ f ( x0) + f ( x1)] +

h 2

 [ f ( x1)+ f ( x2)] + ... +

h 2

 [ f ( xn-1) + f ( xn)]

 [ f ( x0) + 2 f ( x1) + 2 f ( x2) + ... + 2 f ( xn-1) + f ( xn)] n −1

 ( f 0 + 2

∑  f 

1

+ f n)

(P.6.11)

i =1

dengan  f r = f  ( xr ) , r  = 0, 1, 2, ..., n.  

278

Metode Numerik

Program 6.1 Kaidah Trapesium

 procedure  trapesium(a, b : real; n: integer; var  I : real); { Menghitung integrasi f(x) di dalam selang [a, b] dan jumlas pias adalah n dengan menggunakan kaidah trapesium. K.Awal : nilai a, b, dan n sudah terdefinisi K.Akhir: I adalah hampiran integrasi yang dihitung dengan kaidah segi-empat. } var h, x, sigma: real; r : integer;  begin h:=(b-a)/n; {lebar pias} x:=a; {awal selang integrasi}

I:=f(a) + f(b); sigma:=0; for  r:=1 to n-1 do  begin

x:=x+h; sigma:=sigma + 2*f(x); end  ; I:=(I+sigma)*h/2; { nilai integrasi numerik} end  ;

6.3.3 Kaidah Titik Tengah Pandang sebuah pias berbentuk empat persegi panjang dari  x = x0 sampai  x = x1 dan titik tengah absis x = x0 + h/2 (Gambar 6.6). Luas satu pias adalah  x1

∫  f ( x)dx ≈ h f ( x  + h/2) ≈ h f ( x 0

1/2)

P.6.12)

 x0

Persamaan (P.6.12) ini dikenal dengan nama kaidah titik-tengah.

Bab 6 Integrasi Numerik

279

 y

 y = f ( x)

h

 x0

 x0+h/2

 x1

 x

Gambar 6.6 Kaidah titik tengah

Kaidah titik-tengah gabungan adalah (Gambar 6.7): b

∫  f ( x)dx



a

 x1

 x2

 x n

 x0

 x1

 xn −1

∫  f ( x)dx  + ∫  f ( x)dx  + ... + ∫  f ( x)dx

≈ hf ( x1/2) + hf ( x3/2) + hf ( x5/2) + hf ( x7/2) + ... + hf ( xn-1/2) ≈ h( f 1/2 + f 3/2 +... + f n-1/2) ≈ h

n−1

∑  f 

i+1/2 

P.6.13)

i =0

yang dalam hal ini,  xr +1/2 = a + (r +1/2)h) dan  f r +1/2 = f ( xr +1/2)

280

r = 0,1,2,..,n-1

Metode Numerik

 y

 y = f ( x)

... b

a  x1/2  x3/2  x5/2

...

 x n-3/2  xn -1/2

 x

Gambar 6.7 Kaidah titik-tengah gabungan

Program 6.2  KaidahTitik-tengah

 procedure titik_tengah(a, b : real; n: integer; var I : real); { menghitung integrasi f(x) dalam selang [a, b] dengan jumlah pias sebanyak n. K.Awal : harga a, b, dan n sudah terdefinisi

K.Akhir: I adalah hampiran integrasi yang dihitung dengan kaidah titik-tengah } var

h, x, sigma : real; r : integer;  begin

h:=(b-a)/n; {lebar pias} x:= a+h/2; {titik tengah pertama} sigma:=f(x); for  r:=1 to n-1 do  begin

x:=x+h; sigma:=sigma + f(x) end  ; I:=sigma*h; end  ;

Bab 6 Integrasi Numerik

{ nilai integrasi numerik}

281

6.3.4 Galat Metode Pias Sekarang akan kita hitung berapa besar galat hasil integrasi untuk masing-masing metode. Misalkan  I adalah nilai integrasi sejati dan  I ' adalah integrasi secara numerik maka galat hasil integrasi numerik didefenisikan sebagai  E  = I  – I  '

(P.6.14)

Untuk penurunan galat, kita tinjau galat integrasi di dalam selang [0, h], h

 I  =

∫  f ( x)dx  

(P.6.15)

0

 y

galat  y = f ( x)

h

0

h

 x

Gambar 6.8 Galat kaidah trapesium (bagian yang diarsir)

Untuk setiap kaidah akan kita turunkan galatnya berikut ini.

282

Metode Numerik

6.3.4.1

Galat Kaidah Trapesium

Galat untuk satu buah pias (Gambar 6.8) adalah h

h

∫ 

 E  =

 f ( x) dx  -

2

0

 ( f 0 + f 1)

Uraikan  f ( x) ke dalam deret Taylor di sekitar x0 = 0  f ( x) = f 0 + xf 0' +

1 2

 x2 f 0" +

1 6

 x3 f 0"' + ...

Uraikan  f 1 = f ( x1) = f (h) ke dalam deret Taylor di sekitar x0 = 0  f 1 = f ( x1) = f (h) = f 0 + hf 0' +

1 2

h2 f 0" + ...

Maka, h

 E  =

∫ 

[ f 0 + xf 0' +

0

1 2

 x2 f 0" +

1 6

 x3 f 0"' + ... ]dx -

h 2

 f 0 -

h 2

 [ f 0 + hf 0' + 1 2

h2 f 0" + ...]

h

=  xf 0 + 1/2 x2 f 0' + 1/6 x3 f 0"'+..]  - 1/2 hf 0 - 1/2h f 0 - 1/2 h2 f 0' - 1/4 h3 f 0"' - ... 0

= (hf o + 1/2 h2 f  '0 + 1/6 h3 f  "0 + ...) - (hf 0 + 1/2 h2 f  ' 0 + 1/4 h3 f 0"'+ ...) = -

≈ -

1 12 1 12

h3 f 0 " + ... h3 f  "(t ) ,

0 < t  < h 

(P.6.16)

≈ O(h3) Jadi, h

∫  0

 f ( x) dx ≈

h 2

Bab 6 Integrasi Numerik

( f 0 + f 1) + O(h3)

(P.6.17)

283

Persamaan (P.6.17) ini menyatakan bahwa galat kaidah trapesium sebanding dengan h3. Pernyataan ini dibuat dengan andaian bahwa  f(x) menerus dalam selang [0, h]. Jika tidak, maka galat tidak sebanding dengan h3 [NAK93]. Untuk n buah pias, galat keseluruhan (total) adalah h3

 E tot ≈ -

 ( f 0" + f 1" + f 2" + ... + f  "n-1)

12

yang dapat disederhanakan dengan teorema nilai antara untuk penjumlahan menjadi n −1

h3

 E tot  ≈ -

∑  f  " i

12

i =1

h3

≈ -n

12

 f  "(t )

, a < t  < b 

(P.6.18)

Mengingat h =

b − a n

maka  E tot  ≈ -n

h3 12

≈ - n ≈ -

h3 12

 f  "(t )

b − a h 3  f  "(t ) 12 n (b - a) f  "(t )

(P.6.19)

≈ O(h2) Dengan demikian, b

∫  f ( x)dx a

=

h 2

n−1

 ( f 0 + 2

∑  f  + f  ) + O(h ) i

n

2

(P.6.20)

i =1

Jadi, galat total integrasi dengan kaidah trapesium sebanding dengan kuadrat lebar pias (h). Semakin kecil ukuran ukuran h, semakin kecil pula galatnya, namun semakin banyak jumlah komputasinya. 284

Metode Numerik

Contoh 6.1 3. 4

[GER85] Hitung integral

∫ e dx  dengan kaidah trapesium. Ambil  x

h = 0.2. Perkirakan

1. 8

 juga batas- batas galatnya. Gunakan 5 angka bena. Penyelesaian:

Fungsi integrand -nya adalah  x

 f ( x) = e

Jumlah pias adalah

n = (b-a)/h = (3.4 - 1.8)/0.2 = 8

Tabel data diskritnya adalah sebagai berikut: r

x r 

f(x ) r 

r

x r 

f(x ) r 

0

1.8

6.050

5

2.8

16.445

1

2.0

7.389

6

3.0

20.086

2

2.2

9.025

7

3.2

24.533

3

2.4

11.023

8

3.4

29.964

4

2.6

13.464

 Nilai integrasinya, 3. 4

∫ 

e x dx



h



0.2

2

1. 8

 ( f 0 + 2 f 1 + 2 f 2+ ... + 2 f 6 + 2 f 7 + f 8)

2

[[6.050 + 2(7.389) + 2(9.025) +....+ 2(16.445)

+ 2(20.086) + 2(24.533) + 29.964]

≈  23.994  Nilai integrasi sejatinya adalah 3. 4

∫ 

 x

e x dx  = e

1. 8

= 1.8 = e 3.4 - e1.8 = 29.964 - 6.050 = 23.914  x = 3. 4  x

Galat kaidah trapesium:  E  = -

h2 12

(b - a) f  "(t ) ,

Bab 6 Integrasi Numerik

1.8 < t   < 3.4

285

Karena  f ( x) = e x,  f  '( x) = e x , dan  f  "( x) = e x maka  E  = -

1

 0.22 (3.4 - 1.8) e x , 1.8 < t   < 3.4

12

 x

Karena fungsi  f ( x) = e menaik secara monoton di dalam selang [1.8, 3.4], maka kita dapat menentukan batas-batas galatnya:

e1.8 (min ) = −0.0323(min)  (0.2)  (3.4 - 1.8) ×   E  = 12 e 3.4 (max ) = −0 .1598(max) 1

2

atau -0.0323
View more...

Comments

Copyright © 2017 KUPDF Inc.
SUPPORT KUPDF