Приятно вести обсуждение с солидной персоной. Во-первых, Математика 11.2 находит интеграл в замкнутом виде:
Вышеуказанная формула не эффективна, ибо интеграл выражается через гипергеометрические функции комплексного аргумента, т.е. суммы степенных рядов (а то и через аналитические продолжения этих сумм).
Ваш подход, который состоит в нахождении асимптотических выражений, приводит к несравнимо более простым формулам. Однако ваши выкладки надо поддержать оценками, ибо
Код:
J1 := int((1-(1/2)*t^2+(1/24)*t^4)*ln(t)/t, t = B .. 3);
J1 := -(1/96)*ln(B)*B^4+(1/384)*B^4+(1/4)*ln(B)*B^2-(1/8)*B^2-(1/2)*ln(B)^2-(45/32)*ln(3)+117/128+(1/2)*ln(3)^2
evalf(eval(J1, B = 0.1e-1));
-10.63131040
J2 := -sin(B)*ln(B)/B+cos(B)*(ln(B)-1)/B^2;
J2 := -sin(B)*ln(B)/B+cos(B)*(ln(B)-1)/B^2
evalf(eval(J2, B = 3));
-0.06252599460
а Математика производит
Код:
NIntegrate[Cos[t]*Log[t]/t, {t, 1/100, Infinity}]
-10.8486
Это же числовой результат получается и аналитически. Исполненные коды могу представить через Dropbox.