/* --------------------------------- */
/* save as fa.h */
/* --------------------------------- */
double f(
double x,
double y)
{
return( (1/(x*x + y*y + 1)) );
}
/* --------------------------------- */
char feq[] = "1/(x*x + y*y + 1)";
/* --------------------------------- */
/* --------------------------------- */