/* ---------------------------------- */
/* save as fa.h                      */
/* --------------------------------- */
double M(
double x,
double y,
double z)
{
      return(   (3*z) );
}
char  Meq[] = "+ 3*z";
/* --------------------------------- */
double N(
double x,
double y,
double z)
{
      return(   (4*x) );
}
char  Neq[] = "+ 4*x";
/* --------------------------------- */
double P(
double x,
double y,
double z)
{
      return(    (2*y) );
}
char  Peq[] = " + 2*y";
/* --------------------------------- */
/* --------------------------------- */
double f(
double x,
double y)
{
      return( (9-x*x-y*y) );
}
char  feq[] = "9-x**2-y**2";
/* --------------------------------- */
/* --------------------------------- */
double u(
double y)
{
 return(       -sqrt(9-y*y) );
}
char  ueq[] = "-sqrt(9-y**2)";
/* --------------------------------- */
double v(
double y)
{
 return(       sqrt(9-y*y)  );
}
char  veq[] = "+sqrt(9-y**2)";
/* --------------------------------- */
/* --------------------------------- */
double g(
double t)
{
        return( 3*cos(t));
}
char  geq[] =  " 3*cos(t)";
/* --------------------------------- */
double h(
double t)
{
        return( 3*sin(t));
}
char  heq[] =  " 3*sin(t)";
/* --------------------------------- */
double k(
double t)
{
        return( 0);
}
char  keq[] =  "0";
/* --------------------------------- */
/* --------------------------------- */