> plot(x^4+x^3-3*x^2+x,x=0..1);

> int(x^4+x^3-3*x^2+x,x=0..1);

> solve(x^4+x^3-3*x^2+x=0);

> interval1:=int(x^4+x^3-3*x^2+x,x=0..sqrt(2)-1);

> interval2:=int(x^4+x^3-3*x^2+x,x=sqrt(2)-1..1);

> interval1+interval2;

> interval1-interval2;

> simplify(interval1-interval2);

> evalf(interval1-interval2);

>