PROGRAM HITUNG LUAS

program hitung_luas;
uses crt;
var
n: byte;
luas,p,l,t : real;
begin
clrscr;
writeln(' ============================');
writeln('|                            |');
writeln('|   PROGRAM MENGHITUNG LUAS  |');
writeln('|    1. PERSEGI PANJANG      |');
writeln('|    2. SEGITIGA             |');
writeln('|    3. LINGKARAN            |');
writeln('|    4. TRAPESIUM            |');
writeln('|    5. SELESAI              |');
writeln('|                            |');
writeln(' ============================');
writeln;
write('Masukkan pilihan anda : ');readln(n);
 case n of
 1: begin
    clrscr;
    writeln('Menghitung Luas Persegi Panjang');
    write('Masukkan panjang : ');readln(p);
    write('Masukkan lebar   : ');readln(l);
    luas := p*l;
    writeln('Luas persegi panjang adalah : ',luas:0:2);
    readkey;
    end;
 2: begin
    clrscr;
    writeln('Menghitung Luas Segitiga');
    write('Masukkkan tinggi : ');readln(t);
    write('Masukkan alas    : ');readln(p);
    luas:= (1/2)*p*t;
    writeln('Luas segitiga adalah : ',luas:0:2);
    readkey;
    end;
 3: begin
    clrscr;
    writeln('Menghitung Luas Lingkaran');
    write('Masukkan jari-jari : ');readln(p);
    luas:=pi*sqr(p);
    writeln('Luas lingkaran adalah : ',luas:0:2);
    readkey;
    end;
 4: begin
    clrscr;
    writeln('Menghitung Luas Trapesium');
    write('Masukkan besar sisi A : ');readln(p);
    write('Masukkan besar sisi B : ');readln(l);
    write('Masukkan tinggi       : ');readln(t);
    luas:= (p+l)*(1/2)*t;
    writeln('Luas trapesium adalah : ',luas:0:2);
    readkey;
    end;
 5: halt;
 else writeln('Pilihan anda tidak ada.');readkey;
 end;
end.

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}