code 31507;
   procedure LNGFULSYMMATVEC(LR, UR, LC, UC, A, B, C);
   value LR, UR, LC, UC, B; integer LR, UR, LC, UC;
   array A, B, C;
   begin real D, DD;
      for LR:= LR step 1 until UR do 
      begin LNGSYMMATVEC(LC, UC, LR, A, B, 0, 0, D, DD);
         C[LR]:= D + DD
      end 
   end LNGFULSYMMATVEC;
        eop