syms x y
f = x^2 + y^2;
F = int(f,x, [0 1])
G = int(F, y, [0, 1])