integral metode simpson dengan program pascal

Hai rek, hehe,kemaren ada tugas kampus matkul program komputer, jadi tak share aja deh cara buat nyari integral dengan metode simpson, gimana yah caranya?check this out :)

uses crt;

var

a,b,h,t,s,xi,p : real;

i,n:integer;

function f(x:real):real;

begin

f:=1/sqrt(2*3.14)*exp(-0.5*x*x);

end;

begin

write('batas atas: ');

readln(a);

write('batas bawah: ');

readln(b);

write('jumlah iterasi: ');

readln(n);

if n mod 2 = 1 then halt;

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

t:=0;

s:=0;

for i:= 1 to n-1 do

 begin

 xi:=a+i*h;

  if i mod 2 = 0 then

   t:= t+f(xi)

 end;

for i:= 1 to n do

 begin

 xi:=a+i*h;

  if i mod 2 = 1 then

   s:= s+f(xi)

 end;

p:=(h/3)*(f(a)+f(b)+2*t+4*s);

writeln('Nilai integral adalah : ',p:1:3);

readln;

end.

Gitu sintagnya. Sintag diatas digunakan untuk menyelesaikan soal dalam mencari integral dari

Hasilnya adalah 0,475 kan?

Untuk soal lainnya tinggal mengganti fungsi pada sintag diatas.

Selamat mencoba, semoga sukses :)