Численные методы решения нелинейных уравнений
Корни нелинейных уравнений
Пусть дано нелинейное уравнение
где — функция, определенная и непрерывная на некотором промежутке. В некоторых случаях на функцию могут быть наложены дополнительные ограничения, например, непрерывность первой и второй производных, что специально оговаривается. Функция может быть задана в виде алгебраического многочлена или трансцендентной функции (тогда ей соответствует алгебраическое или трансцендентное уравнение).
Требуется найти корни нелинейного уравнения (3.1), т.е. числа которые путем подстановки их в (3.1) превращают уравнение в верное числовое равенство. Числа называются также нулями функции .
На практике часто бывает выгодно уравнение (3.1) заменить равносильным ему уравнением (уравнения равносильны, если имеют одинаковые корни):
 (3.2)
где функции — более простые, чем функция . Тогда при задании уравнения в виде (3.1) нулями функции являются точки пересечения с осью (рис. 3.1,д), а при задании в виде (3.2) — абсциссы точек пересечения функций и (рис. 3.1,б).
Замечания
1. Если — алгебраический многочлен, то уравнение (3.1) называется также алгебраическим n-й степени:
 (3.3)
где — действительные числа, коэффициенты уравнения.
 2. На практике встречаются задачи нахождения корней уравнения , левая часть которого задана сеточной функцией (рис. 3.2).
Число есть корень уравнения (3.1) кратности , если при вместе с функцией обращаются в нуль ее производные до (k-1)-го порядка включительно, т.е. , а . Корень кратности к = 1 называется простым. На рис 3.1,с простыми корнями являются , a корни — кратные.
В соответствии с классическим результатом Галуа алгебраическое уравнение (3.1) при не имеет решения в замкнутом (формульном) виде. Сеточные уравнения вообще не имеют формульных решений. Поэтому корни алгебраических , трансцендентных и сеточных уравнений, как правило, определяются приближенно с заданной точностью.
Решение осуществляется в два этапа:
Первый этап. Находятся отрезки , внутри каждого из которых содержится один простой или кратный корень (см. рис. 3.1). Этот этап называется процедурой отделения корней. По сути на нем осуществляется грубое нахождение корней .
Второй этап. Грубое значение каждого корня уточняется до заданной точности одним из численных методов, в которых реализуются последовательные приближения. Порядок (скорость) сходимости метода определяется так же, как в методе простых итераций.
Отделение корней уравнения
Для отделения действительных корней полезно определять заранее число корней, а также верхнюю и нижнюю границы их расположения. Для этого используется ряд теорем.
Теорема 3.1 (о числе корней алгебраического уравнения (3.3)). Алгебраическое уравнение (3.3) n-й степени имеет ровно корней, действительных или комплексных, при условии, что каждый корень считается столько раз, какова его кратность.
Теорема 3.2 (о свойстве парной сопряженности комплексных корней уравнения (3.3)). Если — корень алгебраического уравнения (3.3) кратности , то число также является корнем той же кратности.
Следствие. Алгебраическое уравнение нечетной степени имеет по крайней мере один действительный корень.
Теорема 3.3 (об оценке модулей корней уравнения (3.3)). Пусть
где — коэффициенты уравнения . Тогда модули всех корней уравнения удовлетворяют неравенству
 (3.4)
т.е. корни уравнения расположены в кольце.
Следствие. Числа и являются соответственно нижней и верхней границами положительных корней алгебраического уравнения: . Аналогично числа и служат нижней и верхней границами отрицательных корней уравнения: .
Приведем полезные теоремы, используемые для более точного установления границ действительных корней алгебраических уравнений.
Теорема 3.4 (теорема Лагранжа о верхней границе положительных корней уравнения (3.3)). Пусть и — первый отрицательный коэффициент в последовательности — наибольшая из абсолютных величин отрицательных коэффициентов. Тогда за верхнюю границу положительных корней уравнения (3.3) может быть принято число
![R=1+\sqrt[\LARGE{n-i}]{\frac{C}{a_n}}\,.](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAJAAAAA/BAMAAAD9IlywAAAAG1BMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABp4cHsAAAACXRSTlMA/UWfgslrKBFeBYyVAAAChklEQVRYw92Xz3OaQBTHVxCWo+sP8AhqrUdocOzRtGuaI9jG5mhsND3iJJ32iJM4+bez6y4GGiAddg+ZvnHWGXE+vF/fxwOA/8P2F2mbVwf1N2m7qQ76JCu0iSxQ8NZAylQSSF+ybzgYAVsEpEbMMc+fDYWi1NjX1UdghFNxkGqSwxMK7ZY55JLDEUr2hh4Gisn5Tgh0UEi9Rc9HIdAlPWpNSY1d6wpzIM0ycBo0Uba4QrYUVF+KgOoHNzSLHL9E/AHaQSF6m0T2XQDkxFwh3sRYR+kr1xjjFzmDxakGe+6a1/mZuTJEvUEn+st/r8SnQVHQpEUdN5vOr2Fyf+aruvB5Ou7nTCF5RTB5U6QtzA7D1ThesB+++IWjn9bRc8tB0RxYYEwehss/VuHErgUAoqgcBLsK91kPCkHOvEdmHbHewUZ5ID1QJ6yJtc1vtwB06s/MYx9gfJ4HUqcqjhnoh14kizCin7zQIMafv2FMLsLYGCVlpAqBJ892dLBDnFqWgVLdSm+oxfkg2kah/Upo3Izg+RmS20bIfiXZyX9p6fY5kPe8jZA9KS8/t52VrxDaPjWEGsDD2ZIqCKFpDmhG5/TZS5BxkRBX/7QvGIt25eUIpgcM7LUk7TRGh1dO2BAJ1pUBIhLQlzJARAIqz9p4fC0Aupomja1e6i0B0DY4blm2YgqAdg2+HJH66ZbIrtblo59IXW2IrKFN8CGRes217x7Pqj72TbYcAaUJnPPAX+wrNgOJiDW2YSrh8KbvahVBsAN5Y9/6DydAs3dRVY1kxsE2XitVNfKQHgdrgFcVQV4/yuzJ/aq73+la0ouRM5YE2rYlgWSs6ccdSIoJKTXzFvrm3osL7QlsIGerMQpdQAAAAABJRU5ErkJggg==) (3.5)
Теорема 3.5 (о нижних и верхних границах положительных и отрицательных корней алгебраического уравнения). Пусть — верхняя граница положительных корней уравнения , — верхняя граница положительных корней уравнения , — верхняя граница положительных корней уравнения , — верхняя граница положительных корней уравнения .Тогда положительные корни и отрицательные корни уравнения (3.3) удовлетворяют неравенствам
 (3.6)
Теорема 3.6 (теорема Декарта о количестве действительных корней алгебраических уравнений). Число положительных корней (с учетом их кратностей) алгебраического уравнения равно числу перемен знаков в последовательности коэффициентов (коэффициенты, равные нулю, не учитываются) многочлена или меньше этого числа на четное число. Число отрицательных корней (с учетом их кратностей) алгебраического уравнения равно числу перемен знаков в последовательности многочлена или меньше этого числа на четное число.
Теорема 3.7 (теорема Гюа о необходимом условии действительности всех корней алгебраического уравнения). Если алгебраическое уравнение (3.3) имеет все действительные корни, то квадрат каждого некрайнего коэффициента больше произведения двух его соседних коэффициентов.
Следствие. Если при каком-нибудь выполнено неравенство , то уравнение (3.3) имеет по крайней мере одну пару комплексных корней.
Для отделения корней применяется следующая теорема.
Теорема 3.8. Если функция , определяющая уравнение , на концах отрезка принимает значения разных знаков, т.е. , то на этом отрезке содержится по крайней мере один корень уравнения. Если же непрерывна и дифференцируема и ее первая производная сохраняет знак внутри отрезка , то на находится только один корень уравнения.
Способы отделения корней
В вычислительной практике обычно используются следующие способы отделения корней:
1) средствами машинной графики: функция представляется на дисплее и приближенно определяются отрезки, которым принадлежат точки ; 2) средствами математического анализа с помощью исследования функций и построения графиков (см. рис. 3.1,д); 3) формированием простых функций и таких, что получается равносильное уравнение в виде (3.2), и дальнейшим построением графиков этих функций (см. рис. 3.1,б).
▼ Примеры 3.1-3.3
Пример 3.1. Определить число положительных и отрицательных корней, а также их границы для уравнения 
▼ Решение
В данной задаче . Согласно теореме 3.1 уравнение имеет пять корней. Поскольку , то по следствию из теоремы 3.2 уравнение имеет по крайней мере один действительный корень. Оценим модули корней по теореме 3.3. Так как то или , т.е. все корни лежат внутри данного кольца. По следствию из теоремы 3.3 это означает, что положительные корни удовлетворяют неравенству , а отрицательные — неравенству . Применим теоремы 3.4 и 3.5 для уточнения приведенных результатов. Найдем верхнюю границу положительных корней. Так как — первый отрицательный коэффициент в последовательности , то , а — наибольшая из абсолютных величин отрицательных коэффициентов. Следовательно, . Найдем нижнюю границу положительных корней. Составим уравнение: или (старший коэффициент должен быть положительным). Для этого уравнения , поэтому . Отсюда . Уточним границы отрицательных корней. Составим уравнение:  или  . Для этого уравнения , поэтому . Составим уравнение или . Для этого уравнения , поэтому . Отсюда находим: . Заметим, что данный результат совпадает с полученным ранее. Исследуем структуру корней уравнения. Так как квадрат каждого некрайнего коэффициента больше произведения двух его соседних коэффициентов, то по теореме 3.7 необходимое условие действительности всех корней уравнения выполняется. На основе теоремы 3.6 определим число положительных и отрицательных корней. Выписываем коэффициенты многочлена . Так как число перемен знака , то число положительных корней равно трем или меньше на четное число, т.е. равно 1. Далее выписываем коэффициенты многочлена . Так как число перемен знаков , то число отрицательных корней равно двум или меньше на четное число, т.е. их вообще нет.
Пример 3.2. Отделить корни кубического уравнения . Пример 3.3. Отделить корни уравнения третьего порядка .
 Метод половинного деления
Пусть дано уравнение и отделен простой корень , т.е. найден такой отрезок , что и на концах отрезка функция имеет значения, противоположные по знаку . Отрезок называется начальным интервалом неопределенности, потому что известно, что корень ему принадлежит, но его местоположение с требуемой точностью не определено.
Процедура уточнения положения корня заключается в построении последовательности вложенных друг в друга отрезков, каждый из которых содержит корень уравнения. Для этого находится середина текущего интервала неопределенности , и в качестве следующего интервала неопределенности из двух возможных выбирается тот, на концах которого функция имеет разные знаки (рис. 3.5).
Процесс завершается, когда длина текущего интервала неопределенности становится меньше заданной величины , задающей точность нахождения корня. В качестве приближенного значения корня берется середина последнего интервала неопределенности.
Алгоритм метода половинного деления
1. Найти начальный интервал неопределенности одним из методов отделения корней, задать малое положительное число . Положить .
2. Найти середину текущего интервала неопределенности: .
3. Если , то положить , а если , то принять . В результате находится текущий интервал неопределенности .
4. Если , то процесс завершить: . Приближенное значение корня можно найти по формуле .
Если , положить и перейти к п.2.
Замечания
1. Метод имеет линейную, но безусловную сходимость, и его погрешность за каждую итерацию уменьшается в два раза:
Из последнего соотношения можно оценить число итераций для достижения заданной точности 
Отсюда видно, что, например, для достижения точности при необходимо выполнить примерно десять итераций.
2. К достоинствам метода следует отнести то, что он позволяет найти простой корень уравнения любых непрерывных функций при любых значениях таких, что . Недостатки метода — он не обобщается на системы нелинейных уравнений и не может использоваться для нахождения корней четной кратности.
 Метод хорд
Этот метод при тех же предположениях обеспечивает более быстрое нахождение корня, чем метод половинного деления. Для этого отрезок делится не пополам, а в отношении .
Геометрически метод хорд эквивалентен замене кривой хордой, проходящей через точки и (рис. 3.6).
Уравнение хорды имеет вид . Полагая и , получаем
Предположим, что вторая производная сохраняет постоянный знак, и рассмотрим два случая: (рис. 3.7,д) и (рис. 3.7,б). Случай сводится к рассматриваемому, если уравнение записать в форме: . Первому случаю (см. рис. 3.7,д) соответствует формула (3.7), а второму случаю (3.8) (см. рис. 3.7,б):
 (3.7)
 (3.8)
В первом случае остается неподвижным конец , а во втором случае конец .
Замечание. Для выявления неподвижного конца используется условие , где или . Если неподвижен конец , применяется формула (3.7), а если конец , — формула (3.8).
Пример 3.7. Найти корень уравнения методом хорд с точностью .
▼ Решение
Рассмотрим задачу нахождения корня на отрезке (см. пример 3.2). Так как , a на отрезке , то и, следовательно, имеем второй случай (см. рис. 3.7,б). Положим . Тогда по формуле (3.8) получаем Так как , то положим и продолжим процесс: Так как , то положим и продолжим процесс: Поскольку , положим  Так как , положим  Поскольку , положим  Так как , положим  Поскольку , положим  Так как , то корень уравнения . Из анализа поведения следует, что сходимость метода хорд линейная, однако более быстрая, чем сходимость метода половинного деления.
 Метод простых итераций
Пусть известно, что корень уравнения лежит на отрезке .
Методика решения задачи
1. Уравнение равносильным преобразованием привести к виду . Это преобразование может быть осуществлено различными путями, но для сходимости нужно обеспечить выполнение условия ( — некоторая константа). При этом задача сводится к нахождению абсциссы точки пересечения прямой и кривой (рис. 3.8).
2. Задать начальное приближение и малое положительное число . Положить .
3. Вычислить следующее приближение:
 (3.9)
4. Если , итерации завершаются и . Если , положить и перейти к п.3.
Замечание. В качестве условия завершения итераций при известном значении может быть использовано неравенство
Проблемы сходимости и единственности численного решения, являющиеся главными при использовании этого метода, решаются и исследуются с помощью понятия о сжимающем отображении и теоремы о достаточном условии сходимости метода.
Отображение (функция) называется сжимающим в области с коэффициентом , если для любых двух из выполнено неравенство
 (3.10)
Теорема 3.9 (о сходимости метода простых итераций и единственности получаемого численного решения)
Пусть выполнены условия: 1. Нелинейное уравнение имеет решение . 2. Отображение является сжимающим в области с некоторым коэффициентом . Тогда:
а) решение является единственным решением в области ;
б) последовательность , определяемая по отображению на основе итерационного процесса, сходится к решению со скоростью геометрической прогрессии, т.е. при выборе из условия , где — некоторое малое число, справедливо неравенство
Теорема 3.9 утверждает, что при выполнении условий 1,2 существует окрестность такая, что если взять в этой окрестности и вычислять по формуле (3.9), то в результате с любой наперед заданной точностью можно вычислить , соответствующее искомому (единственному) корню. Но так как эта окрестность неизвестна, то можно взять произвольное . Если при этом вычисляется последовательность , сходящаяся к некоторому значению , то в силу теоремы . Если сходимость отсутствует, то надо взять другое и повторить расчет.
Теорема 3.10 (о достаточном условии сходимости метода простых итераций). Пусть выполнены условия: 1. Функция имеет производные для всех . 2. Существует число , такое, что для всех .
Тогда отображение является сжимающим в с коэффициентом сжатия х и последовательность , определяемая на основе итерационного процесса, сходится к решению , то есть при .
▼ Доказательство
В силу условия 1 справедлива теорема Лагранжа о конечных приращениях:  для  , где — точка, лежащая между и . Но тогда с учетом условия 2 и замены ее мажорантой получим т.е. отображение — сжимающее, и, таким образом, последовательность сходится к некоторому . Но тогда в силу непрерывности в имеет место равенство . Геометрическая интерпретация процесса сходимости и расходимости в зависимости от выполнения или невыполнения достаточного условия сходимости представлена на рис. 3.9. Из рис. 3.9 видно, что при и при (см.рис. 3.9,д,б) итерационные последовательности сходятся к . Причем в первом случае реализуется односторонняя (монотонная) сходимость, а во втором — двусторонняя (немонотонная). При (см.рис. 3.9,в,г) процесс расходится несмотря на то, что точка очень близка к . Процессы на рис. 3.9,а,б соответствуют сжимающему отображению. Преобразование уравнения к равносильному виду может быть выполнено неоднозначно. Рассмотрим универсальные практические приемы равносильного преобразования , дающие сжимающие отображения. 1. Можно заменить уравнение на равносильное , где . Тогда, принимая правую часть этого уравнения за и раскрывая , получаем условие  (3.11) Таким образом, можно найти константу на отрезке так, чтобы удовлетворялись неравенства (3.11). При этом надо стремиться получить такую постоянную , которая бы больше отличалась от нуля, и тогда будет реализовы-ваться более быстрая сходимость. Можно показать, что абсолютная погрешность метода оценивается по формуле (эта оценка выражается через два последовательных приближения)  (3.12) или по формуле (эта оценка выражается через два приближения и )  (3.13) Из оценки следует условие выхода из итерационного процесса. Если значение неизвестно, то выход осуществляется по условию 2. Уравнение заменяется равносильным (где знак в правой части выбирается из условия ):  при  , 3. Можно выразить из уравнения так, чтобы для полученного уравнения выполнялось условие сходимости в окрестности искомого корня. Замечание к доказательству. Выполнение условия не гарантирует близости к точному решению (на рис. 3.10 достигается область, где условие окончания выполняется, но корень расположен достаточно далеко от этой области).
▼ Примеры 3.8-3.10
Пример 3.8. Методом простых итераций с точностью уточнить корень трансцендентного уравнения , причем искомый корень . Пример 3.9. Найти решение уравнения методом простых итераций с точностью и . Пример 3.10. Найти корни уравнения методом простых итераций с точностью .
▼ Решение
1. Преобразуем уравнение к виду . В примере 3.3 была выполнена операция отделения корней и получены отрезки . Можно показать, что на отрезках функция удовлетворяет условию . На отрезке используем другой вид уравнения: . Также легко проверить, что функция удовлетворяет достаточному условию сходимости на отрезке . 2. В качестве начальных приближений выберем: – точку  на отрезке ![[-4;-2]](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAEcAAAAWBAMAAABgTJx+AAAAGFBMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAABoAtTLAAAACHRSTlMAgT/6xBpqnOAbTdYAAADPSURBVCjPndNRDoMgDADQDhF/nSx+k7kDbNQDkMn+WfQAJrvA7v81jDFCMcasHybUF9IWADgUUu/+Zlj6bxVmxljkyuco4ibexfptEiRjNLjikSDWxUiOOTqKeE+QgU5RVGcGTmpdvx0gRcx41Nyj9lpaE3ceDd8Qcb0gdp6ighqyuaYlAfBRCxJzLtfX5ulWVAMITYfJrO3aeOQ3J5LuoDEgXDC3l69zC4Xd9dZiMvEMsewDJBHb9OymMPSu/IvEeAAVsIUuz/3ra8tjz+AH+JkYlh6oY1sAAAAASUVORK5CYII=) ; – точку  на отрезке ![[0,\!5;-2]](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAEgAAAAWBAMAAACRR8fzAAAAG1BMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABp4cHsAAAACXRSTlMAg0P9yhAqZah03qxQAAABIUlEQVQoz43Tz0+DMBQH8O+AMo6dEuPRiXqeFOeOzWCbRxIvHlmWeJ7CEo/EX3+3fdBm0Biyl8AL8Enz+gWAkyq6H35+HqvTFGB3zSXL4VviVR1nLaou13THSb6zvqnS2CBnjjAnJMRPzzg1S680CiR2e7q3vO0vFMSIpEauhFsTktZE3iPcmUYRJ6fQ6lpNf3FE/gHPXTQjlE1ijJLeYmVtI4mv3N90DctyC1GyvGnbCVWhRl/BQn5Bm6GaGkTRtGgn28HdB4NMjWtwjTzKiR1Illz1Y30AZnfjOcpC7SrYI81Hi04ESZaaMBG+xeqdgD2Fa/UCO2EKIUwEcN7VfB51NL1fGjX1ovt2CH3q/juAfK4z5v+hm83w51stT/sN/gCTACkCCfC0MwAAAABJRU5ErkJggg==) ; – точку  на отрезке ![[2;4]](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAACgAAAAWBAMAAABAoMODAAAAGFBMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAABoAtTLAAAACHRSTlMAgkL/vizUYTOBEmkAAACzSURBVBjTbdG5DsIwDAbgP1WOtVAJ1kDFDrhS16CmewdeAMIDVDDw+vRKSKV4yPHJcSwZSEZxXd/Ly7BsAaEDMeAwobDOhHceeaeqxTLyWBq4pYB0Ac8o8hn3jUel0eRi5MwEHII6PnYiESGrodph38RY9lNFYSIUz+Xv99GdPMoum1K5tXTz2IL1Yq4Qmpcv6zSv1/ghIrCxJTRU/fvEnBAOHr8pNAlUOsZdux7H/ZEe2w9WNhaFQe952AAAAABJRU5ErkJggg==) ;. В поставленной задаче  . 3,4. Выполним расчеты по формуле с начальными значениями и и по формуле с начальным значением . Результаты расчетов занесены в табл. 3.6–3.8. В результате получены приближенные значения корней: Обратим внимание на сильное различие в числе итераций, потребовавшихся для нахождения корней (табл. 3.6) и (табл. 3.7), с помощью одной и той же формулы. Заметим, что в окрестности корня значения модуля производной функции равны: С другой стороны, в окрестности корня имеем: Анализ результатов показывает, что чем меньше значения модуля производной , тем быстрее сходимость.
Метод Ньютона для решения нелинейных уравнений
Метод Ньютона (метод касательных, или метод линеаризации) является одним из наиболее популярных численных методов. Он быстро сходится (имеет квадратичную сходимость) и допускает различные модификации, приспособленные для решения векторных задач и сеточных уравнений. Однако этот метод эффективен при весьма жестких ограничениях на характер функции 
 1) существование второй производной функции  на множестве  ; 2) удовлетворение первой производной условию  для всех  ; 3) знакопостоянство  для всех  .
Поэтому его желательно использовать совместно с другими методами, например методом половинного деления, чтобы достигнуть диапазона , где указанные условия начинают выполняться.
Геометрическая интерпретация метода Ньютона состоит в следующем. Задается начальное приближение . Далее проводится касательная к кривой в точке (рис. 3.11), т.е. кривая заменяется прямой линией. В качестве следующего приближения выбирается точка пересечения этой касательной с осью абсцисс. Процесс построения касательных и нахождения точек пересечения с осью абсцисс повторяется до тех пор, пока приращение не станет меньше заданной величины .
Получим расчетную формулу метода Ньютона. Вместо участка кривой (точка соответствует ) возьмем участок — касательную, проведенную в точке ). Для этого отрезка справедливо конечное соотношение:
где — угол наклона касательной в точке к оси абсцисс. Разрешая это соотношение относительно , получаем . Повторяя процесс, находим общую формулу:
 (3.14)
Подчеркнем, что если отбросить итерационный индекс, то (3.14) записывается в виде нелинейного уравнения
 (3.15)
которое, однако, на не равносильно исходному, а является таковым только в одной точке при . Поэтому данный метод не служит разновидностью метода простых итераций.
Применим теперь для вывода формулы (3.14) метод линеаризации. Положим, что итерационный процесс имеет вид
 (3.16)
где — поправка к k-му приближению, которую необходимо найти. Предполагая, что имеет непрерывную вторую производную, разложим по формуле Тейлора относительно точки 
где . Учитывая, что (это соответствует нахождению точки пересечения с осью абсцисс), и оставляя только линейную (относительно ) часть разложения (отсюда и название — метод линеаризации), записываем линейное относительно уравнение
Отсюда выражается поправка . Подставляя в (3.16), получаем (3.14).
Теорема 3.11 (о достаточных условиях сходимости метода Ньютона). Пусть выполняются следующие условия:
1. Функция определена и дважды дифференцируема на . 2. Отрезку принадлежит только один простой корень , так что . 3. Производные на сохраняют знак, и . 4. Начальное приближение удовлетворяет неравенству (знаки функций и в точке совпадают).
Тогда с помощью метода Ньютона (3.14) можно вычислить корень уравнения с любой точностью.
Замечания
1. Метод Ньютона характеризуется вторым порядком сходимости вблизи корня и первым порядком — вдали от него. Данную оценку проверяем для двух случаев, когда находится далеко от корня (это возможно на первых итерациях) и когда располагается вблизи .
Первый случай. Пусть отрезок не мал. Тогда на основе теоремы Лагранжа получим (где )
 (3.17)
Преобразуем с использованием (3.17) итерационную формулу (3.14) (где ):
В силу монотонности последовательности имеем
Таким образом, вдали от корня получаем линейную сходимость
 , где  .
Второй случай. Пусть теперь отрезок мал, т.е. итерации выполняются вблизи корня. Тогда, полагая, что указанные выше предположения 1–4 теоремы 3.10 выполнены, разложим функцию в окрестности корня относительно точки , учитывая члены до второго порядка. Получим
где ( — малая величина). Из данного соотношения выражаем 
Но первые два слагаемых в правой части в соответствии с (3.14) равны . Тогда будем иметь
Из последнего соотношения следует оценка погрешности (k+1)-го приближения через погрешность k-го приближения:
 (3.18)
где принято . Оценка (3.18) свидетельствует о квадратичной сходимости метода касательных вблизи корня. С вычислительной точки зрения это означает, что на каждом приближении количество верных цифр результата удваивается.
2. Для нахождения комплексных корней уравнения можно также использовать (3.14) в форме (где — комплексная переменная)
Метод Ньютона может применяться не только для нахождения простых корней, но для определения кратных корней, т.е. когда на отрезке содержащем корень, не выполняется условие (условие 2 теоремы 3.11). Наряду с теоремой 3.11 удобно также пользоваться следующей теоремой.
Достаточные условия сходимости метода Ньютона
Теорема 3.12 (достаточные условия сходимости метода Ньютона). Пусть: а) дана функция , где — открытый интервал, и ; б) для некоторого выполняется условие при всех из .
Если уравнение имеет решение , то существует некоторое , такое, что если , то последовательность , задаваемая формулой (3.14), существует и сходится к . Более того, справедлива оценка
Здесь обозначение означает, что функция непрерывна по Липшицу с константой на множестве , то есть для любых из .
Замечания
1. Требование теоремы 3.12, чтобы производная  |