expand((e^x+x)^5);

sol: ode2('diff(y, x) = x+y, y, x);

ic1(sol, x=0, y=1);

jacobian([r*cos(theta), r*sin(theta)], [r, theta]);

determinant(%);

trigsimp(%);