Ocena użytkownikóww: ***** / 3
Nadesłany przez Tomasz Lubiński, 27 maja 2007 01:00
Kod przedstawiony poniżej przedstawia główną część rozwiązania problemu.
Pobierz pełne rozwiązanie.Jeżeli nie odpowiada Ci sposób formatowania kodu przez autora skorzystaj z pretty printer'a i dostosuj go automatycznie do siebie.
Metoda trapezow - Delphi/metoda_trapezow.dpr:
// Całkowanie numeryczne - metoda trapezow
// www.algorytm.org
// (c)2007 by Tomasz Lubinski
program metoda_trapezow;
{$APPTYPE CONSOLE}
uses
SysUtils;
//funkcja dla ktorej obliczamy calke
function func(x: Real): Real;
begin
result := x*x + 3;
end;
var
xp, xk, dx, calka: Real;
n, i: Integer;
begin
writeln('Podaj poczatek przedzialu calkowania');
readln(xp);
writeln('Podaj koniec przedzialu calkowania');
readln(xk);
writeln('Podaj dokladnosc calkowania');
readln(n);
dx := (xk - xp) / (n * 1.0);
calka := 0;
for i:=1 to n-1 do
calka := calka + func(xp + i * dx);
calka := calka + ((func(xp) + func(xk)) / 2);
calka := calka *dx;
writeln('Wartosc calki wynosi w przyblizeniu ' + FloatToStr(calka));
readln;
end.
