eq1=Eq(0, -sympy.sqrt(x**2) + sympy.sqrt((x + 4)**2))
sympy.solveset(eq1,x)
это должно работать
OUTPUT
{-2}
Кстати, у симпози есть отличный сайт, который бы дал вам ответ, вместе с очень интуитивно понятной оболочкой внутри браузера https://docs.sympy.org/latest/tutorial/solvers.html
вот что говорит симпатия по этому вопросу:
«Обратите внимание, что есть еще одна функция, называемая решением, которая также может использоваться для решения уравнений. Синтаксис решения (уравнения, переменные) Однако рекомендуется использовать вместо нее решение набора.»