format(’long’)
A=[ 0.119000000000000 – 3.470000000000000 0.000000178000000
– 0.524000000000000 0.005350000000000;
438000.000000000000000 0.866000000000000 0.000000325000000
– 0.005660000000000 – 0.000001940000000;
5600.000000000000000 0.035000000000000
– 0.000007740000000 0.021800000000000 14100.000000000000000;
458000.000000000000000 – 0.860000000000000 0.000000008000000
0.570000000000000 – 0.000715000000000;
– 16900.000000000000000 – 6120.000000000000000
– 0.000005120000000 – 261000.000000000000000
– 0.000000645000000;]
cond (A)
B=[– 3.869649822000000; 438000.860338384980000;
19700.056792259998000; 457999.709285008020000;
– 284020.000005765000000]
x=A
\
B
Результат решения тестовой задачи в пакете MATLAB:
x =
1.000000000000000
1.000000000003840
1.000074582711069
0.999999999999908
1.000000000000041
Для тестовой задачи Вохминцева решатель СЛАУ в пакете
MATLAB не обеспечил даже шести верных значащих цифр для тре-
тьей неизвестной! Аналогичный результат решения тестовой задачи
Вохминцева получился для решателя СЛАУ в пакете Maple и в пакете
Mathcad. Поскольку пакет Mathcad использует программы из библио-
теки стандартных программ Intel MKL на языке Си, то этот результат
относится и к этой библиотеке. Решатель СЛАУ из математической
библиотеки стандартных программ IMSL на языке Си также не обес-
печил требуемой точности при решении тестовой задачи Вохминцева.
Аналогичный результат решения тестовой задачи Вохминцева полу-
чился для решателей СЛАУ из математических библиотек стандарт-
ных программ CUBLAS+MAGMA. Даже решатель СЛАУ в пакете
Mathematica для тестовой задачи Вохминцева не обеспечил требуемой
точности в 15 верных значащих цифр для всех элементов вектора
решений СЛАУ, но выдал предупреждение о возможных неверных
значащих цифрах в решении. Решение тестовой задачи Вохминцева
с точностью более чем 15 верных значащих цифр для всех неизвест-
ных СЛАУ было получено только при использовании вычислений
32 ISSN 0236-3933. Вестник МГТУ им. Н.Э. Баумана. Сер. “Приборостроение”. 2011. № 4
