/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
#define LOOP 2*100
/* ---------------------------------- */
double h(
double r,
double k)
{
return( (4 - r*r) );
}
char heq[] = "4 - r**2";
/* ---------------------------------- */
/* ---------------------------------- */
double t(
double k)
{
return( (2) );
}
char teq[] = "2";
/* ---------------------------------- */
double s(
double k)
{
return( (0) );
}
char seq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */
double bk = PI/2.; char bkeq[] = "PI/2";
double ak = 0.; char akeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */