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

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

•  TDictionary Custom Sort  5 898

•  Fast Watermark Sources  5 692

•  3D Designer  8 465

•  Sik Screen Capture  6 030

•  Patch Maker  6 475

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

•  ListBox Drag & Drop  5 330

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

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

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

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

•  Canvas Drawing  5 229

•  Рисование Луны  4 950

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

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

•  Paint on Shape  2 427

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

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

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

•  Пазл Numbrix  2 247

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

•  Игра HIP  1 876

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

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

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

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

•  Проверка числового ввода  1 981

•  HEX View  2 289

•  Физический маятник  1 955

 
скрыть

  Форум  

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

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