|
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.
|