Hello, Is the following desired behavior? (%i1) m: matrix([1,2,3]); (%o1) [ 1 2 3 ] (%i2) m . transpose(m); (%o2) 14 One would expect multiplying a 1xn matrix by a nx1 matrix to return a 1x1 matrix rather than just the value. --Mike