// pdedef232.c created by  form_eqn23  from user web input

#include "pdedef232.h" // has define nx, ny, nz
#include <math.h>

void coefs(double coef[]) 
{
  double xmin=0.0;
  double xmax=1.0;
  double ymin=0.1;
  double ymax=1.1;
  double zmin=0.2;
  double zmax=1.2;
  coef[0] = xmin;
  coef[1] = xmax;
  coef[2] = nx;
  coef[3] = ymin;
  coef[4] = ymax;
  coef[5] = ny;
  coef[6] = zmin;
  coef[7] = zmax;
  coef[8] = nz;
} // end ccoefs

double ccxx(double x, double y, double z) 
{
  return x;
} // end ccxx

double ccyy(double x, double y, double z) 
{
  return y;
} // end ccyy

double cczz(double x, double y, double z) 
{
  return 2.0*z;
} // end cczz

double ccxy(double x, double y, double z) 
{
  return 3.0*y*z;
} // end ccxy

double ccxz(double x, double y, double z) 
{
  return 4.0*x*z;
} // end ccxz

double ccyz(double x, double y, double z) 
{
  return 5.0*x*z;
} // end ccyz

double ccx(double x, double y, double z) 
{
  return 6.0*x*x;
} // end ccx

double ccy(double x, double y, double z) 
{
  return 7.0*y*y;
} // end ccy

double ccz(double x, double y, double z) 
{
  return 8.0*x*y;
} // end ccz

double ccc(double x, double y, double z) 
{
  return 9.0*x*y*z;
} // end ccc

double ub(double x, double y, double z) // ub(x, y, z) for boundary 
{
  return x*x*y+2.0*y*y*z+3.0*x*z*z;
  
} // end ub

double f(double x, double y, double z) // RHS f(x, y, z)
{
  return 2.0*x*y + 4.0*y*z + 12.0*x*z + 26.0*y*z*x + 24.0*x*z*z +
         6.0*x*x*(2.0*x*y + 3.0*z*z) + 7.0*y*y*(x*x  + 4.0*y*z) +
         8.0*x*y*(2.0*y*y + 6.0*x*z) + 
         9.0*x*y*z*(x*x*y + 2.0*y*y*z + 3.0*x*z*z);
} // end f

// end  pdedef232.c