I know the power of Maple for exact symbolic and numeric computing. With Maple, we can calculate symbolic formulas expressed with scalar values such as x, y, z, and so on. For example, >> solve(x+z+1=0,x); x=-z-1 However, I want to calculate formulas including matrices algebraically, too. I’ll show an example of my ideal result; >> solve(A*x=b,x); x=A^-1*b or x=(A^T*A)^-1*A^T*b In the example above, I assume A is a matrix, and x, b are vectors. “A^-1” means the inverse of the matrix A and “A^T” means the transpose of the matrix A. I want to keep “A” as a matrix and get a result like the equation which is often written in well-known textbooks of linear algebra.