Thank you for your response. However, I am not trying to solve the system of ODEs. RAther, I am trying to modify the coefficient matrix C, more specifically, trying to use Gauss Elimination so that I can get L(y)*w(y) = 0 where L(y) is a differential operator. Do you think this is possible given the way you have set up the coefficient matrix? The final matrix would look something like this: C_modified := Matrix(3,3,[> [P(y), Q(y), R(y)],[0, A(y), B(y)],[0, 0, L(y)]]); V := [u(y);v(y);w(y)];

Thank you for your response. However, I am not trying to solve the system of ODEs. RAther, I am trying to modify the coefficient matrix C, more specifically, trying to use Gauss Elimination so that I can get L(y)*w(y) = 0 where L(y) is a differential operator. Do you think this is possible given the way you have set up the coefficient matrix? The final matrix would look something like this: C_modified := Matrix(3,3,[> [P(y), Q(y), R(y)],[0, A(y), B(y)],[0, 0, L(y)]]); V := [u(y);v(y);w(y)];