syms x y f = exp(-x*y); F = int( int(f, x, 0, x), y , 0, y); F = matlabFunction(F); fsurf(F, [-1, 1, -1, 1])