за программировать

1.  Ввод размерности n; m=0, Enorm21=Enorm22=0.0;
2.  i=1, 
3.  j=1
4.  Ввод A[i,j]
5.  если  (i<>j)&abs(A[i,j])>max_el , то max_el:=abs(A[i,j]); p:=i; q:=j;
6.   Enorm21=Enorm21+A[i,j]*A[i,j];
7.  j=j+1
8.  если j<=n пер. на п.4
9.  i=i+1
10.  если i<=n, то пер. на п.3
11. L:=Arctan( (2*A[p,q])/(A[p,p]-A[q,q]) )/2
12.  T1=a[p,p]*cos(l)*cos(l)+2*a[p,q]*sin(l)*cos(l)+a[q,q]*sin(l)*sin(l);
13.  T2=a[p,p]*sin(l)*sin(l)-2*a[p,q]*sin(l)*cos(l)+a[q,q]*cos(l)*cos(l);
14.  A[p,q]:=(a[q,q]-a[p,p])*cos(l)*sin(l)+a[p,q]*(cos(l)*cos(l)-sin(l)*sin(l));
15.  A[q,p]=A[p,q]; a[p,p]=T1; a[q,q]=T2
16.  i=1
17. Если (i<>p) and (i<>q), то пер.на п.22
18. T1=a[i,p]*cos(l)+a[i,q]*sin(l);
19.  A[i,q]:=a[i,p]*sin(l)-a[i,q]*cos(l);
20. A[i,p]=T1; A[p,i]:=A[i,p];
21.  A[q,i]:=A[i,q];
22.  i=i+1
23.  если i<=n, то пер. на п.17
24. max_el=A[2,1]
25. p=2; q=1
26. i=2
27. j=1
28. если  A[i,j]>max_el , то max_el:=A[i,j]; p:=i; q:=j;
29.  j=j+1
30. если j<=i, то  пер. на п.28
31. i=i+1
32.  если i<=n, то пер. на п.27
33. m=m+1
34. Если max_el>1.0E-12, то пер. на п.11
35. i:=1
36. Вывод A[i,i] - на экран
37. Enorm22=Enorm22+ A[i,i]* A[i,i]
38.  i=i+1
39.  если i<=n пер. на п.36
40. Вывод abs(Enorm21-Enorm22)
41. Вывод m – количество итераций