We propose a new method for constructing integrable systems of natural form. In this method, integrals of motion are solutions of an overdetermined system of algebraic and partial differential equations obtained from the compatibility condition for Poisson tensors polynomial in the momenta and from the condition that the bi-Lagrangian distribution corresponding to the integrals of motion is invariant under the action of the recursion operator.