/* ---------------------------------- */
/* save as fa.h                   */
/* --------------------------------- */
double M(
double x,
double y,
double z)
{
      return( (y*sin(x)) );
}
char  Meq[] = "+ y*sin(x)";
/* --------------------------------- */
double N(
double x,
double y,
double z)
{
      return( (y*y*z) );
}
char  Neq[] = "+ y^2*z";
/* --------------------------------- */
double P(
double x,
double y,
double z)
{
      return( (x+3*z) );
}
char  Peq[] = " + (x+3*z)";
/* --------------------------------- */
double u(
double x,
double y)
{
 return(     (-1));
}
char  ueq[] = "-1";
/* --------------------------------- */
double v(
double x,
double y)
{
 return(      (1) );
}
char  veq[] = "+1";
/* --------------------------------- */
double s(
double y)
{
 return(       -1 );
}
char  seq[] = "-1";
/* --------------------------------- */
double t(
double y)
{
 return(       1  );
}
char  teq[] = "+1";
/* --------------------------------- */
/* --------------------------------- */