METODE SIMPSON KOMPOSIT DALAM BAHASA PASCAL

program simpson_komposit;
uses crt;
var
a,b,h,jtrap,jsimp,xi,isimp : real;
i,n:integer;

function f(x:real):real;
begin
f:=exp((2/3)*ln(x));
end;

begin
clrscr;
writeln('PROGRAM HITUNG INTEGRAL MENGGUNAKAN METODE SIMPSON KOMPOSIT');
write('Masukkan batas bawah    : ');readln(a);
write('Masukkan batas atas     : ');readln(b);
write('Masukkan banyak partisi : ');readln(n);

if n mod 2 = 1 then halt;

h:= (b-a)/(n);

jtrap:=0;
jsimp:=0;

for i:= 1 to n-1 do
 begin
 xi:=a+i*h;
  if i mod 2 = 0 then
   jtrap:= jtrap+f(xi)
 end;
for i:= 1 to n do
 begin
 xi:=a+i*h;
  if i mod 2 = 1 then
   jsimp:= jsimp+f(xi)
 end;

isimp:=(h/3)*(f(a)+f(b)+2*jtrap+4*jsimp);
writeln('Nilai integral adalah : ',isimp:1:5);
readkey;
end.

{coded by : Sylvert Tahalea and Christian Halim}
{dosen pembimbing : Francis Y Rumlawang, S.Si, M.Si}