/* --------------------------------- */
/* save as fa.h */
/* --------------------------------- */
#define LOOP 2*100
/* --------------------------------- */
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 v(
double y)
{
return( sqrt(9-y*y) );
}
char veq[] = "+sqrt(9-y**2)";
/* --------------------------------- */
double u(
double y)
{
return( -sqrt(9-y*y) );
}
char ueq[] = "-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";
/* ---------------------------------- */
/* ---------------------------------- */
double by = +3.; char byeq[] = "+3";
double ay = -3.; char ayeq[] = "-3";
/* ---------------------------------- */
/* ---------------------------------- */
double b = 2*PI; char beq[] = "2*PI";
double a = 0.; char aeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */