Eu preciso fazer um modelo de programação linear. Aqui estão as desigualdades que estou usando (por exemplo):

6x + 4y <= 24
x + 2y <= 6
-x + y <= 1
y <= 2

Preciso encontrar a área descrita por essas desigualdades e sombrear em um gráfico, bem como acompanhar os vértices das linhas delimitadoras dessa área e desenhar a linha delimitadora em uma cor diferente. Veja o gráfico abaixo para um exemplo do que estou procurando.

.

Estou usando o Python 3.2, numpy e matplotlib. Existem módulos melhores para programação linear em Python?