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

/* ---------------------------------- */
/* ---------------------------------- */
double M(
double x,
double y)
{
        return( (2*x+y*y*y));
}
char  Meq[] =  "(2*x + y**3)";
/* ---------------------------------- */
double N(
double x,
double y)
{
        return(  (3*x*y*y  + 4));
}
char  Neq[] =  "(3*x*y**2 + 4)";
/* ---------------------------------- */
/* ---------------------------------- */
double g(
double t)
{
        return( (t));
}
char  geq[] =  "t";
/* ---------------------------------- */
double h(
double t)
{
        return( (t+1));
}
char  heq[] =  "t+1";
/* ---------------------------------- */
/* ---------------------------------- */