Недавно добавленные исходники

•  DeLiKaTeS Tetris (Тетрис)  4 445

•  TDictionary Custom Sort  6 476

•  Fast Watermark Sources  6 267

•  3D Designer  9 217

•  Sik Screen Capture  6 605

•  Patch Maker  6 991

•  Айболит (remote control)  6 991

•  ListBox Drag & Drop  5 861

•  Доска для игры Реверси  97 024

•  Графические эффекты  7 192

•  Рисование по маске  6 486

•  Перетаскивание изображений  5 357

•  Canvas Drawing  5 736

•  Рисование Луны  5 445

•  Поворот изображения  4 974

•  Рисование стержней  3 533

•  Paint on Shape  2 801

•  Генератор кроссвордов  3 668

•  Головоломка Paletto  2 954

•  Теорема Монжа об окружностях  3 761

•  Пазл Numbrix  2 481

•  Заборы и коммивояжеры  3 166

•  Игра HIP  2 129

•  Игра Go (Го)  2 063

•  Симулятор лифта  2 438

•  Программа укладки плитки  2 111

•  Генератор лабиринта  2 581

•  Проверка числового ввода  2 261

•  HEX View  2 590

•  Физический маятник  2 202

 
скрыть

  Форум  

Delphi FAQ - Часто задаваемые вопросы

| Базы данных | Графика и Игры | Интернет и Сети | Компоненты и Классы | Мультимедиа |
| ОС и Железо | Программа и Интерфейс | Рабочий стол | Синтаксис | Технологии | Файловая система |



Вычисление определённого интеграла методом Симпсона



{ **** UBPFD *********** by delphibase.endimus.com ****
>> Вычисление определённого интеграла методом Симпсона

A, B - границы интегрирования
Eps - заданная относительная точность вычисления
F - подинтегральная функция

Зависимости: нет
Автор:       Dimka Maslov, mainbox@endimus.ru
Copyright:   Dimka Maslov
Дата:        26 ноября 2003 г.
***************************************************** }

type
  TDoubleFunc = function(X: Double): Double;

function Integral(A, B, Eps: Double; F: TDoubleFunc): Double;

  function InternalCalc(A, B: Double; F: TDoubleFunc; N: Integer): Double;
  var
    x, dx: Double;
    i: Integer;
  begin
    dx := (B - A) / N;
    Result := 0;
    x := A;
    for i := 1 to N do
    begin
      Result := Result + dx * (F(x) + 4 * F(x + dx / 2) + F(x + dx)) / 6;
      x := x + dx;
    end;
  end;

var
  N: Integer;
  Prev: Double;
begin
  Result := InternalCalc(A, B, F, 4);
  N := 4;
  repeat
    Prev := Result;
    N := N shl 1;
    Result := InternalCalc(A, B, F, N);
  until (Result = 0) or (Abs((Result - Prev) / Result) < Eps);
end;

Пример использования:

function F(X: Double): Double;
begin
  Result := X * X * X;
end;

procedure TForm1.Button1Click(Sender: TObject);
begin
  Label1.Caption := FloatToStr(Integral(-10, 10, 0.00001, F));
end;




Похожие по теме исходники

Вычисление даты Пасхи

Сортировка методом Хоара

Метод Симпсона