/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
double a_n(
double n)
{
return( pow(-1.,n)*(n*n+4)/pow(2.,n) );
}
char a_neq[] = "(-1)**n * (n**2 +4)/ 2**n";
/* ---------------------------------- */
double a_npls1(
double n)
{
return( pow(-1.,n)*((n+1)*(n+1)+4)/pow(2.,(n+1)) );
}
char a_npls1eq[] = "(-1)**n * ((n+1)**2+4)/ 2**(n+1)";
/* ---------------------------------- */
/* ---------------------------------- */