Собственно вопрос по работе в пакете 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. Если есть решения в других матпакетах, можно и так =)