VecMatrix( int Ans[], int Vec[], int *M, int hi1, int hi2 )
/* Ans[0..hi2) = Vec[0..hi1) x M[0..h1)[0..hi2)   -- L.Allison */
 { int j, k, Ansj;

   for(j=0; j < hi1; j++)
    { Ansj=0;
      for(k=0; k < hi2; k++)
      /* Ansj += Vec[k] x M[k][j] */
	 Ansj += Vec[k] * M[k*hi2+j];
      Ans[j]=Ansj;
    }
 }/*VecMatrix*/

/* Multiply a Vector and a Matrix */