function g=grad(x,y)
%function g=grad(x,y)
g=[-20*x+4*y*cos(x*y)-2+4*x^3
       20*y+4*x*cos(x*y)];
