Frobenius theorem states that a set of differential operators is completely integrable if corresponding smooth distribution is involutive. We know, however, the distribution gives a set of linear first-order partial differential equations for each differential operator acting on one single function. So my question is, whether the theorem can apply to vector-valued function and how if it does? For example, how to solve {du_i/dx_j=du_j/dx_i for i,j=1..n}?
http://en.wikipedia.org/wiki/Frobenius_theorem_(differential_topology)