辛普森不能解一组方程

Question

我有一个模式，一个最大值和最小值X（Xmin和Xmax），以及一个百分数（percentage）。在

我想使用以下函数来找到理论对数正态分布的μ和{}：

累积分布函数：

模式呢

我决定把它作为一组方程来求解。在这个例子中，我使用mode、percentage、Xmin和{}的任意值

from sympy import *

def CalcScaleParamSET (mode, percentage, Xmin, Xmax):

    s = Symbol('s', Real=True)
    mu = Symbol('mu', Real =True)

    eqn1 = (1/2+1/2*erf((log(Xmax)-mu)/sqrt(2)*s)-1/2-1/2*erf((log(Xmin)-mu)/sqrt(2)*s)) -percentage
    eqn2 = (exp(mu-s**2)) - mode

    eqns = [eqn1, eqn2]

    ans = linsolve(eqns, s, mu)

    print(ans)

CalcScaleParamSET(2, 0.95, 1, 4)

这将返回一个EmptySet()，这意味着sympy无法找到方程的解。我不知道我做错了什么，是数学问题还是代码问题？我编写了一个matlab脚本来执行同样的任务，它能够得到一个单独的数值μ和{}。在