integrate_volume.txt   given  the  function  f(x,y,z)
and max and min values,  do numerical integration

  xmax  ymax  zmax
   (     (     (
q= |     |     |  f(x,y,z) dz   dy   dx
   )     )     )
  xmin  ymin  zmin

 |    |     |   z integrate  |    |      |
 |    |	    +----------------+    |      |
 |    |                           |      |
 |    |      y integrate          |      |
 |    +---------------------------+      |
 |                                       |
 |    x integrate                        |
 +---------------------------------------+

 z integrate =   (f(x,y,zmax)-f(x,y,zmin))

 y integrate =  ((f(x,ymax,zmax)-f(x,ymax,zmin))+
                 (f(x,ymin,zmax)-f(x,ymin,zmin)))

 x integrate = (((f(xmax,ymax,zmax)-f(xmax,ymax,zmin))+
                 (f(xmax,ymin,zmax)-f(xmax,ymin,zmin)))-
                ((f(xmin,ymax,zmax)-f(xmin,ymax,zmin))+
                 (f(xmin,ymin,zmax)-f(xmin,ymin,zmin)))) = q

Do numerical check for

 intergral xmin..xmax, ymin..ymax zmin..zmax f(x,y,z) dx dy dz
 f(x,y,z) = x*z+y*y+2*x
 xmin = 0.1;
 xmax = 1.0;
 ymin = 0.0;
 ymax = 1.1;
 zmin = 0.2;
 zmax = 0.9;

  double intfxyz() //  (((f(xmax,ymax,zmax)-f(xmax,ymax,zmin))+
                   //    (f(xmax,ymin,zmax)-f(xmax,ymin,zmin)))-
                   //   ((f(xmin,ymax,zmax)-f(xmin,ymax,zmin))+
                   //    (f(xmin,ymin,zmax)-f(xmin,ymin,zmin))))
    double q;
    q = (((f(xmax,ymax,zmax)-f(xmax,ymax,zmin))+
          (f(xmax,ymin,zmax)-f(xmax,ymin,zmin)))-
         ((f(xmin,ymax,zmax)-f(xmin,ymax,zmin))+
          (f(xmin,ymin,zmax)-f(xmin,ymin,zmin))));
    return q;
  } // intfxyz