Навигация
Главная
Поиск
Форум
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
Подключение Mic... 65535
Приложение «Про... 65535
Оператор выбора... 65535
Организация зап... 64646
Создание потоко... 62075
Модуль Forms 61245
Создание отчето... 61201
ТЕХНОЛОГИИ ДОСТ... 57660
Пример работы с... 55172
Имитационное мо... 53052
Реклама
Сейчас на сайте
Гостей: 10
На сайте нет зарегистрированных пользователей

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

Калькулятор на Delphi с переводом в другую систему исчисления + Блок схемы
Моделирование процесса обработки заданий пакетным режимом работы с квант...
Метод половинного деления для нахождения корня уровнения на Turbo Pascal...

Реклама



Подписывайся на 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 Комментариев · 17249 Прочтений · Для печати

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


Комментарии
Александр 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...
Шаблон для новост...

Случайные загрузки
XPButtons
RbControls
Print Grid
Delphi 2005 для W...
Клавиатурный трен...
Введение в станда...
Игра Car [Исходни...
100 компонентов о...
Добавление к ссы...
ATComponents
EMSQuickImport
Задача о 8ми ладьях
Моделирование дви...
Таймер и секундомер
Пример OpenGL гра...
Blobs [Исходник н...
Х. М. Дейтел, П. ...
Язык программиров...
CwstatusBar
Измерение тактово...

Топ загрузок
Приложение Клие... 100384
Delphi 7 Enterp... 83555
Converter AMR<-... 20051
GPSS World Stud... 11313
Borland C++Buil... 11240
Borland Delphi ... 8182
Turbo Pascal fo... 6987
Visual Studio 2... 4970
Калькулятор [Ис... 4418
FreeSMS v1.3.1 3516
Случайные статьи
Монетизация интерн...
Дополнительные нас...
Область Selection
Где найти хорошие ...
Проверка состоятел...
Автоматическое р...
Убрать с экрана пе...
Алгоритм unordered...
Объединение SLAAC ...
Оперативно доступн...
Функция GetPalette...
Файл R.java
Создание общей дир...
Повторяющиеся дейс...
Передача сообщений...
Доставка цветов
Windows Server vNe...
Определение вторич...
но вместо прав адм...
Технологии SQL
В свойствах File S...
DVD-риппинг
Разберемся теперь ...
OpenGL. Шесть куби...
Простой дизайн
Статистика



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


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