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

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

•  TDictionary Custom Sort  5 777

•  Fast Watermark Sources  5 583

•  3D Designer  8 165

•  Sik Screen Capture  5 887

•  Patch Maker  6 364

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

•  ListBox Drag & Drop  5 215

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

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

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

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

•  Canvas Drawing  5 114

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

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

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

•  Paint on Shape  2 335

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

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

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

•  Пазл Numbrix  2 183

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

•  Игра HIP  1 803

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

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

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

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

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

•  HEX View  2 209

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

 
скрыть

  Форум  

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

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



Delphi Sources

Вычисление определенного интеграла методом трапеций с заданной точностью



Автор: Romkin

{ **** UBPFD *********** by delphibase.endimus.com ****
>> Вычисление определенного интеграла методом трапеций с заданной точностью

Просто расчет площади под функцией, параметры: a,b - пределы интегрирования, a<=b
eps - допустимая погрешность, практически гарантируется, что расхождение результата
с истинным значением интеграла не превосходит по модулю указанную величину.
Только не переборщите :-))
intF - подинтегральная функция. Естественно, желательно задавать функции,
интегрируемые в смысле Римана. Объявление смотри в примере.
Примечание: Несобственные интегралы не считаем :-)
Проверок на переполнение нет, да и вообще нет проверок...

Зависимости: Ну какие могут быть зависимости?
Автор:       Romkin, romkin@pochtamt.ru, Москва
Copyright:   Romkin 2002
Дата:        19 ноября 2002 г.
***************************************************** }

unit intfunc;

interface

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

function TrapezeInt(a, b: Double; eps: Double; IntF: TIntFunc): Double;

implementation

function TrapezeInt(a, b: Double; eps: Double; IntF: TIntFunc): Double;
var
  //S - площадь на предыдущей итерации,
  //x - текущее значение аргумента
  //base - высота трапеции
  //n - число трапеций, удваивается на каждой итерации
  S, x, base: Double;
  i, n: Integer;
begin
  //Сначала приближение одной трапецией
  base := b - a;
  Result := (IntF(a) + IntF(b)) / 2 * base;
  eps := eps / 10; //Вообще говоря, величина делителя зависит от функции
  n := 1;
  repeat
    S := Result;
    base := base / 2;
    n := n * 2;
    //Новая площадь вычисляется на основе старой
    Result := Result / 2;
    //Ниже - просто вычисляем площади новых трапеций
    for i := 1 to n div 2 do
    begin
      x := a + base * (i * 2 - 1);
      Result := Result + IntF(x) * base;
    end;
  until abs(S - Result) <= eps;
end;

end.

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

uses intFunc;

function intSin(x: Double): Double;
begin
  Result := sin(x);
end;

procedure TForm1.Button1Click(Sender: TObject);
begin
  Edit1.Text := FloatToStr(TrapezeInt(0, Pi, 0.00001, intSin));
  //результат у меня получился 1.99999990195429 - с запасом
  //Точный ответ - 2.0
end;




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

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

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




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

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