function out = integrand(x, y)
out = y*sin(x) + x*cos(y);