в примерах на EXCEL 1 страница

под редакцией профессора В.Ф.Корнюшко

Москва

1999 г.

ББК 32.97

УДК 681.3

БУРЛЯЕВ В.В.

Численные методы в примерах на EXCEL.

.

Под редакцией проф. Корнюшко В.Ф.

Рецензент - д.т.н., профессор Бахвалов Л.А.

Пособие предназначено для самостоятельного изучения дисциплины “Численные методы расчетов” при подготовке к выполнению лабораторных работ. Оно должно дать студенту основные понятия о численных методах вычислительной математики с использованием современных компьютеров и доступных программных средств.

Основное внимание уделено тщательно подобранным примерам, позволяющим наиболее ярко проиллюстрировать те или иные особенности каждого метода. Все примеры выполнены на одном из самых мощных современных программных средств - табличном процессоре EXCEL, входящим в состав широко распространенного пакета MICROSOFT OFFICE.

Пособие охватывает все темы раздела учебной программы указанной дисциплины. Кроме методов, входящих в учебную программу, в пособии описаны алгоритмы и вычислительные процедуры встроенных в EXCEL специальных подпрограмм и функций, позволяющих реализовать те или иные численные методы, например, матричные вычисления, линейный регрессионный анализ, метод сопряженных градиентов, линейное программирование и т.п.

СОДЕРЖАНИЕ

1. РЕШЕНИЕ НЕЛИНЕЙНОГО УРАВНЕНИЯ с одним неизвестным. 5

1.1 Отделение корней.................................................................................................. 5

Пример 1.1.................................................................................................................... 5

1.2 Уточнение корней: метод итераций.................................................................... 6

Пример 1.2.................................................................................................................... 7

1.3 Уточнение корней: метод Ньютона..................................................................... 8

Пример 1.3.................................................................................................................... 9

1.4. Уточнение корней: метод бисекции (деления отрезка пополам)................ 10

Пример 1.4.................................................................................................................. 10

1.5 Уточнение коней: подпрограмма EXCEL “ Подбор параметра ”................... 12

2. СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ................. 13

2.1. Матричный метод............................................................................................... 13

Пример 2.1.................................................................................................................. 14

2.2. Метод приближенных вычислений.................................................................. 15

Пример 2.2.................................................................................................................. 16

2.3. Метод Гаусса – Зайделя...................................................................................... 18

Пример 2.3.................................................................................................................. 18

3. РЕШЕНИЕ СИСТЕМ НЕЛИНЕЙНЫХ УРАВНЕНИЙ.................................. 20

3.1. Выбор начальных приближений....................................................................... 20

Пример 3.1.................................................................................................................. 20

3.2 Метод Ньютона.................................................................................................... 21

Пример 3.2.................................................................................................................. 22

3.3. Метод итераций.................................................................................................. 23

Пример 3.3.................................................................................................................. 24

4. ЧИСЛЕННЫЕ МЕТОДЫ ОДНОМЕРНОЙ ОПТИМИЗАЦИИ..................... 25

4.1. Метод дихотомии................................................................................................ 25

Пример 4.1.................................................................................................................. 26

4.2. Метод золотого сечения..................................................................................... 27

Пример 4.2.................................................................................................................. 28

4.3. Встроенная подпрограмма EXCEL “ Поиск решения ”.................................... 29

5. МНОГОМЕРНЫЕ ЗАДАЧИ ОПТИМИЗАЦИИ................................................ 30

5.1. Безусловная оптимизация: метод покоординатного спуска.......................... 30

Пример 5.1.................................................................................................................. 31

5.2. Безусловная оптимизация: метод наискорейшего спуска.............................. 32

Пример 5.2.................................................................................................................. 33

5.3. Безусловная оптимизация: подпрограмма EXCEL “Поиск решения”.......... 35

5.4. Условная оптимизация: метод штрафных функций....................................... 35

Пример 5.3.................................................................................................................. 37

5.5. Условная оптимизация: подпрограмма EXCEL “Поиск решения”............... 38

5.6. Условная оптимизация: линейное программирование.................................. 39

Пример 5.4.................................................................................................................. 39

6. МЕТОД НАИМЕНЬШИХ КВАДРАТОВ........................................................... 43

Пример 6.1.................................................................................................................. 44

Пример 6.2.................................................................................................................. 45

7. ВЫЧИСЛЕНИЕ ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ..................................... 48

Пример 7.1.................................................................................................................. 49

8. ЧИСЛЕННОЕ РЕШЕНИЕ ОБЫКНОВЕННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ......................................................................................................................................... 51

8.1. Метод Эйлера...................................................................................................... 51

Пример 8.1.................................................................................................................. 51

8.2. Метод Рунге-Кутта четвертого порядка........................................................... 52

Пример 8.2.................................................................................................................. 53

8.3. Метод прогноза и коррекции: метод Адамса................................................... 53

Пример 8.3.................................................................................................................. 54

9. РЕШЕНИЕ СИСТЕМ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ. 55

9.1. Задача Коши........................................................................................................ 55

Пример 9.1.................................................................................................................. 56

9.2. Краевая задача: метод стрельбы........................................................................ 57

Пример 9.2.................................................................................................................. 57

9.3. Краевая задача: метод прогонки........................................................................ 57

Пример 9.3.................................................................................................................. 58

10. ЧИСЛЕННОЕ РЕШЕНИЕ УРАВНЕНИЙ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ 60

Пример 10.1................................................................................................................ 61

1. РЕШЕНИЕ НЕЛИНЕЙНОГО УРАВНЕНИЯ с одним неизвестным.

Уравнение с одним неизвестным можно записать в каноническом виде

f(x) = 0

Решение уравнения заключается в нахождении корней, т.е. таких значений х, которые обращают уравнение в тождество. В зависимости от того, какие функции входят в уравнение, разделяют два больших класса уравнений - алгебраические и трансцендентные. Функция называется алгебраической, если для получения значения функции по данному значению х нужно выполнить арифметические операции и возведение в степень. К трансцендентным функциям относятся показательная, логарифмическая, тригонометрические прямые и обратные и т.п.

Найти точные значения корней можно лишь в исключительных случаях. Как правило, используются методы приближенного вычисления корней с заданной степенью точности Е. Это означает, что если установлено, что искомый корень лежит внутри интервала [a,b], где a - левая граница, а b - правая граница интервала, и длина интервала (b-a) <= E, то за приближенное значение корня можно принять любое число, находящееся внутри этого интервала.

Процесс нахождения приближенных значений корней разбивается на два этапа: 1) отделение корней и 2) уточнение корней до заданной степени точности. Рассмотрим эти этапы подробнее.

1.1 Отделение корней.

Любой корень уравнения считается отделенным на отрезке [a,b], если на этом отрезке исследуемое уравнение не имеет других корней.

Отделить корни - это значит разбить всю область допустимых значений х на отрезки, в каждом из которых содержится только один корень. Эту операцию можно провести двумя способами - графическим и табличным. Если функция f(x) такова, что можно легко построить качественный график ее изменения, то по этому графику достаточно грубо находятся два числа, между которыми лежит одна точка пересечения функции с осью абсцисс. Иногда с целью облегчения построения, целесообразно представить исходное каноническое уравнение в виде f₁(x) = f₂(x), затем построить графики этих функций, причем абсциссы пересечения графиков и служат корнями данного уравнения.

При наличии компьютера наиболее распространен табличный способ отделения корней. Он заключается в табулировании функции f(x) при изменении х от некоторого значения х_нач до значения х_кон с шагом dx. Задача заключается в том, чтобы найти в этой таблице такие два смежных значения х, для которых функция имеет разные знаки. Предположим, что такие два значения a и b=a+dx найдены, т.е. f(a)*f(b)<0. Тогда согласно теореме Больцано-Коши внутри отрезка [a,b], если функция f(x) непрерывна, существует точка с, в которой f(c)=0. EXCEL позволяет легко реализовать оба способа отделения корней. Рассмотрим их на примере.

Пример 1.1.

Требуется отделить корни уравнения

e^x-10x = 0

Для этого надо протабулировать функцию f(Х) = exp(Х) - 10*Х, записанную по правилам EXCEL, и построить ее график при изменении Х от какого-то Х_нач до Х_кон с шагом dХ. Пусть эти значения сначала будут таковы: Х_нач = 0, Х_кон = 5, dХ = 0,5. Если в этих пределах изменения Х нам не удастся отделить ни одного корня, тогда надо будет задать новые начальное и конечное значения х и, может быть, изменить шаг.

Для построения таблицы целесообразно воспользоваться специальной подпрограммой ТАБЛИЦА. Для этого на новом рабочем листе в ячейке B1 введем текст: ОТДЕЛЕНИЕ КОРНЕЙ. Затем в ячейку А2 введем текст: x, а в смежную ей ячейку В2 - текст: f(x). Далее оставим ячейку А3 пустой, но в ячейку В3 введем формулу исследуемой функции по правилам EXCEL, а именно

=EXP(A3)-10*A3

Затем заполним числовой ряд изменений X в строках А4:A14 от 0 до 5 с шагом 0,5.

Выделим блок ячеек А3:B14. Теперь дадим команду меню Данные- Таблица. Результаты табулирования будут помещены в блок ячеек В4:В14. Для того чтобы сделать их более наглядными, нужно отформатировать блок В4:B14 так, чтобы отрицательные числа окрашивались в красный цвет. В этом случае легко найти два смежных значения X, для которых значения функции имеют разные знаки. Их и надо принять за концы интервала отделения корней. В нашем случае таких интервалов, как видно из таблицы два - [0;0,5] и [ 3,5;4].

Далее следует построить график нашей функции, выделив блок А4:B14 и вызвав Мастер Диаграмм. В результате получим на экране диаграмму изменения f(X), из которой видны следующие интервалы отделения корней [0;1] и [3;4].

Если изменять теперь числовые значения х в блоке А4:A14 то значения функции в ячейках B4:B14и график будут изменяться автоматически.

1.2 Уточнение корней: метод итераций.

Для уточнения корня методом итераций должно быть задано:

1) уравнение f(X) = 0, причем f(X) должно быть задано в виде формулы,

2) числа a - левая граница и b - правая граница интервала, внутри которого лежит один корень,

3) число Е - заданная точность получения корня.

Сам метод можно разбить на два этапа:
а) переход от канонического вида записи уравнения f(X)=0 к итерирующему виду X = g(X),
б) вычислительная итерирующая процедура уточнения корня.

Перейти от канонического вида уравнения к итерирующему можно различными способами, важно лишь чтобы при этом выполнялось достаточное условие сходимости метода: çg’(X)ç<1 на [a,b], т.е. модуль первой производной итерирующей функции должен быть меньше 1 на интервале [a,b]. Причем чем меньше этот модуль, тем больше скорость сходимости.

Вычислительная процедура метода состоит в следующем. Выбираем начальное приближение, обычно равное Х₀ = (a+b)/2. Затем вычислим X₁=g(X₀) и D= X₁- X₀. Если модуль D <= E, то X₁ является корнем уравнения. В противном случае переходим ко второй итерации: вычисляем Х₂=g(X₁) и новое значение D=X₂- X₁. Опять проводим проверку на точность и при необходимости продолжаем итерации. Если g(X) выбрано правильно и удовлетворяет достаточному условию сходимости, то эта итерирующая процедура сойдется к корню. Следует отметить, что от знака g’(X) зависит характер сходимости: при g’(X)>0 сходимость будет монотонной, т.е. с увеличением итераций D будет приближаться к Е монотонно (не меняя знака), в то время как при g’(X)<0 сходимость будет колебательной, т.е. D будет приближаться к Е по модулю, меняя знак на каждой итерации.

Рассмотрим реализацию метода итераций на EXCEL на примере.

Пример 1.2

Уточним методом итераций значение корней, отделенных в примере 2.1. Итак пусть f(X)= exp(X) - 10*X, для первого корня a=0 и b=0,5. Пусть Е=0,00001. Как выбрать итерирующую функцию? Например, так g(X)=0,1*exp(X). На интервале [a,b] çg’(X)ç<1 и достаточное условие сходимости выполняется. Кроме того, эта производная >1 на интервале [a,b] и характер сходимости будет монотонный.

Запрограммируем метод итераций для этого примера на том же рабочем листе, где мы проводили отделение корней. В ячейку А22 внесем число, равное 0. В ячейку В22 запишем формулу =0,1*EXP(A22), а в ячейку С22 формулу =А22- В22. Таким образом 22 строка содержит данные по первой итерации. Чтобы получить в строке 23 данные по второй итерации, скопируем содержимое ячейки В22 в ячейку А23, записав в А23 формулу =В22. Далее надо скопировать формулы ячеек В22 и С22 в ячейки В23 и С23. Для получения данных всех остальных итераций надо выделить ячейки А23,В23,С23 и скопировать их содержимое в блок А24:C32. После этого следует проанализировать изменение D = Х - g(X) в столбце С, найти D<0,00001 по модулю и выбрать соответствующее ему значение Х из столбца А. Это и есть приближенное значение корня.

Для большей наглядности можно построить диаграмму для метода итераций. Выделяя блок А22:С32 и используя Мастер диаграмм, получим три графика изменения Х,g(X) и D в зависимости от номера итераций, для чего на шаге 3 из 5 выберем формат 2, а на шаге 4 из 5 построения диаграммы нужно отвести ноль столбцов для меток оси Х. Теперь хорошо виден монотонный характер сходимости D.

Для уточнения второго корня этого уравнения на интервале [3,5;4], нужно выбрать другую итерирующую функцию, такую чтобы ее первая производная была по модулю меньше единицы. Выберем g(X)= LN(X)+LN(10). В ячейку А22 внесем новое Х0=3,75, а в ячейку В22 - новую формулу =LN(A22)+LN(10). Скопируем формулу из В22 в блок В23:В32 и сразу получим новые данные и перестроенную диаграмму. Определим приближенное значение второго корня.

1.3 Уточнение корней: метод Ньютона.

Для уточнения корня методом Ньютона должно быть дано:

1) уравнение f(X) = 0, причем f(X) должно быть задано в виде формулы,

2) числа a - левая граница и b - правая граница интервала, внутри которого лежит один корень,

3) число Е - заданная точность получения корня,

4) функция f(X) должна быть дважды дифференцируемой, причем формулы f’(X) и f”(X) должны быть известны.

Метод состоит в итерационных вычислениях последовательности

X_i+1 = X_i - f(X_i)/f’(X_i), где i=0,1,2,...,

исходя из начального приближения Х₀, принадлежащего интервалу [a,b] и удовлетворяющего условию f(X₀)*f”(X₀)>0. Достаточные условия сходимости метода заключаются в том, что первая и вторая производные исследуемой функции должны сохранять знак на интервале [a,b]. В качестве начального приближения выбирают обычно или a, или b, в зависимости от того, кто из них соответствует формуле выбора Х₀.

Метод Ньютона допускает простую геометрическую интерпретацию. Если через точку с координатами (X_i;f(X_i)) провести касательную к кривой f(X), то абсцисса точки пересечения этой касательной с осью 0Х и есть очередное приближение корня Х_i+1.

Метод Ньютона можно рассматривать как некоторую модификацию метода итераций, дающую наилучшую итерирующую функцию g(X) на каждом шаге итерации. Проведем следующие преобразования с исходным каноническим уравнением f(X)=0. Умножим левую и правую его части на некоторое число l, отличное от нуля. Затем прибавим слева и справа по Х. Тогда будем иметь

Х = g(X) = Х +l*f(X).

Дифференцируя g(X), получим g’(X) = 1 + l*f’(X). Из достаточного условия сходимости метода итераций çg’(X)ç<1. Потребуем, чтобы на i-том шаге итерации сходимость была самой быстрой, т.е. çg’(X_i)ç =0. Тогда l=-1/ f’(X_i) и мы пришли к методу Ньютона.

Вычислительная процедура метода состоит в следующем. Выбираем начальное приближение X₀, обычно равное a или b. Затем вычислим X₁= X₀ - f(X₀)/f’(X₀) и D= X₁- X₀. Если модуль D <= E, то X₁ является корнем уравнения. В противном случае переходим ко второй итерации: вычисляем Х₂ и новое значение D=X₂- X₁. Опять проводим проверку на точность и при необходимости продолжаем итерации. Если X₀ выбрано правильно, а функция удовлетворяет достаточному условию сходимости, то эта итерирующая процедура быстро сойдется к корню.

Пример 1.3.

Уточним методом Ньютона значение корня, отделенного в примере 1.1. Итак пусть f(X)= exp(X) - 10*X, для первого корня a=0 и b=0,5. Пусть Е=0,00001. Формулы для первой и второй производной f(X) таковы

f’(X) = exp(X) - 10 и f”(X) = exp(X).

Очевидно, что X₀ = a = 0, т.к. f(0)*f”(0) = 1 >0.

Запрограммируем метод Ньютона для этого примера на том же рабочем листе, где мы проводили отделение корней. В ячейку А42 внесем число, равное Х₀=0. В ячейку В42 запишем формулу =EXP(A42)-10*А42, в ячейку С42 формулу =EXP(A42)-10, а в ячейку D42 формулу =А42- В42/C42. Затем в ячейку Е42 запишем формулу =А42-D42. Таким образом 42 строка содержит данные по первой итерации.

Чтобы получить в строке 43 данные по второй итерации, скопируем содержимое ячейки D42 в ячейку А43, записав в А43 формулу =D42. Далее надо скопировать формулы ячеек В42, С42, D42, E42 в ячейки В43, С43, D43, E43. Для получения данных всех остальных итераций надо выделить ячейки в 43 строке и скопировать их содержимое в блок А44:Е47. После этого следует проанализировать изменение D в столбце E, найти D<0,00001 по модулю и выбрать соответствующее ему значение Х из столбца А. Это и есть приближенное значение корня. При правильно введенных формулах метод Ньютона сходится за 3 или 4 итерации. Поэтому строить диаграмму для этого метода нет необходимости.

1.4. Уточнение корней: метод бисекции (деления отрезка пополам).

Для уточнения корня методом бисекции должно быть дано:

1) уравнение f(X) = 0, причем f(X) должна быть задана в виде формулы,

2) числа a - левая граница и b - правая граница интервала, внутри которого лежит один корень,

3) число Е - заданная точность получения корня.

Напомним, что на концах интервала функция f(X) имеет разные знаки. Вычислительная процедура метода состоит в том, что на каждом шаге итерации на интервале [a,b] выбирают промежуточную точку с так, чтобы она являлясь серединой интервала, т.ет с=(a+b)/2. Тогда интервал разделится этой точкой на два равных отрезка [a,c] и [c,b], длины которых равны (b-a)/2. Из двух полученных отрезков выберем тот, на концах которого функция f(X) принимает значения противоположных знаков. Обозначим его снова как [a,b]. На этом заканчивается первая итерация. Далее новый отрезок [a,b] делим снова пополам и проводим вторую и последующие итерации. Процесс деления отрезка пополам производим до тех пор, пока на каком-либо К-том шаге вновь получающийся отрезок не станет меньше или равным величине точности Е. Значение шага К легко рассчитать из формулы

(b-a)/2^k <=E,

где a и b - начальные значения левой и правой границ интервала.

Метод бисекций сходится для любых непрерывных функций, в том числе и недифференцируемых.

Пример 1.4.

Уточним методом бисекции значение корня, отделенного в примере 1.1. Итак пусть f(X)= exp(X) - 10*X, для первого корня a=0 и b=0,5. Пусть Е=0,00001.

Запрограммируем метод бисекции для этого примера на том же рабочем листе, где мы проводили отделение корней. В ячейки А52 и В52 надо внести числовые значения a и b,в ячейку С52 - формулу =(А52+В52)/2. Далее в ячейку D52 внесем формулу =EXP(A52)-10*A52, в ячейку Е52 - формулу =EXP(C52)-10*C52, в ячейку F52 - формулу =D52*E52, и, наконец, в ячейку G52 запишем формулу =B52- A52. В строке 52 мы сформировали первую итерацию. На второй итерации значения в ячейках А53 и В53 зависят от знака числа в ячейке F52. Если F52>0, то значение А53 равно С52. В противном случае оно должно быть равно А52. В ячейке В53 наоборот: если F52<0, то значение В53 равно С52, иначе В52.

Разрешить это затруднение поможет встроенная функция EXCEL, которая носит название ЕСЛИ. Сделаем текущей ячейку А53. В строке формул, рядом с зеленой галочкой щелкнем на кнопке с изображением f(x). Так вызывается Мастер Функций. В появившемся диалоге выберем в поле Категории Функции категорию Логические, а в поле Имя Функции - имя ЕСЛИ. На втором шаге диалога заполним три свободных поля следующим образом: в поле Логическое_выражение внесем “F52>0” (разумеется без кавычек!), в поле Значение_если_истина внесем С52, а в поле Значение_если_ложь - А52. Щелкнем по кнопке Закончить. Вот и все.

То же самое надо проделать с ячейкой В53. Только Логическое выражение будет “F52<0”, Значение_если_истина будет С52, а Значение_если_ложь соответственно В52.

Далее надо скопировать формулы в блоке ячеек С52:G52 в блок С53:G53. После этого вторая итерация будет проведена в строке 53. Для получения следующих итераций достаточно скопировать формулы из строки 53 в блоке А53:E53 в блок А54: E68. Затем, как обычно, следует найти с столбце Е такую строку, где значение D будет меньше Е. Тогда число в столбце С в этой строке и есть приближенное значение корня.

Можно построить диаграмму изменения значений в столбцах А, В и С, начиная с первой и кончая последней итерацией. Для этого нужно выделить блок ячеек А52:С68. За дальнейшими инструкциями обратитесь к примеру 1.2.

1.5 Уточнение коней: подпрограмма EXCEL “ Подбор параметра ”.

EXCEL обладает большим набором средств, позволяющих решить те или иные вычислительные задачи. Для решения нелинейного уравнения предусмотрена подпрограмма Подбор параметра. Продемонстрируем действие этой подпрограммы на предыдущем примере.

Уточним значение корня, отделенного в примере 1.1. Итак пусть f(X)= exp(X) - 10*X. Найдем корень, лежащий на интервале [0; 0,5]. Оставим пустой ячейку А70. В ячейку В70 запишем формулу =EXP(A70)-10*A70. Выберем команду меню Сервис - Подбор параметра. Откроется диалог Подбор параметра, в котором в поле Установить в ячейке запишем В70, в поле Значение занесем 0 (ноль), в поле Изменяя ячейку укажем А70. Щелкнем по кнопке ОК и появится новый диалог, в котором будет показан результат выполнения операции. В окне Состояние подбора решения будет показано найденное значение. Теперь если щелкнуть на кнопке ОК, в ячейку А70 будет внесено найденное значение корня, а в ячейку B70 - значение функции.

Для того, чтобы найти другой корень, лежащий на интервале [3,5; 4] необходимо изменить начальное приближение, которое в нашей таблице находится в ячейке A70. Запишем в эту ячейку одну из границ интервала, например, 4, и снова выполним процедуру подбора параметра. Содержимое клеток A70 и B70 изменится, теперь в этих клетках появятся координаты большего корня.

2. СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ.

В общем виде система линейных алгебраических уравнений записывается так: a₁₁x₁+a₁₂x₂+... +a_1nx_n = b₁

a₂₁x₁+a₂₂x₂+... +a_2nx_n = b₂

....................................................

a_n1x_n+a_n2x₂+... +a_nnx_n = b_n

Совокупность коэффициентов этой системы запишем в виде квадратной матрицы A из n строк и n столбцов

a₁₁ a₁₂ ... a_1n

a₂₁ a₂₂... a_2n

...........................

a_n1 a_n2... a_nn

Используя матричное исчисление, исходную систему уравнений можно записать в виде

А*Х = В,

где Х - вектор- столбец неизвестных размерностью n, а В - вектор- столбец свободных членов, тоже размерностью n.

Эта система называется совместной, если она имеет хотя бы одно решение, и определенной, если она имеет одно единственное решение. Если все свободные члены равны нулю, то система носит название однородной.

Необходимым и достаточным условием существования единственного решения системы является условие DET=0, где DET - определитель матрицы А. На практике при вычислениях на компьютере не всегда удается получить точное равенство DET нулю. В том случае, когда DET близко к нулю, системы называются плохо обусловленными. При их решении на компьютере малые погрешности в исходных данных могут привести к существенным погрешностям в решении. Условие DET~0 является необходимым для плохой обусловленности системы, но не достаточным. Поэтому при решении системы на ЭВМ требуется оценка погрешности, связанной с ограниченностью разрядной сетки компьютера.