sympy really wants to simplify by pulling terms from sqrt , which makes sense. I think you should do what you want manually, i.e. Get the simplification you want without calling sqrt , and then release it using Symbol with the LaTex \sqrt wrapper. For instance:
from sympy import * init_printing(use_latex='mathjax')
This is really not perfect, but I looked through the documents for about an hour and simply could not find support for what you want directly. The sympy library is more about actual symbolic manipulations, less about printing, so I can hardly blame them.
source share