/* --------------------------------- */
/* save as fa.h */
/* --------------------------------- */
double f(
double x,
double y)
{
return( (x*x + y*y - 1) );
}
char feq[] = "x**2 + y**2 - 1";
/* --------------------------------- */
double g(
double x,
double y)
{
return( ( y - sin(x) ) );
}
char geq[] = "y - sin(x)";
/* --------------------------------- */