read(m);
if m=1 then Assign(f,'Matrix.txt') else Assign(f,'Matrix2.txt'); writeln('pogreshost'); Read(eps); Reset(f); for i:=1 to n do Begin for j:=1 to n do Begin read(f,a[i,j]); a1[i,j]:=a[i,j]; End; End; norm1:=norma(a); DECOMP(a,ip,jp,eps,det,norm1,ifsolve); if (ifsolve and (abs(det)>eps)) then Begin for i:=1 to n do b[i]:=0; for k:=1 to n do Begin b[k]:=1; SOLVE(a,ip,jp,x,b); b[k]:=0; for i:=1 to n do Ao[i,k]:=x[i]; End; for i:=1 to n do Begin for j:=1 to n do Begin s:=0; for k:=1 to n do s:=s+Ao[i,k]*a1[k,j]; if i=j then s:=abs(1-s) else s:=abs(s-0); write(s:10); End; Writeln; End; writeln('obratnaia matrica'); for i:=1 to n do Begin for j:=1 to n do write(Ao[i,j]:10:4); Writeln; End; norm2:=norma(Ao); cond:=norm1*norm2; writeln('cond(A)=',cond:9); End else writeln('systema ne imeet ed. reshenia '); Readln; Readln; End.
Результаты машинных вычислений
Вывод В ходе выполнения лабораторной работы изучены и освоены численные методы линейной алгебры: прямые методы обращения матриц. Разработана программа на языке Pascal, использующая 2 стандартные подпрограммы, и реализующая метод обращения матриц описанный в пункте 2.4. Результаты работы программы, а так же вычислительная ошибка представлены в пункте 6,(CЛАУ является хорошо обусловленной, поскольку Сond(A)
|
).
