SymPy и mpmath дают "TypeError: не может создать MPF" при использовании функции егГ() в пределах solveset()

Question

У меня есть 4 входных переменных (поплавки):

• Xmax
• Xmin
• процент
• режим

и я хочу, чтобы решить следующую (довольно долго) уравнение с:

> (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage == 0 

Я хочу использовать mpmath и SymPy решить уравнение, но это дает мне следующее сообщение об ошибке:

TypeError: cannot create mpf from 0.707106781186547*(-s**2 - 0.287682072451781)/s

Мой код выглядит следующим образом:

from mpmath import erf, log, sqrt 
from sympy import Symbol, solveset, S 

percentage = 0.95 
mode = 2 
Xmin = 1. 
Xmax = 1.5 
s = Symbol('s') 

eqn = (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage) 

solveset(eqn, s, domain=S.Reals) 

MPF является поплавок тип, созданный mpmath.

Я думаю, что я сузил проблему функции ERF(), она возвращает

EmptySet() 

когда я бегу

solveset(log(Xmax) - (log(mode) + s ** 2), s, domain=S.Reals) 

Я не могу понять, что попробовать дальше, любая помощь будет оценили!

Первоначально я думал, что это проблема с математикой, но уравнение успешно решено в Matlab, поэтому проблема, вероятно, исходит из sympy или mpmath.

Answer 1

Потеря mpmath импорта erf и использование версии sympy устраняет вашу ошибку.

from sympy import Symbol, solveset, S, erf, log, sqrt 

percentage = 0.95 
mode = 2 
Xmin = 1. 
Xmax = 1.5 
s = Symbol('s', real=True) 

eqn = (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage) 

solveset(eqn, s) 

Обратите также внимание:

• вам не придется импортировать log и sqrt из mpmath. Это не повлияет на ваш результат здесь, чтобы получить их от .
• вы можете указать реальный домен на переменной s, это позволит вам сделать это на вызовах solveset.

Пример дальнейшего использования в тестах упаковки here, если вам это нужно.