syms x
expr = x^2;
F = int(expr,[0 1])