We prove a Lusin type theorem for a certain class of linear partial differential operators G(D), reducing to [1, Theorem 1] when G(D) is the gradient. Moreover, we describe the structure of the set {G(D)f = F}, under assumptions of non-integrability on F, in terms of lower dimensional rectifiability and superdensity. Applications to Maxwell type system and to multivariable Cauchy–Riemann system are provided.