In this paper, we study the degenerate parabolic systemwhere X = {X 1 , . . . , X m } is a system of smooth real vector fields satisfying Hörmander's condition and the coefficients a αβ ij are measurable functions and their skew-symmetric part can be unbounded. After proving the L 2 estimates for the weak solutions, the higher integrability is proved by establishing a reverse Hölder inequality for weak solutions.
MSC: Primary 35K65; secondary 35K40; 35B65