Metode Simpson 1/3 pada Metode Numerik

Algoritma Program Pascal untuk Simpson 1/3 pada Metode Numerik

program Simpson_1_per_3;

uses wincrt;

var

   i, j, n, it, max : integer;

   jumsig1, jumsig2, h, a, b, luas: real;

label 1, 2, 3;

function f(x:real):real;

begin

     f:= 3*x*x + 2*x + 4;

end;

begin

1:

  write('Masukkan banyaknya selang yang anda inginkan: ');readln(n);

  write('Masukkan titik awal                         : ');readln(a);

  write('Masukkan tiitk akhir                        : ');readln(b);

  write('Masukkan banyak trapesium yang dihitung     : ');readln(max);

  writeln;

  writeln('  n        h         T(f,h)');

  repeat

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

        jumsig1:= 0; jumsig2:= 0;

        for i:= 1 to n-1 do

            begin

                 jumsig1:= jumsig1 + f(a + 2*i*h);

            end;

        for i:= 1 to n do

            begin

                 jumsig2:= jumsig2 + f(a + (2*i-1)*h);

            end;

        luas:= (h/3)*(f(a)+f(b)+4*jumsig2+2*jumsig1);

        if n < 10 then writeln('  ',n,'    ',h:0:6,'    ',luas:0:6);

        if n >= 10 then writeln('  ',n,'   ',h:0:6,'    ',luas:0:6);

        n:= 2*n;

        it:= it + 1;

  until it = max;

end.


Preview Program: