Метод искусственного базиса

Реализовать программу на любом языке (кроме Pascal), выполняющая задачу ЛП методом искусственного базиса
F(x): 3x1 + 2x2 -> ma
5x1 + 3x2 <=
3x1 - 10x2 = 1
x1, x2 >=
Т.е. после преобразования получается такая функция
F(x): 3x1 + 2x2 + (0x3) - Mx4 -> ma
5x1 + 3x2 + x3 + (0x4) =
3x1 - 10x2 + (0x3) + x4 = 1
x1, x2, x3, x4 >= 0