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

#include "pdedef23.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 sin(x);
} // end ccxx

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

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

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

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

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

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

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

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

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

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

double f(double x, double y, double z) // RHS f(x, y, z)
{
  return  (-sin(x))*sin(x + y + z) - 
          sin(y)*sin(x + y + z) -
          cos(z)*sin(x + y + z) -
          sin(x + y + z) + 
          18.0*cos(x + y + z);
} // end f

// end  pdedef23.c