// pdedef232.h created by  form_eqn23  from user web input
#define nx 5
#define ny 6
#define nz 7
void coefs(double coef[]); 
double ccxx(double x, double y, double z); 
double ccyy(double x, double y, double z);
double cczz(double x, double y, double z);
double ccxy(double x, double y, double z);
double ccxz(double x, double y, double z);
double ccyz(double x, double y, double z);
double ccx(double x, double y, double z);
double ccy(double x, double y, double z);
double ccz(double x, double y, double z);
double ccc(double x, double y, double z);
double ub(double x, double y, double z); // ub(x, y, z) for boundary 
double f(double x, double y, double z); // RHS f(x, y, z)