怎么解方程···········？

# w=1800
# a+b=w
# a>b
# a:w=b:a约等于0.618
# 求a的值
#导入sbmpb
from sympy import *

#定义变量
a = Symbol('a')
b = Symbol('b')
print(solve([a+b-1800,a/1800-0.618,b/a-0.618],[a,b]))

由于浮点数计算的精度问题会导致方程无解

二元一次方程式的求解只需要两条方程，例如a:w(1800)=0.618,b:a=0.618可以确定a与b的值

我们将其代入代码进行计算

>>> from sympy import *
>>> a = Symbol('a')
>>> b = Symbol('b')
>>> print(solve([a/1800-618/1000,b/a-618/1000],[a,b]))
{a: 1112.40000000000, b: 687.463200000000}

此时已经得出唯一解a=1112.4,b=687.4632,并且a+b=1799.863,与a+b=1800有一定的误差；故将此唯一解代回a+b=1800方程时会否定该解，造成无解的现象。

所以说只需要删除三条方程中的任意一条即可。

？？你不是解了吗

因题目要求中b/a约等于0.618.

from sympy import *

a = Symbol('a')

b = Symbol('b')

sol = solve([a/1800-0.618, a+b-1800], [a, b])

print(sol)

#{a: 1112.40000000000, b: 687.600000000000}

#验证

x,y=sol[a],sol[b]

print(y/x)

#0.618122977346278，符合条件要求。

哈哈