我想在没有任何模块（NumPy、Sympy...等）的情况下求解这个方程
Px + Qy = W
（例如 5x + 6y = 55）
谢谢。
【问题讨论】：
x(y) = ... 还是y(x) = ...？
标签： python-3.x
这是一种非常粗略的方法，但您可以使用蛮力技术，正如我在您的问题下的评论中所说的那样。它可能可以优化很多，只给出int 输出，但总体显示方法：
import numpy as np
# Provide the equation:
print("Provide a, b and c to evaluate in equation of form {ax + by - c = 0}")
a = float(input("a: "))
b = float(input("b: "))
c = float(input("c: "))
x_range = int(input("x-searching range (-a, a): "))
y_range = int(input("y-searching range (-b, b): "))
error = float(input("maximum accepted error from the exact solution: "))
x_range = np.arange(-x_range, x_range, 1)
y_range = np.arange(-y_range, y_range, 1)
for x in x_range:
for y in y_range:
if -error <= a * x + b * y - c <= error:
print("Got an absolute error of {} or less with numbers x = {} and y = {}.".format(error, x, y))
a = 1、b = 2、c = 3、x_range = 10、y_range = 10、error = 0.001 的示例输出：
Got an error of 0.001 or less with numbers x = -9 and y = 6.
Got an error of 0.001 or less with numbers x = -7 and y = 5.
Got an error of 0.001 or less with numbers x = -5 and y = 4.
Got an error of 0.001 or less with numbers x = -3 and y = 3.
Got an error of 0.001 or less with numbers x = -1 and y = 2.
Got an error of 0.001 or less with numbers x = 1 and y = 1.
Got an error of 0.001 or less with numbers x = 3 and y = 0.
Got an error of 0.001 or less with numbers x = 5 and y = -1.
Got an error of 0.001 or less with numbers x = 7 and y = -2.
Got an error of 0.001 or less with numbers x = 9 and y = -3.
我正在使用numpy，但不是用于求解方程本身的内置函数，只是为了创建一个数组。当然，这可以在没有它的情况下完成。
【讨论】：
使用 python 求解方程有数千种方法。
其中之一是：
def myfunc (x=None, y=None):
return ((55-6*y)/5.0) if y else ((55-5*x)/6.0)
print(myfunc(x=10)) # OUTPUT: 0.833333333333, y value for x == 10
print(myfunc(y=42)) # OUTPUT: -39.4, x value for y == 42
您只需在函数内部定义求解方程所需的步骤。
在我们的示例中，如果我们有 y 值，我们将 6*y 减去 55 然后除以 5.0（我们添加 .0 以得到浮点数），否则（意味着我们有 x）我们从55 中减去5*x，然后除以6.0
同理，可以概括：
def myfunc (x=None, y=None, P=None, Q=None, W=None):
if not W:
return P*x + Q*y
elif not x:
return (W-Q*y)/float(P)
elif not y:
return (W-P*x)/float(Q)
elif not P:
return (W-Q*y)/float(x)
elif not Q:
return (W-P*x)/float(y)
print(myfunc(x=10, P=5, Q=6, W=55)) # OUTPUT: 0.833333333333, y value for x == 10
print(myfunc(y=42, P=5, Q=6, W=55)) # OUTPUT: -39.4, x value for y == 42
查看this QA 了解解决此问题的其他有趣方法
【讨论】：