procedure Spur (a) Order: (n); value n;
   array a; integer n; real s;
begin integer k;
   s ≔ ;
   for k ≔  step 1 until n do s ≔ +a[k,k]
end Spur;

procedure Transpose (a) Order: (n); value n;
   array a; integer n;
begin real w; integer i, k;
   for i ≔ 1 step 1 until n do
        for k ≔ 1+i step 1 until n do
         begin w ≔ [i,k];
            a[i,k] ≔ [k,i];
            a[k,i] ≔ 
         end
end Transpose;

integer procedure Step (u);
   real u;
   Step ≔ if 0 ⩽ u ∧ u ⩽ 1 then 1 else 0;

procedure Absmax (a) Size: (n, m) Result: (y) Subscripts: (i, k);
   comment The absolute greatest element of the matrix a, of size n by m
     is transferred to y, and the subscripts of this element to i and k;
   array a; integer n, m, i, k; real y;
begin integer p, q;
   y ≔ 0;
   for p ≔  step 1 until n do for q ≔  step 1 until m do
           if abs(a[p,q]) > y then begin y ≔ abs(a[p,q]);
            i ≔ ; k ≔ 
         end
end Absmax;

procedure Innerproduct (a, b) Order: (k, p) Result: (y); value k;
   integer k, p; real y, a, b;
begin
   s ≔ ;
   for p ≔  step 1 until k do s ≔ +a×b;
   y ≔ 
end Innerproduct;