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