Maxima Function
rowop (M, i, j, theta)
If M is a matrix, return the matrix that results from doing the row operation R_i <- R_i - theta * R_j. If M doesn't have a row i or j, signal an error.
R_i <- R_i - theta * R_j