function value = f ( x, y )

  m = 10;

  value = - sin ( x ) .* ( sin (     x.^2 / pi ) ).^( 2 * m ) ...
          - sin ( y ) .* ( sin ( 2 * y.^2 / pi ) ).^( 2 * m );

  return
end
