informatyki i systemów mikroprocesorowych

program metoda_trapezow;
uses crt;
var
 a,b:real;
 n:integer;
function funkcja(x:real):real;
begin
funkcja:=sin(x);
end;
function poletrapez(a,b:real;n:integer):real;
var
 p:real;
 k,xi:real;
 i:integer;
begin
 p:=0;
 k:=(b-a)/n;
 xi:=a+k;
 while xi<b do
 begin
  p:=p+funkcja(xi);
  xi:=xi+k;
 end;
 poletrapez:=(Fun(a)+Fun(b)+2*p)*k/2;
end;
begin
 clrscr;
 write('Podaj dolna granice przedzialu: ');
 readln(a);
 write('Podaj gorna granice przedzialu: ');
 readln(b);
 write('Podaj liczbe iteracji: ');
 readln(n);
 writeln('Pole=',poletrapez(a,b,n));
 repeat until keypressed;
end.