/* --------------------------------- */
/* save as fb.h */
/* --------------------------------- */
#define LOOP 2*300
/* ---------------------------------- */
double f(
double x,
double y)
{
return( sqrt(pow(x*x+y*y,3) + 1) );
}
char feq[] = "((x**2 + y**2)**3) + 1)**(1/2)";
/* --------------------------------- */
/* --------------------------------- */
double v(
double x)
{
return( (sqrt(1-x*x)) );
}
char veq[] = "+(1-x**2)**(1/2)";
/* --------------------------------- */
double u(
double x)
{
return( (-sqrt(1-x*x)) );
}
char ueq[] = "-(1-x**2)**(1/2)";
/* ---------------------------------- */
/* ---------------------------------- */
double bx = +1.; char bxeq[] = "+1";
double ax = -1.; char axeq[] = "-1";
/* ---------------------------------- */
/* ---------------------------------- */
/* --------------------------------- */
/* --------------------------------- */
/* Polar form : *(r) dkdr */
double h(
double r,
double k)
{
return( sqrt(pow(r,6) + 1) );
}
char heq[] = "(r**6 + 1)**(1/2)";
/* --------------------------------- */
/* ----------- dkdr ---------------- */
double t(
double r)
{
return( (2*PI) );
}
char teq[] = "2*Pi";
/* ---------------------------------- */
double s(
double r)
{
return( (0) );
}
char seq[] = "0";
/* ---------------------------------- */
/* ------------ dkdr ---------------- */
double br = 1.; char breq[] = "1";
double ar = 0.; char areq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */