/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
double f(
double x,
double y)
{
return( (y*y*y*y + 3*y - 4*x*x*x - 5*x - 1) );
}
char feq[] = "y**4 + 3*y - 4*x**3 - 5*x - 1";
/* ---------------------------------- */
double dy(
double x,
double y)
{
return( (12*x*x + 5) / (4*y*y*y + 3) );
}
char dyeq[] = "(12*x**2 + 5) / (4*y**3 + 3)";
/* ---------------------------------- */
/* ---------------------------------- */