Hola de nuevo: Volviendo con el producto de matrices, para ver si funcionaba la solución que amablemente me envió Carlos
Amelia:
con la función load(eigen)  se llama a la funcion inprod(a,b) que hace lo que te interesa fila por columna.

me encontré que sigue produciéndose el mismo problema cuando se utiliza después de haber cargado el paquete "vect"; pienso que ahí está el problema.
Adjunto un documento en el que se muestra que las cosas van bien si no está cargado "vect", pero en el momento en que se carga "vect", las cosa dejan de funcionar para el producto de matrices.
Un saludo
Amelia