Контрольная работа: Метод квадратных корней

sm=0

for k=(i+1):s

sm=sm+T(i,k)*X(k)

sm

end

X(i)=(Y(i)-sm)/T(i,i)

E=A*X-B'

end

else

error('B не соответствует матрице А')

end

else

error('det А = 0')

end

else

B = B*A'

A = A*A'

if det(A)~=0

s=size(A,1);

if size(B',1) == s

T=zeros(s);

T(1,1)=sqrt(A(1,1));

for k=2:s

T(1,k)=A(1,k)/T(1,1)

end

for j=2:s

for i=2:s

if i==j

sm=0