x=var('x')
W=det(matrix([[exp(x),x*exp(x)],[exp(x),(x+1)*exp(x)]]))
print W
W.simplify_rational()