Ошибка "Error in integrate..."

[nickleb]

Здравствуйте. При выполнении скрипта возникает ошибка:

Error in integrate(integrand, 0, 1, subdivisions = subdivisions, rel.tol = rel.tol) :
extremely bad integrand behaviour

Из-за чего? Как её исключить, обработать - "ни в зуб ногой"... Кто-нибудь сталкивался?

[gamm]

Из-за чего? Как её исключить, обработать - "ни в зуб ногой"... Кто-нибудь сталкивался?

можно попробовать просто "продавить" (тыц), но скорее всего функция с особенностями (типа log(0)) и требует учесть асимптотику, а не просто сгущать разбиение. Нужно смотреть на функцию, и читать "Численные методы" Самарского про интегрирование функций с особенностями, например резать диапазон на куски и использовать какие-то квадратуры типа Гауссовой, либо аналитические/асимптотические аппроксимации около особенностей ...

[nickleb]

Из-за чего? Как её исключить, обработать - "ни в зуб ногой"... Кто-нибудь сталкивался?

можно попробовать просто "продавить" (тыц), но скорее всего функция с особенностями (типа log(0)) и требует учесть асимптотику, а не просто сгущать разбиение. Нужно смотреть на функцию, и читать "Численные методы" Самарского про интегрирование функций с особенностями, например резать диапазон на куски и использовать какие-то квадратуры типа Гауссовой, либо аналитические/асимптотические аппроксимации около особенностей ...

gamm, спасибо! разрядил прОфили и "прокатило"... да у меня уже по-своему было просчитано... решил с R сравнить, т.к. там дальше в R шаги - так я с ними уже работал... Самарский, Гулин - у них крайнее издание в 89 ? или есть и попозже издания?

[gamm]

"Попозже" не надо, зватит изданий 70х годов, по которым меня учили. И в R численные методы сделаны абы как, "чтобы было", лучше Матлаб или Питон.

[ Сообщение с мобильного устройства ]

[nickleb]

... в R численные методы ...

[ Сообщение с мобильного устройства ]

https://cran.r-project.org/web/views/Nu ... atics.html