python求解一元一次方程

待求解方程： 61.08x/(61.08x+18.02*(1-x-0.15*x)) = 0.3

from sympy import *


x = symbols('x')
z = solve(61.08*x/(61.08*x+18.02*(1-x-0.15*x)) - 0.3, x)
print(z)
[0.110387581703350]

from sympy import *

''' 待求解方程：61.08x/(61.08x+18.02*(1-x-0.15*x)) = 0.3
其中0.3表示mea的质量分数为0.3
0.15*x中的0.15表示co2的负载为0.15mol/mol。y/x = 0.15
x表示mea的摩尔分数
'''
q = 0.15  #co2的负载
mf = 0.3  #mea的质量浓度
x = symbols('x')
z = solve(61.08*x/(61.08*x+18.02*(1-x-q*x)) - mf, x)

'''保留4位小数'''
mea = round(z[0],4)
co2 = round(mea*q,4)
h2o = round(1-mea-co2,4)
print('mea的摩尔分数为：',mea)
print('co2的摩尔分数为：',co2)
print('h2o的摩尔分数为：',h2o)