Предлагаю попробывать найти оценку области притяжения "изнутри" тривиального решения математической модели, описывающей движение крылатого аппарата в фазовом пространстве системы, заданной уравнениями:
[MATH]sys(diff(x;t)=y;diff(y;t)=-y+z-sign(x+y);diff(z;t)=2*y-z;3;1)[/MATH]

Для построения оценки области притяжения данной системы, существует много методов, например можно воспользоваться методом нелокального сведения Леонова. Данный частотный критерий использует информацию о поведении решений систем более низкого порядка чем исходная. Для данной системы она будет иметь вид:
[MATH]sys(diff(x1;t)=x2;diff(x2;t)=-x2*a-sign(1-x);2;1)[/MATH]