from sympy import *
m,x = symbols('m x') y = 2**(x-m)+2 expr = Eq(log(8*(x-2),2).subs(x, y),x) expr
nsolve(expr.subs(x,-3),m,0)