code 34018;
   real procedure SYMMATVEC(L, U, I, A, B); value L,U,I;
   integer L,U,I; array A,B;
   begin integer K, M;
       M:= if L > I then L else I; K:= M * (M - 1) // 2;
       SYMMATVEC:= VECVEC(L, if I <= U then I-1 else U, K, B, A)
                   + SEQVEC(M, U, K + I, 0, A, B)
   end SYMMATVEC;
       eop