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