Курсовая работа: Метод квадратных корней для симметричной матрицы при решении систем линейных алгебраических уравнений

end

x=zeros(n,1); %создаем столбец точных решений

e=zeros(n,1); % создаем столбец ошибок

x(n)=y(n)/t(n,n); %1.8 вычисляем вектор Х

for i=(n-1):-1:1

c=0;

for k=(i+1):n

c=c+t(i,k)*x(k);

end

x(i)=(y(i)-c)/t(i,i);

e=a*x-f';

end

else

error('Внимание! Ошибка! Размерность матрицы А не соответствует размерности вектора F');

end

else

error('Внимание! Ошибка! Определитель матрицы А равен 0')

end

else

f=f*a';

a=a*a';

if (det(a)~=0) % проверяем, чтобы система имела единственное решение

ifsize(f',1)==n%проверяем соответствует ли мерность матрицы А мерности вектора f

t=zeros(n); %создаем матрицу элементов T и заполняем ее нулями

t(1,1)=sqrt(a(1,1)); % 1.3

for k=2:n

t(1,k)=a(1,k)/t(1,1);

end

for j=2:n