/* abc2.h  just the functions that need to be computed, not the solver  */
#define nd 5
#define xmin -1.0
#defin2 xmax  1.0
#define ymin -1.0
#define ymax  1.0
#define nx      7
#define ny      7
#define ifcheck 5

double a1(double x, double y);
double b1(double x, double y);
double c1(double x, double y);
double d1(double x, double y);
double e1(double x, double y);
double f1(double x, double y);
double u(double x, double y);
double c(double x, double y);