Оганичения на линейные уравнения, Как решать? Опции

[RHnd]

Имеем невырожденную матрицу A 6x6. Произвольно назначаем вектор Y 6x1 так, чтоб корни полинома, образуемого этим вектором, имели действительную часть меньше r0. Требуется найти X: X*A=Y. Решал эту задачу в матлабе просто: R=[желаемые корни]; Y=poly®; X=inv(A)*Y;
Теперь появилось ограничение на вектор X - первые три компоненты (x(1), x(2), x(3)) должны быть одного знака, например, больше нуля - не принципиально. В идеале, все компоненты вектора X должны быть одного знака.
Как решать?

[RHnd]

Вариант с рэндомом был попробован в первую очередь. Получилось, что для real(Ri)<-1 решения находится быстро, а вот для real(Ri)<-10 - решение так и не нашлось.
Кстати, похожая ситуация и просто с нахождением полинома. Т.е. генерируем рандомно коеффициенты полинома, ищем корни. Сгенерировать такой полином, чтоб все корни меньше -10 - не удалось.

Соответственно это задача при строгой постановке - оптимизационная задача программирования нелинейного размерности 16. Тут не до симплекса :-)

Так и есть. Только там не 5 пар корней - комплексные должны быть комплексно сопряжены, коеффициенты полинома - действительные числа.
И как подступаться к такой задаче?