recently in my researches I've come across the following operator

$$L\left(\begin{array}{c} a_1\\ \vdots\\ a_n \end{array}\right)=M_1\left(\begin{array}{c} \partial_{z_1}a_1\\ \vdots\\ \partial_{z_1}a_n \end{array}\right)+\dotso+M_n\left(\begin{array}{c} \partial_{z_n}a_1\\ \vdots\\ \partial_{z_n}a_n \end{array}\right),$$ where $a_1,\dotso,a_n:\mathbb R^n\to\mathbb R$ are smooth functions of the independent variables $z_1,\dotso,z_n$ and $M_1,\dotso M_n$ are real square $n\times n$ matrices with constant coefficients (though none of them is invertible).

Is there a general theory to treat this simple system of linear first order pdes? Is it possible to write explicitly a solution if we were to solve

$$L\left(\begin{array}{c} a_1\\ \vdots\\ a_n \end{array}\right)=\left(\begin{array}{c} b_1\\ \vdots\\ b_n \end{array}\right),$$ for some smooth functions $b_1,\dotso,b_n:\mathbb R^n\to\mathbb R$?

Thank you for the references and the patience.

Kind regards,

Guido

**Edit**

Following the suggestions, I report here the $3\times 3$ system which is the toy model I'm currently studying

$$\left(\begin{array}{ccc} B_{11} & B_{12} & 0\\ -B_{11} & 0 & 0\\ 0 & -B_{12} & 0 \end{array} \right)\left(\begin{array}{c} \partial_{z_1}a_1\\ \partial_{z_1}a_2\\ \partial_{z_1}a_3\\ \end{array}\right)+ \left(\begin{array}{ccc} 0 & 0 & -B_{23}\\ 0 & B_{22} & B_{23}\\ 0 & -B_{22} & 0 \end{array} \right)\left(\begin{array}{c} \partial_{z_2}a_1\\ \partial_{z_2}a_2\\ \partial_{z_2}a_3\\ \end{array}\right)+ \left(\begin{array}{ccc} 0 & 0 & -B_{33}\\ -B_{31} & 0 & 0\\ B_{31} & 0 & B_{33} \end{array} \right)\left(\begin{array}{c} \partial_{z_3}a_1\\ \partial_{z_3}a_2\\ \partial_{z_3}a_3\\ \end{array}\right)= \left(\begin{array}{c} b_1\\ b_2\\ -b_1-b_2 \end{array}\right), $$

$a_i$ and $b_i$ are smooth functions, the $b_i$'s are compactly supported. I'm interested in solving locally the above system. I think it has something to do with the usual theory for the hyperbolic systems of differential equations, however I'm to unexperienced in the field to be sure about the steps to follow. I understand the condition to be imposed on the functions $b_i$ as a necessary condition for the system to be possibly solvable, but it is not clear to me why there should exist a solution at all! This is the most I can say up to now. Thanks again for the patience and the help