Метод вращения Якоби

[Nuc]

Здравствуйте!!! Написал код, для поиска собственных значений и векторов методом вращений Якоби

restart;
with(linalg):
n:=10;
N:=5; 
j:=1.5+0.1*n;
k:=n;
l:=n;   
A:=matrix([[j,0.5*j,0,0.2*l,0],[0.5*j,j,0.3*j,0,0.1*l],[0,0.3*j,10,-0.3*j,0.5*l],[0.2*k,0,-0.3*j,j,-0.1*j],[0,0.1*k,0.5*k,-0.1*j,j]]);
eigenvectors(A);
B:=A: k:=0: flag:=0: eps:=0.001:
while (flag = 0) do 
m:=abs(B[1,2]):
iStr:=1: iCol:=2:
for i from 1 to N do 
for j from 1 to N do 
if (i < j) and (j > 2) then 
if (m < abs(B[i,j])) then 
m:=abs(B[i,j]):
iStr:=i: iCol:=j:
end if: end if: end do: end do:
if abs(m) <= eps then
flag:=1:
else 
phi:=evalf((1/2)*arctan(2*B[iStr,iCol]/(B[iStr,iStr]-B[iCol,iCol]))):
H[k]:=diag(1,1,1,1,1): 
H[k][iStr,iStr]:=cos(phi):
H[k][iStr,iCol]:=-sin(phi): 
H[k][iCol,iStr]:=sin(phi):
H[k][iCol,iCol]:=cos(phi): 
if k = 0 then 
G:=H[0]:
end if:
#print(B);
#print(H[k]);
B:=multiply(transpose(H[k]), multiply(B,H[k])): 
if k <> 0 then 
G:=multiply(G,H[k]):
end if:
k:=k+1:
end if: end do:
#print(k);
G:=multiply(seq(H[i],i=0..k-1)):
for i from 1 to N do 
lambda[i]:=B[i,i]; 
e[i]:=col(-G,i); 
end do;

но препод сказал, что он коряв с точки зрения программирования.
Сможете ответить, в чём проявляется корявость?