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

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