尝试使用 GEKKO 求解非线性方程组。系统有多种解决方案，但 GEKKO 给出了错误的解决方案。我该如何解决这个问题？

Question

我正在尝试求解与此类似的多个非线性方程组：

从图中可以看出，系统有两种解决方案。通过多次反复试验和不同的方法（fsolve 等）来实现这一点，我偶然发现了 GEKKO，因此如下：

from gekko import GEKKO
m = GEKKO()
x,y = [m.Var(1) for i in range(2)]
m.Equations([0.2618 == 0.2094 + x*(1-(1/y))**(1/y),\
             0.9200 == 1-math.e**((-((.5236-.2094)/x)**y))])
m.solve(disp=False)
print(x.value,y.value)

此代码成功解决了系统问题，但输出的解决方案不正确。我不太清楚如何获得另一个（这是我需要的）。有 GEKKO 专家可以提供帮助吗？谢谢！

我尝试了多种不同的方法来求解非线性系统（fsolve、Sympy nsolve 等），但始终遇到我无法解决的被零除错误。解决系统的其他方法没有给出解决方案。我找不到 GEKKO 的正确文档来弄清楚如何获得其他解决方案。