Решение СЛАУ простой итерацией (матрица Гильберта)
На носу сессия, а до сих пор проблемы с зачётами, а именно с прогой, осталось последнюю лабу сдать, а сама лаба какая-то муторная, я и задание толком-то не понимаю, куда уж там до реализации...
Нужно написать на языке С++ решение СЛАУ методом простой итерации. Реализовать шаблон класса матрицы, представляющий возможность обращения к хранимым в ней элементам через оператор скобки () с двумя индексами, а также основные операции с ними - присваивание, умножение, сложение, вычитание, умножение на число. у конструктора два параметра размеры по горизонтали и вертикали, эти же два параметра возвращаются двумя соответствующими const-методами. Индексация элементов матрицы производится с нуля, хранение элементов осуществляется в динамическом массиве. Решаемая слау это слау с матрицей Гильберта ( все значения меньше единицы) :
((вектор)u)-1/10H((вектор)u)=((вектор)f) , h(i,j)=1/(i+j-1), f(i)= 1, если i=1 и 0 если i не равно 1...
Полученная в результате СЛАУ A(вектор(u))=(вектор(f)) решается методом простой итерации, а именно:
((вектор)u(нулевое)) = (вектор(f))
((вектор)u(n-ое)) = ((вектор)u(n-ое минус 1)) - (A((вектор)u(n-ое минус 1)) - (вектор(f)) = ( E - A )((вектор)u(n-ое минус 1)) + (вектор(f))
Выручите, пожалуйста, кто сможет...Сдать всё это нужно до среды(((
Нужно написать на языке С++ решение СЛАУ методом простой итерации. Реализовать шаблон класса матрицы, представляющий возможность обращения к хранимым в ней элементам через оператор скобки () с двумя индексами, а также основные операции с ними - присваивание, умножение, сложение, вычитание, умножение на число. у конструктора два параметра размеры по горизонтали и вертикали, эти же два параметра возвращаются двумя соответствующими const-методами. Индексация элементов матрицы производится с нуля, хранение элементов осуществляется в динамическом массиве. Решаемая слау это слау с матрицей Гильберта ( все значения меньше единицы) :
((вектор)u)-1/10H((вектор)u)=((вектор)f) , h(i,j)=1/(i+j-1), f(i)= 1, если i=1 и 0 если i не равно 1...
Полученная в результате СЛАУ A(вектор(u))=(вектор(f)) решается методом простой итерации, а именно:
((вектор)u(нулевое)) = (вектор(f))
((вектор)u(n-ое)) = ((вектор)u(n-ое минус 1)) - (A((вектор)u(n-ое минус 1)) - (вектор(f)) = ( E - A )((вектор)u(n-ое минус 1)) + (вектор(f))
Выручите, пожалуйста, кто сможет...Сдать всё это нужно до среды(((