function z = fun1(x,y)
% function z = fun1(x,y) computes sin(2*pi*x)*cos(4*pi*y)
z = sin(2*pi*x).*cos(4*pi*y);
