The linear homotopy theory for codifferential operator on Riemannian manifolds is developed in analogy to a similar idea for exterior derivative. The main object is the cohomotopy operator, which singles out a module of anticoexact forms from the module of differential forms defined on a star-shaped open subset of a manifold. It is shown that there is a direct sum decomposition of a differential form into coexact and anticoexat parts. This decomposition gives a new way of solving exterior differential systems. The method is applied to equations of fundamental physics, including vacuum Dirac-Kähler equation, coupled Maxwell-Kalb-Ramond system of equations occurring in a bosonic string theory and its reduction to the Dirac equation.