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