function z = func(x,y)

z = 1./((x-y).^2 +1).* exp(-4*(x+y).^2);
