Question

此方程组用matlab或者Mathematica解，代码该如何写

求助下这个方程组用matlab解，代码该如何写，得出的a,b,c,d,e,f,g,h整数解（温馨提示不只是正整数）是多少。（希望能穷举所有解并附上代码运行得出结果的截图）

%a+b=13
%c-d*e=4
%f+g-h=4
%a+c/f=4
%b-d*g=4
%e+h=5
用matlab或者Mathematica解都可以，能看到完整代码运行成功的截图就行

Answer 1

syms a b c d e f g h
[a,b,c,d,e,f,g,h]=solve(a+b==13,c-d*e==4,f+g-h==4,a+c/f==4,b-d*g==4,e+h==5,[a,b,c,d,e,f,g,h]);
disp([a,b,c,d,e,f,g,h])
结果：
[ 9, 4, -20, -24/5, 5, 4, 0, 0]

追问

说了要求是整数解~ -24/5是整数么？望兄台改进代码，且记得附个图。

追答

>> syms a b c d e f g h
>> assume(a,'integer')
>> assume(b,'integer')
>> assume(c,'integer')
>> assume(d,'integer')
>> assume(e,'integer')
>> assume(f,'integer')
>> assume(g,'integer')
>> assume(h,'integer')
>> [a,b,c,d,e,f,g,h]=solve(a+b==13,c-d*e==4,f+g-h==4,a+c/f==4,b-d*g==4,e+h==5,[a,b,c,d,e,f,g,h]);
disp([a,b,c,d,e,f,g,h])
Warning: Solutions might be lost.
> In symengine
In mupadengine/evalin (line 102)
In mupadengine/feval (line 158)
In solve (line 300)
Warning: Cannot find explicit solution.
> In solve (line 316)
Empty sym: 0-by-8
说明，没有全部是整数的解。

追问

我已验证的一组整数解：
a: 4
b:9
c: 0
d: -1
e: 4
f: 10
g: -5
h: 1，敬请完善代码

追答

6个方程，8个未知数，这是个不定方程组。