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

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

•  TDictionary Custom Sort  6 019

•  Fast Watermark Sources  5 809

•  3D Designer  8 623

•  Sik Screen Capture  6 144

•  Patch Maker  6 582

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

•  ListBox Drag & Drop  5 432

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

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

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

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

•  Canvas Drawing  5 333

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

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

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

•  Paint on Shape  2 497

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

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

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

•  Пазл Numbrix  2 289

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

•  Игра HIP  1 939

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

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

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

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

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

•  HEX View  2 367

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

 
скрыть

  Форум  

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

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