Графический метод решения ЗЛП основан на построении выпуклого многоугольника (симплекса), полностью лежащего в первом квадранте. Начальным приближением является весь первый квадрант, т.е. треугольник с вершинами (0,0), (Inf,0), (0, Inf), где Inf - бесконечность, которую нужно "отрезать".
Идеальный вариант для "отрезания" - неравенство вида ax1 + bx2 <= c, (a>0, b>0, c>0), которое ликвидирует обе сингулярных вершины. В данном случае такого неравенства нет, но есть первое неравенство (a1>0, c1>0), и второе неравенство, которое после умножения на "-1" даёт (b2>0, c2>0), причём точка пересечения M(x,y) соответствующих им прямых лежит в первом квадранте. Это значит, что получен симплекс с вершинами (0,0), (c1/a1,0), (x,y), (0,c2/a2). Вершины следуют в циклическом порядке против часовой стрелки, в таком же порядке их следует хранить.
Каждое новое неравенство будет отрезать на симплексе группу смежных точек, которые ему не удовлетворяют, и заменять их на пару новых. Алгоритм несложен, если к каждой вершине привязать уравнения прямых, соединяющих её с предыдущей и последующей, и сохранять циклический порядок вершин.
Если известны все вершины симплекса, то максимум целевой функции достигается в одной из них. Начальную точку можно выбрать из соображений удобства реализации, направление обхода - в сторону возрастания целевой функции. При любом направлении обхода можно брать первый локальный максимум.