Добавил:
korayakov
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:ЧМ (ЭКТ-3) / Лабы / Всякие сделанные варианты / other / laba3p2
.docA=rand(9)
C=imag(eig(A))
k=0
for i=1: size(A,1)
if C(i)~=0
k=k+1
break
end
end
while( k~=0)
A=rand(9)
C=imag(eig(A))
k=0
for i=1: size(A,1)
if C(i)~=0
k=k+1
break
end
end
end
[Q,R]=qr(A)
for i=1:20
B=R*Q
[Q,R]=qr(B)
end