algorytm.org

Implementacja w Delphi/Pascal



Baza Wiedzy
wersja offline serwisu przeznaczona na urządzenia z systemem Android
Darowizny
darowiznaWspomóż rozwój serwisu
Nagłówki RSS
Artykuły
Implementacje
Komentarze
Forum
Bookmarki






Sonda
Implementacji w jakim języku programowania poszukujesz?

Całkowanie numeryczne - metoda Simpsona - Implementacja w Delphi/Pascal
Ocena użytkownikóww: *****  / 3
SłabyŚwietny
Nadesłany przez Tomasz Lubiński, 18 października 2009 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 Simpsona - Delphi/metoda_simpsona.dpr:
// Calkowanie numeryczne - metoda Simpsona
// www.algorytm.org
// (c)2007 by Tomasz Lubinski

program metoda_simpsona;
{$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, s, x: 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;
        s := 0;
        for i:=1 to n-1 do
        begin
                x := xp + i*dx;
                s := s + func(x - dx / 2);
                calka := calka + func(x);
        end;
        s := s + func(xk - dx / 2);
        calka := (dx/6) * (func(xp) + func(xk) + 2*calka + 4*s);

	writeln('Wartosc calki wynosi w przyblizeniu ' + FloatToStr(calka));
        readln;
end.
Dodaj komentarz