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

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

•  TDictionary Custom Sort  5 893

•  Fast Watermark Sources  5 688

•  3D Designer  8 458

•  Sik Screen Capture  6 025

•  Patch Maker  6 471

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

•  ListBox Drag & Drop  5 324

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

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

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

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

•  Canvas Drawing  5 225

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

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

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

•  Paint on Shape  2 426

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

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

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

•  Пазл Numbrix  2 246

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

•  Игра HIP  1 874

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

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

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

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

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

•  HEX View  2 289

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

 
скрыть

  Форум  

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

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