/* ---------------------------------- */
/* save as fa.h                       */
/* ---------------------------------- */
double f(
double x,
double y)
{
 return(      (x + y) );
}
/* ---------------------------------- */
char  feq[] = "x + y";
/* ---------------------------------- */
/* ---------------------------------- */
double u(
double x)
{
 return( ( 0) );
}
/* ---------------------------------- */
char  ueq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */
double v(
double x)
{
 return(      sqrt(x) );
}
/* ---------------------------------- */
char  veq[] = "sqrt(x)";
/* ---------------------------------- */
/* ---------------------------------- */
double g(
double x,
double y)
{
 return(       x*x*   (x + y) );
}
/* ---------------------------------- */
char  geq[] = "x**2   (x + y)";
/* ---------------------------------- */
/* ---------------------------------- */
double h(
double x,
double y)
{
 return(       y*y*   (x + y) );
}
/* ---------------------------------- */
char  heq[] = "y**2   (x + y)";
/* ---------------------------------- */
/* ---------------------------------- */
double gplsh(
double x,
double y)
{
 return(   (x*x + y*y)*   (x + y) );
}
/* ---------------------------------- */
char  gplsheq[] = "(x**2 + y**2) (x + y)";
/* ---------------------------------- */
/* ---------------------------------- */