1:

procedure rotate3x3( a: array of array of  double; i,j,k,l: integer; tau,s: double);					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:

procedure rotate3x3( a: array of array of  double; i,j,k,l: integer; tau,s: double);

    var g,h: double;

  begin

  g:=a[i][j];
  h=a[k][l];
  a[i][j]=g-s*(h+g*tau);
  a[k][l]=h+s*(g-h*tau);
  end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:

procedure rotate3x3( a: array of   double; i,j,k,l: integer; tau,s: double);


    var g,h: double;

  begin

  g:=a[i,j];
  h:=a[k,l];
  a[i,j]:=g-s*(h+g*tau);
  a[k,l]:=h+s*(g-h*tau);
  end;					

1:
2:
3:
4:
5:
6:

type
  TmyArray = array of array of double;

{...}

procedure rotate3x3(a : TmyArray; ...)