Cos(x-pi/2), огромное спасибо!
И здорово, что Вам хватило энтузиазма это проверить.
Я тоже Вас поздравляю и желаю приятного пребывания на форуме. Ваши сообщения всегда читаю с огромным интересом.
Расскажу (принося извинения за оффтоп), откуда эта идея.
(Оффтоп)
Когда-то давно надо было вычислить последовательность значений функции Бесселя
, где
фиксировано (и невелико, в пределах сотни), а целый порядок
меняется от
до
. При таких больших порядках
быстро стремится к нулю. И по этой причине рекуррентные формулы мало что давали. Срабатывал эффект малых разностей. При малейшей неточности появлялась «примесь»
, а Нейман с ростом порядка, наоборот, растёт. Очень быстро развивалась неустойчивость, и последовательность «улетала». (Есть асимптотические формулы для больших значений порядка, но по каким-то причинам их решили не применять.)
Тогда возникла мысль пустить рекуррентную формулу в обратном направлении. Беря
с некоторым запасом, в качестве начальных значений задаём
и
почти произвольно. Когда мы доходим до
, то теперь уже всякая примесь Неймана (бывшая в начальных значениях) волшебным образом улетучивается, и при
получается последовательность химически чистых бесселей
с точностью до постоянного множителя. Этот постоянный множитель находится сравнением
, вычисленного таким вот способом, и «библиотечного», и на него потом корректируется вся последовательность. Получилось быстро и эффективно.
Вообще, я никогда не претендовал на знание чего-то выходящего за рамки второго-третьего курса, но довольно часто (как в этом примере) этого хватает.
Эх... вот бы ещё и в реальной жизни научиться так решать проблемы...
Да, да! Задав в будущем предельно счастливую опорную точку и построив от неё прошлое.