We construct examples of robust homoclinic orbits for systems of ordinary differential equations coupled with the Preisach hysteresis operator. Existence of such orbits is demonstrated for the first time. We discuss a generic mechanism that creates robust homoclinic orbits and a method for finding them. An example of a homoclinic orbit in a population dynamics model with hysteretic response of the prey to variations of the predator is studied numerically.
