/* ---------------------------------- */
/* 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 ds(
double x)
{
return( sqrt(1. + pow(fx_x(f,x,H),2)) );
}
/* ---------------------------------- */
/* ---------------------------------- */