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

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

Моделирование работы участка термической обработки шестерен на GPSS + По...
Меры близости на векторах в 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 Комментариев · 18163 Прочтений · Для печати

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


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

Случайные загрузки
Matrix2D
WordReport
PrevInst
ADVstatusbar
Архив значков
Ранги для форума
Размещение элемен...
Основы программир...
IconCut [Исходник...
Пример работы с ф...
ЯЗЫК ПРОГРАММИРОВ...
PDJXPPack
Карта сайта
RAS
База игр
Фильтры изображений
Стелтинг Стивен, ...
начисление процен...
MPTools
Удаление своего EXE

Топ загрузок
Приложение Клие... 100490
Delphi 7 Enterp... 88213
Converter AMR<-... 20084
GPSS World Stud... 13851
Borland C++Buil... 12154
Borland Delphi ... 8712
Turbo Pascal fo... 7058
Visual Studio 2... 5007
Калькулятор [Ис... 4936
FreeSMS v1.3.1 3547
Случайные статьи
Сборка
Работа в режиме HT...
Вывод текста в гра...
String length mism...
Декларативная и пр...
Дизайнерские ресур...
Реализация перемещ...
Что нужно для двои...
Об этой книге XXXII
Структура распрост...
write(X)
HTPC
Разборка приставки...
Общие подходы к пр...
Содержание
Казино Вулкан Слотс
Классы-оболочки вс...
Программа вычитает...
TEST (ПРОВЕРИТЬ)
Установка дистанци...
Описание GPSS-прог...
Функции для выполн...
Поиск максимальног...
Объектные модели M...
Глава 10
Статистика



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


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