function z = fgrad(delta,x,y,gradx,grady)

z = -func(x + delta*gradx,y + delta*grady);
