function [ fxy_min, t ] = min1 ( n )

  if ( nargin < 1 )
    n = 101;
  end

  tic;

  fxy_min = Inf;

  for j = 1 : n
    for i = 1 : n
      x = pi * ( i - 1 ) / ( n - 1 );
      y = pi * ( j - 1 ) / ( n - 1 );
      fxy = f1 ( x, y );
      fxy_min = min ( fxy_min, fxy );
    end
  end

  t = toc;

  return
end
