/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
#define LOOP 2*50
/* ---------------------------------- */
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";
/* ---------------------------------- */
/* ---------------------------------- */
double b = 1.; char beq[] = "1";
double a = 0.; char aeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */