Well, Frobenius, offhand, applies to a differential system generated by a family of $1$-forms. But let's be clever. In standard coordinates $(y^1,\dots,y^m)$, we can write $$\omega = \iota_X dy^1\wedge\dots\wedge dy^m$$
for some vector field $X$ with $X(p)\ne 0$. By the flowbox theorem, we can choose local coordinates $(z^1,z^2,\dots,z^m)$ so that $X = \partial/\partial z^1$, so $\omega = f\,dz^2\wedge\dots\wedge dz^m$ for some nonzero function $f$. Since $d\omega = 0$, we infer that $\dfrac{\partial f}{\partial z^1} = 0$, and so $f=f(z^2,\dots,z^m)$. Now set
$$x^2=\int f(z^2,\dots,z^m)\,dz^2\,,$$ so that
$$dx^2 = f\,dz^2 + \sum_{j=3}^m \left(\int \dfrac{\partial f}{\partial z^j}\,dz^2\right)\,dz^j\,.$$
It follows that if we set $x^1=z^1, x^3=z^3,\dots, x^m=z^m$, then $\omega = dx^2\wedge\dots\wedge dx^m$.