Contoh Program Pascal untuk Metode Simpson 3/8

Metode Integrasi Simpson 3/8

Metode integrasi simpson 3/8 adalah salah satu metode integrasi numerik guna menghitung suatu luasan kurva yang menggunakan fungsi polinomial dengan orde tiga yang melalui empat titik. Metode Simpson 1/3 biasanya lebih disukai karena mencapai ketelitian orde 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. Namun, metode ini tidak begitu baik karena adanya kesalahan yang cukup besar. Untuk itu, kedua metode dapat digabung, yaitu sejumlah genap pias menggunakan metode Simpson 1/3 dan kelipatan 3 pias sisanya menggunakan metode Simpson 3/8.

Berikut ini contoh program metode integrasi simpson 3/8 dalam bahasa pemrograman Pascal

Baca juga: Program Pascal untuk Metode Trapesium

program coba;

uses crt,sysutils;

var bb,ba,h,int,k:real;

i,n:integer;

t1,t2:tdatetime;

hr,min,sec2,ms2:word;

function f(x:real):real;

begin

f:=x*x;

end;

procedure garis;

begin

writeln('===========================================');

end;

begin

clrscr;

write('Masukkan batas bawah: '); readln(bb);

write('Masukkan batas atas: '); readln(ba);

write('Masukkan banyak segmen (kelipatan 3): '); readln(n);

garis;

t1:=now;

writeln(timetostr(t1)); //waktu awal

garis;

h:=(ba-bb)/n;

int:=f(bb)+f(ba);

for i:=1 to (n-1) do

begin

    k:=bb+i*h;

    if i mod 3 = 0 then

    int:=int+2*f(k)

    else

    int:=int+3*f(k);

end;

int:=int*h*0.375;

writeln('Hasil = ',int:5:2);

garis;

t2:=now;

writeln(timetostr(t2)); //waktu akhir

DecodeTime(t2-t1,hr,min,sec2,ms2);

writeln('Processing time = ',sec2*1000+ms2,' milisekon');

readkey;

end.