/* ---------------------------------- */
/* save as fa.h                      */
/* --------------------------------- */
/* ---------------------------------- */
double M(
double x,
double y)
{
        return( (x*x+y));
}
char  Meq[] =  "(x**2+y)";
/* ---------------------------------- */
double N(
double x,
double y)
{
        return(  (x*y*y));
}
char  Neq[] =  "(x*y**2)";
/* ---------------------------------- */
/* ---------------------------------- */
double N_x_mns_M_y(
double x,
double y)
{
 return(  (y*y) - 1);
}
/* --------------------------------- */
char  N_x_mns_M_y_eq[] = "(y**2)-1";
/* ---------------------------------- */
/* ---------------------------------- */
double u(
double x)
{
        return( (x*x));
}
char  ueq[] =   "x**2";
/* ---------------------------------- */
double v(
double x)
{
        return( (-x));
}
char  veq[] =   "-x";
/* ---------------------------------- */
/* ---------------------------------- */