Нужен способ интегрирования для осцилирующей функциии (на си).

А слабовато было той же Математикой получить аналитическое решение? -Sin[Pi*Cos[x]^2]/(2*Pi), всего то...

Кстати руками не намного сложнее

Конечно же руками данную конкретную функцию не проблемма взять. Просто мне надо расчитать 31*31*6=5766 похожих функций, но значительно более навороченых ( каждая из которых будет минимум в 7 раз длинее вышеуказанной). А конкретную функцию я получил посути посредством отбрасывания всяких частей. Просто если кто скажет мне способ как посчитать численно вышеприведённую функцию, то я смогу посчитать и те что мне надо Если очень надо могу привести точный вид тех функций, но смысла я не вижу. Пример с математикой я привел чтобы показать что mathematica данные функции численно интерирует плохо.

Я никак не могу увидеть в чем сложность численного интегрирования этой функции :0

Да, ваша mathematica меня тоже удивляет. Ответ на отрезке [0,Pi] будет твердый ноль. Так что не должно быть никаких значений там.

Впрочем, может быть имеется в виду несколько другая функция? Например

cos(pi*cos(x^2))*sin(x)*cos(x)?

Тогда все несколько интереснее. Так что?

Имеется именно та функция которую я написал. A в mathematica я интегрировал не через Integrate[...,{x,0,2*Pi}], а через NIntegrate[...,{x,0,2*Pi}], т. е. использовал именно численные алгоритмы. Оно выдаёт ответ порядка 10^-17 ругась при этом. Причём что само интересно это то что и когда я численно интегрирую ( трапециями) ответ получается тоже порядка 10^-17.

Тфу ты. Ещё утром не проснулся. В предыдущем посте пределы интегрирования {x,0,Pi}. А по поводу совпадения порядков - ничего там интересного нету.

Вообще, для численного интегрирования быстро осциллирующих функций есть специальные методы, в хороших книгах по численному анализу они описаны. А готовые программы можно на www.netlib.org поискать.

> 10^-17 a eto i est null (mashinny,blinn)