/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
double f(
double x)
{
return(x + 1.);
}
char feq[] = "x + 1.";
/* ---------------------------------- */
double Df(
double x)
{
return(1.);
}
char Dfeq[] = "1.";
/* ---------------------------------- */
double ArcLength(
double x)
{
return( sqrt(1. + pow(Df(x),2)) );
}
/* ---------------------------------- */
/* ---------------------------------- */