/* --------------------------------- */
/* save as fa.h */
/* --------------------------------- */
#define LOOP 2*50
/* ---------------------------------- */
double f(
double x,
double y)
{
return( (x*x+y*y+1) );
}
char feq[] = "x**2+y**2+1";
/* ---------------------------------- */
/* ---------------------------------- */
double by(
double x)
{
return( (2-2*x) );
}
char byeq[] = "2-2*x";
/* ---------------------------------- */
double ay(
double x)
{
return( ( 0) );
}
char ayeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */
double bx = 1.; char bxeq[] = "1";
double ax = 0.; char axeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */