Enable an indexed summation of matrices

- on AddWithLimits.expand() distribute sum to each matrix entry
- Sum.doit() recognizes matrix summands and distributes
- Matrix.doit() maps to matrix entries
- includes 1 new test