СЛАУ методом простой итерации - ошибка переполнения

Дано задание - написать программу, реализующую метод простой итерации для решения СЛАУ. Программу в общем-то написал, проект прилагаю. Проблема в том, что в некоторых случаях при попытке расчёта возникает ошибка "Floating point overflow". Привожу пример. Ввожу коэффициенты, которые мне даны в качестве примера по заданию (правый крайний столбец - то, что во втором StringGrid:

-0,77	-0.04	0,21	-0,18	-1,24
0,45	-1,23	0,06	0	0,88
0,26	0,34	-1,11	0	-0,62
0,05	-0,26	0,34	-1,12	1,17

Начальное приближение - беру 1, точность 0,0001. Программа выдаёт следующий результат: 2,0482; 0,0858; 1,0646; -0,6499. Проверяю в MathCAD - сходится.
Но вот берём какой-то другой набор, например, системку из 3 уравнений:

5	-2	11	4
3	-1	-2	-5,5
7	3	2	16

MathCAD даёт ответ: 0,222; 4,138; 1,015.
Моя программка выкидывает ошибку и выводит в memo неверный результат.
Я так понял, что массив x0 (см. исходники) после некоторого числа итераций заполняется в итоге этими значениями, большими по модулю. Почему так происходит (и не всегда, а время от времени), я, честно говоря, пока не понял. Возможно, что-то надо поправить в алгоритме. Но с другой стороны, программа же смогла расчитать первое уравнение (см. табл.) и не только его. Короче, я в задумчивости.

Уважаемые форумчане, может быть, кто-нибудь сможет поправить написанное?

Вложения

slau.zip (11.2 Кб, 180 просмотров)

0

Вопрос снят, разобрался сам

0