Nxx писал(а):
Ну посколько интегралы не берутся, введите две функции:
И ищите
стандартно методом Ньютона как предел последовательности
Функция выпукла, так что последовательность мгновенно сойдётся (за где-нибудь пять-шесть итераций), если удачно выбрать начальное приближение. Оно должно быть достаточно близким к единице, но всё же не слишком близким, иначе объём вычислений увеличится раза в два или (жутко даже сказать) в три. Скорее всего,
вполне подойдёт.
Другое дело, что пакет может обидеться на попытку заставить его считать эти интегралы, поскольку логарифм в нуле всё-таки не определён. Тогда самый рабоче-крестьянский выход -- это попросту отойти маленько от нуля и интегрировать, скажем, начиная с точки
; соответственно, и
получится с точностью до четвёртого знака. Если захочется всё же большей точности -- ну тут могут быть разные варианты развития. Во всяком случае, уточнять имеет смысл только саму
, а вот уточнение
практического влияния на метод Ньютона не окажет.