/* ---------------------------------- */
/* save as fa.h */
/* ---------------------------------- */
#define LOOP 2*50
/* ---------------------------------- */
double M(
double x,
double y)
{
return( (y-sin(x)));
}
char Meq[] = "(y-sin(x))";
/* ---------------------------------- */
double N(
double x,
double y)
{
return( (cos(x)));
}
char Neq[] = "(cos(x))";
/* ---------------------------------- */
/* ---------------------------------- */
double N_x_mns_M_y(
double x,
double y)
{
return( (-sin(x))-(1));
}
char N_x_mns_M_y_eq[] = "(-sin(x))-(1)";
/* ---------------------------------- */
/* ---------------------------------- */
double v(
double x)
{
return( (2*x/PI));
}
char veq[] = "2*x/Pi";
/* ---------------------------------- */
double u(
double x)
{
return( (0));
}
char ueq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */
double bx = PI/2.; char bxeq[] = "PI/2";
double ax = 0.; char axeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */