# include <stdio.h>

int main ( )
{
  float a;
  float b;
  float c;
  float d;
  float dx;
  float dy;
  float fxy;
  float fxymin;
  int i;
  int j;
  int m;
  int n;
  float x;
  float xmin;
  float y;
  float ymin;

  a = -10.0;
  b = +10.0;

  c = -10.0;
  d = +10.0;

  m = 11;
  n = 11;

  dx = ( b - a ) / ( m - 1 );
  dy = ( d - c ) / ( n - 1 );
//
//  Initialize the minimum by computing the value at (a,c).
//
  x = a;
  y = c;

  fxymin = x * x + y * y + 2 * x + 8 * y;
  xmin = x;
  ymin = y;
//
//  Now search for the minimum.
//
  for ( i = 0; i < m; i++ )
  {
    x = a + i * dx;
    for ( j = 0; j < n; j++ )
    {
      y = b + j * dy;
      fxy = x * x + y * y + 2 * x + 8 * y;

      printf ( "F(%f,%f) = %f\n", x, y, fxy );

      if ( fxy < fxymin )
      {
        fxymin = fxy;
        xmin = x;
        ymin = y;
      }
    }
  }

  printf ( "F(%f,%f) = %f\n", xmin, ymin, fxymin );

  return 0;
}