function y=fdes(x)
y = 1- cos(x(1)) * cos(x(2)) ;