Given a system of vector fields on a smooth manifold that spans a plane field of constant rank, we present a systematic method and an algorithm to find submanifolds that are invariant under the flows of the vector fields. We present examples of partition into invariant submanifolds, which further gives partition into orbits. We use the method of generalized Frobenius theorem by means of exterior differential systems.
