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

if (i==j)

c=0;

for k=1:(i-1)

c=c+t(k,i)^2;

end

t(i,i)=sqrt(a(i,i)-c);

else

if (i<j)

c=0;

for k=1:(i-1)

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

end

t(i,j)=(a(i,j)-c)/t(i,i);

end

end

end

end

y=zeros(n,1);

y(1)=f(1)/t(1,1);

for i=2:n

c=0;

for k=1:(i-1)

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

end

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

end

x=zeros(n,1);

x(n)=y(n)/t(n,n);

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