/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
double f(
double x)
{
return(x*x + 2.);
}
char feq[] = "x**2 + 2.";
/* ---------------------------------- */
double VCircularDisk(
double x)
{
return(PI * pow(f(x),2) ); // * dx
}
/* ---------------------------------- */
/* ---------------------------------- */