这是求解线性方程组的简单示例，也是对许多方程使用for循环的示例.
Here is a simple example of solving the system of linear equations and the example of using for loop for many equations.
import numpy as np from gekko import GEKKO m = GEKKO(remote=False) # Random 3x3 A = np.random.rand(3,3) # Random 3x1 b = np.random.rand(3) # Gekko array 3x1 x = m.Array(m.Var,(3)) # solve Ax = b eqn = np.dot(A,x) for i in range(3): m.Equation(eqn[i]==b[i]) m.solve(disp=False) X = [x[i].value for i in range(3)] print(X) print(b) print(np.dot(A,X))具有正确的输出.结果为X(np.dot(A，X)== b)-正确！
with the correct output. With the result X (np.dot(A,X)==b) - correct!
[[-0.45756768428], [1.0562541773], [0.10058435163]] [0.64342498 0.34894335 0.5375324 ] [[0.64342498] [0.34894335] [0.5375324 ]]在最近的Gekko 0.2rc6中，还引入了用于线性编程的axb()函数.使用此功能可能会解决相同的问题，但是我不确定如何获得正确的结果.
In the recent Gekko 0.2rc6 there is also introduced axb() function for linear programing. This might be the same problem solved with this function, but I am not sure how to get the correct result.
m = GEKKO(remote=False) # Random 3x3 A = np.random.rand(3,3) # Random 3x1 b = np.random.rand(3) # Gekko array 3x1 x = m.Array(m.Var,(3)) # solve Ax = b m.axb(A,b,x=x) m.solve(disp=False) X = [x[i].value for i in range(3)] print(X) print(b) print(np.dot(A,X))但是似乎我错过了一些事情，因为输出不是解决方案???结果X(np.dot(A，X)== b)-不正确！
but it seems I missed something because the output is not the solution??? With the result X (np.dot(A,X)==b) - is not correct!
[[0.2560342704], [0.7543346092], [-0.084190799732]] [0.27262652 0.61028723 0.74616952] [[0.4201021 ] [0.5206979 ] [0.39195592]] 推荐答案m.axb()当前在Gekko v0.2.8中存在一个错误.它已在下一版Gekko中修复，或者您可以从GitHub获取预发行版本: github/BYU-PRISM/GEKKO/blob/master/gekko/gekko.py ，bin文件夹中有APM可执行文件.该修复程序可用于Gekko v1.0.0.
The m.axb() currently has a bug in Gekko v0.2.8. It is fixed in the next version of Gekko or you can get the pre-release version from GitHub: github/BYU-PRISM/GEKKO/blob/master/gekko/gekko.py with APM executable in the bin folder. The fix is available with Gekko v1.0.0.