code 34011;
   real procedure MATVEC(L, U, I, A, B); value L,U,I;
   integer L,U,I; array A,B;
   begin integer K; real S;
       S:= 0;
       for K:=L step 1 until U do S:= A[I,K] * B[K] + S;
       MATVEC:= S
   end MATVEC;
       eop