jackal
|
Carpal Tunnel
|
|
|
|
Рег.: 21.02.2003
|
Сообщений: 3692
|
Из: \\ResidentEvil
|
Рейтинг: 195
|
|
Найти позорный алгоритм для расчета двойного интеграла
02.05.2008 11:25
|
|
|
Найти позорный алгоритм для расчета двойного определенного интеграла от заданной функции, на языке Си.
Область определения пусть квадрат, будем считать что функция определена всюду в нем и нигде не уходит в бесконечность
Погугить чтото не получилось легко, есть ли место где подобные алгоритмы можно скачать?
зы. Я понимаю что написать такое - минут 10 но честно говоря жалко времени (не мне нужно, попросили найти такой), вдруг просто есть сайт какой где лежат решения для стандартных задач.
|
|
|
Tarass
|
парень в форуме
|
|
|
|
Рег.: 26.01.2003
|
Сообщений: 7913
|
|
Рейтинг: 697
|
|
Re: Найти позорный алгоритм для расчета двойного интеграла
[re: jackal]
02.05.2008 12:55
|
|
|
|
halyavin
|
кфмн
|
|
|
|
Рег.: 14.12.2005
|
Сообщений: 916
|
Из: Moscow
|
Рейтинг: 622
|
|
Re: Найти позорный алгоритм для расчета двойного интеграла
[re: jackal]
02.05.2008 13:30
|
|
|
И в чем смысл экономить 10 минут ожидая ответ час?
|
|
KI
|
|
|
|
|
Рег.: 14.11.2003
|
Сообщений: 7908
|
Из: МГУ ГЗ V-806л
|
Рейтинг: 3333
|
|
Re: Найти позорный алгоритм для расчета двойного интеграла
[re: halyavin]
02.05.2008 13:31
|
|
|
|
ksa
|
Умка
|
|
|
|
Рег.: 04.10.2006
|
Сообщений: 14535
|
Из: где-то на белом свете
|
Рейтинг: 7761
|
|
Re: Найти позорный алгоритм для расчета двойного интеграла
[re: jackal]
02.05.2008 15:20
|
|
|
А вы собственно зачем интересуетесь? (c) Почтальон Печкин
Для треугольника есть кубатруная формула 1/3 умножить на площадь треугольника умножить на сумму трех значений функции в серединах сторон.
Обычно этого достаточно с заданной точностью. Если нет, то можно разбить стороны пополам и посчитать по 4 треугольникам с площадью в 4 раза меньше, и т.д. цикл организуется пока разница между двумя соседними разбиениями не станет меньше заданной точности.
Только аккуратно надо с носителем функции (самый первый треугольник должен быть размера порядка его), а то если носитель маленький, то 2 крупных разбиения дадут 0 и ответ выдаст 0. Хотя подобрав функцию типа быстро осцилирующего синуса можно обмануть любую самую извратную программу вычисления интегралов.
Удачи.
|
|