code 34418;
   procedure LNGSYMMATVEC(L, U, I, A, B, C, CC, D, DD);
   value L, U, I, C, CC;
   integer L, U, I; real C, CC, D, DD; array A, B;
   begin integer K, M;
      M:= if L > I then L else I; K:= M * (M - 1) // 2;
      LNGVECVEC(L, if I <= U then I - 1 else U,
      K, B, A, C, CC, C, CC);
      LNGSEQVEC(M, U, K + I, 0, A, B, C, CC, D, DD)
   end LNGSYMMATVEC;
        eop