Навигация
Главная
Поиск
Форум
FAQ's
Ссылки
Карта сайта
Чат программистов

Статьи
-Delphi
-C/C++
-Turbo Pascal
-Assembler
-Java/JS
-PHP
-Perl
-DHTML
-Prolog
-GPSS
-Сайтостроительство
-CMS: PHP Fusion
-Инвестирование

Файлы
-Для программистов
-Компонеты для Delphi
-Исходники на Delphi
-Исходники на C/C++
-Книги по Delphi
-Книги по С/С++
-Книги по JAVA/JS
-Книги по Basic/VB/.NET
-Книги по PHP/MySQL
-Книги по Assembler
-PHP Fusion MOD'ы
-by Kest
Professional Download System
Реклама
Услуги

Автоматическое добавление статей на сайты на Wordpress, Joomla, DLE
Заказать продвижение сайта
Программа для рисования блок-схем
Инженерный калькулятор онлайн
Таблица сложения онлайн
Популярные статьи
OpenGL и Delphi... 65535
Форум на вашем ... 65535
HACK F.A.Q 65535
Бип из системно... 65535
Гостевая книга ... 65535
Invision Power ... 65535
Содержание сайт... 65535
Организация зап... 65535
Вызов хранимых ... 65535
Программируемая... 65535
Эмулятор микроп... 65535
Подключение Mic... 65535
Создание потоко... 65535
Приложение «Про... 65535
Оператор выбора... 65535
Создание отчето... 65444
Модуль Forms 65205
Пример работы с... 64484
ТЕХНОЛОГИИ ДОСТ... 61861
Имитационное мо... 57771
Реклама
Сейчас на сайте
Гостей: 5
На сайте нет зарегистрированных пользователей

Пользователей: 13,085
новичок: nikita222dolban
Новости
Реклама
Выполняем курсовые и лабораторные по разным языкам программирования
Подробнее - курсовые и лабораторные на заказ
Delphi, Turbo Pascal, Assembler, C, C++, C#, Visual Basic, Java, GPSS, Prolog, 3D MAX, Компас 3D
Заказать программу для Windows Mobile, Symbian

Моделирование круглосуточного интернет кафе на GPSS + Отчет
База данных студентов на Turbo Pascal (Списки) + Пояснительная записка
Медиа плейер на Delphi + Пояснительная записка

Реклама



Подписывайся на YouTube канал о программировании, что бы не пропустить новые видео!

ПОДПИСЫВАЙСЯ на канал о программировании
Вычисление определенного интеграла методом трапеций с заданной точностью
{ **** UBPFD *********** ****
>> Вычисление определенного интеграла методом трапеций с заданной точностью

Просто расчет площади под функцией, параметры: 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;
Опубликовал Kest December 27 2008 23:44:21 · 3 Комментариев · 17977 Прочтений · Для печати

• Не нашли ответ на свой вопрос? Тогда задайте вопрос в комментариях или на форуме! •


Комментарии
Александр June 30 2010 16:38:20
Спасибо. Очень помогло в написании аналогичной программы только на C# smiley
Sem February 09 2011 03:36:02
intSin - выдает ошибку. Почему? Как правильно?
Sem_21@mail.ru
Kest February 09 2011 11:52:09
2Sem

intsin должна быть с аргументом, например intsin(0)
Добавить комментарий
Имя:



smiley smiley smiley smiley smiley smiley smiley smiley smiley
Запретить смайлики в комментариях

Введите проверочный код:* =
Рейтинги
Рейтинг доступен только для пользователей.

Пожалуйста, залогиньтесь или зарегистрируйтесь для голосования.

Нет данных для оценки.
Гость
Имя

Пароль



Вы не зарегистрированны?
Нажмите здесь для регистрации.

Забыли пароль?
Запросите новый здесь.
Поделиться ссылкой
Фолловь меня в Твиттере! • Смотрите канал о путешествияхКак приготовить мидии в тайланде?
Загрузки
Новые загрузки
iChat v.7.0 Final...
iComm v.6.1 - выв...
Visual Studio 200...
CodeGear RAD Stud...
Шаблон для новост...

Случайные загрузки
Иллюстрированный ...
Пользовательская...
Черный круг двига...
База англоязычных...
Midi
C++ для начинающих
Encrypt Decrypt
Шифрование по алг...
Proeffectimage
Измерение тактово...
Анимированное поя...
iChat v.7.0 Final...
Экранная лупа
БД сеть компьютер...
Crypt32
Dreamsoft Progres...
EditButton
Tank [Исходник на...
Abbrevia
PrevInst

Топ загрузок
Приложение Клие... 100470
Delphi 7 Enterp... 87060
Converter AMR<-... 20078
GPSS World Stud... 12822
Borland C++Buil... 11835
Borland Delphi ... 8604
Turbo Pascal fo... 7039
Visual Studio 2... 4999
Калькулятор [Ис... 4800
FreeSMS v1.3.1 3542
Случайные статьи
Генеалогическое де...
Определение размер...
Введение в создани...
Как зарегистрирова...
Работа с нижней по...
Исключения
Добавляем к титлам...
Что собой представ...
Выравнивание текст...
Написанный вручную...
Процедура GetFillP...
Практичность
Разработка общего ...
Анонимность пользо...
В какую объединенн...
Решение проблемы п...
Большой босс
Что означает равен...
Эппель обрабатывал...
метод DelOrders
Загрузить/скачать ...
Маски VLSM, концеп...
Иерархия классов о...
Параллельная магис...
Удаление из файла ...
Статистика



Друзья сайта
Программы, игры


Полезно
В какую объединенную сеть входит классовая сеть? Суммирование маршрутов Занимают ли таблицы память маршрутизатора?