Metode Trapesium pada Metode Numerik

Algoritma Program Pascal untuk Metode Trapesium pada Metode Numerik

program Trapesium;

uses wincrt;

var

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

   jumsig, 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)      E(f,h)');

  repeat

        h:= (b-a)/n;

        jumsig:= 0;

        for i:= 1 to n-1 do

            begin

                 jumsig:= jumsig + f(a+i*h);

            end;

        luas:= (h/2)*(f(a)+f(b)+2*jumsig);

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

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

        n:= 2*n;

        it:= it + 1;

  until it = max;

end.


Preview Program: