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

•  DeLiKaTeS Tetris (Тетрис)  3 927

•  TDictionary Custom Sort  6 030

•  Fast Watermark Sources  5 820

•  3D Designer  8 640

•  Sik Screen Capture  6 152

•  Patch Maker  6 591

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

•  ListBox Drag & Drop  5 442

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

•  Графические эффекты  6 796

•  Рисование по маске  5 966

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

•  Canvas Drawing  5 342

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

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

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

•  Paint on Shape  2 507

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

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

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

•  Пазл Numbrix  2 295

•  Заборы и коммивояжеры  2 985

•  Игра HIP  1 945

•  Игра Go (Го)  1 858

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

•  Программа укладки плитки  1 931

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

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

•  HEX View  2 376

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

 
скрыть

  Форум  

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

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



Delphi Sources

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



{ **** 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;




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

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

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

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




Copyright © 2004-2025 "Delphi Sources" by BrokenByte Software. Delphi World FAQ

Группа ВКонтакте