ANTICHAT - Показать сообщение отдельно - Студентам с лабораторными сюда

Взятие интеграла методом Симпсона

A, B - интервал интегрирования
N - число точек на интервале
Func - функция, от которой берется интеграл.
Возвращаемое значение - значение интеграла

PS: сделано на Delphi (указывайте среду программирования)

PHP код:


		
			
type

  TFunction = function(X: Extended; Arg: Pointer): Extended;



function Simpson(A, B: Extended; N: Cardinal; Func: TFunction; Arg: Pointer):

  Extended;

var

  h: Extended;

  X: Extended;

  K: Extended;

  I: Integer;

begin

  Assert(N > 0);

  h := 0.5 * (B - A) / N;

  Result := Func(A, Arg);

  X := A + h;

  for I := 1 to 2 * N - 1 do

  begin

    if I mod 2 = 0 then

      K := 2

    else

      K := 4;

    Result := Result + K * Func(X, Arg);

    X := X + h;

  end;

  Result := Result + Func(B, Arg);

  Result := h * Result / 3;

end;