Решение в wxMaxima

[IGeorgeI]

Собственно вопрос по работе в пакете Maxima, ибо не нашел более подходящего раздела.
Есть некая несовместная система линейных уравнений. Необходимо найти набор значений переменных при котором отклонение полученных левых частей от правых будет минимально. Собственно сложность возникла на почве того что переменных по условию могут быть только неотрицательными. Кто-нибудь может подсказать на каком этапе и как забить это условие.
P.S. Пробовал через assume, он его игнорирует.
Код под спойлером.

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

eq: matrix(
[0.65, 0.23, 17, 84, 0.94, 1.04, 0.93, 11.6],
[0.85, 0.25, 0.12, 0.9, 0.8, 0.9, 0.87, 11.6],
[95, 25, 13, 171, 99, 405, 175, 1675],
[257, 75, 9, 207, 0, 129, 68.7, 2670],
[10, 6, 40, 20, 543, 62.6, 50, 930],
[25, 10, 1, 31, 0, 77, 32, 335],
[8.3, 1.4, 0.4, 4, 3.2, 5.9, 3, 95],
[2, 0.4, 0.5, 2.9, 0.2, 12.9, 6, 55]);
teq:transpose(eq);
e1:teq[1];
e2:teq[2];
e3:teq[3];
e4:teq[4];
e5:teq[5];
e6:teq[6];
e7:teq[7];
b:teq[8];

fin:matrix(
[e1.e1,e1.e2,e1.e3,e1.e4,e1.e5,e1.e6,e1.e7,e1.b],
[e2.e1,e2.e2,e2.e3,e2.e4,e2.e5,e2.e6,e2.e7,e2.b],
[e3.e1,e3.e2,e3.e3,e3.e4,e3.e5,e3.e6,e3.e7,e3.b],
[e4.e1,e4.e2,e4.e3,e4.e4,e4.e5,e4.e6,e4.e7,e4.b],
[e5.e1,e5.e2,e5.e3,e5.e4,e5.e5,e5.e6,e5.e7,e5.b],
[e6.e1,e6.e2,e6.e3,e6.e4,e6.e5,e6.e6,e6.e7,e6.b],
[e7.e1,e7.e2,e7.e3,e7.e4,e7.e5,e7.e6,e7.e7,e7.b]
);
X:[x1,x2,x3,x4,x5,x6,x7,-1];
A:matrix(
[e1.e1,e1.e2,e1.e3,e1.e4,e1.e5,e1.e6,e1.e7],
[e2.e1,e2.e2,e2.e3,e2.e4,e2.e5,e2.e6,e2.e7],
[e3.e1,e3.e2,e3.e3,e3.e4,e3.e5,e3.e6,e3.e7],
[e4.e1,e4.e2,e4.e3,e4.e4,e4.e5,e4.e6,e4.e7],
[e5.e1,e5.e2,e5.e3,e5.e4,e5.e5,e5.e6,e5.e7],
[e6.e1,e6.e2,e6.e3,e6.e4,e6.e5,e6.e6,e6.e7],
[e7.e1,e7.e2,e7.e3,e7.e4,e7.e5,e7.e6,e7.e7]
);

solve(
[35619.14400000001*x7+74255.211*x6+14863.251*x5+70513.36500000001*x4+3988.472*x3+21972.782*x2+75873.035*x1=863905.9,
10154.5314*x7+20959.4842*x6+5737.976199999999*x5+20256.305*x4+1254.7*x3+6388.2354*x2+21972.782*x1=251215.568,
4945.414400000001*x7+9033.598000000002*x6+23024.456*x5+6348.158*x4+2140.4244*x3+1254.7*x2+3988.472*x1=83604.092,
46246.203*x7+99746.18000000001*x6+27882.06*x5+80532.22*x4+6348.158*x3+20256.305*x2+70513.36500000001*x1=869624.3400000001,
44487.3702*x7+74109.95760000001*x6+304661.8036*x5+27882.06*x4+23024.456*x3+5737.976199999999*x2+14863.251*x1=671150.184,
85428.15019999999*x7+190716.8716*x6+74109.95760000001*x5+99746.18000000001*x4+9033.598000000002*x3+20959.4842*x2+74255.211*x1=1108110.504,
38915.3118*x7+85428.15019999999*x6+44487.3702*x5+46246.203*x4+4945.414400000001*x3+10154.5314*x2+35619.14400000001*x1=534409.88],
[x1,x2,x3,x4,x5,x6,x7]);


Теория для решения взята отсюда

P.P.S. Задача для себя, не из учебы.
P.P.P.S. Если есть решения в других матпакетах, можно и так =)