/* ---------------------------------- */
/* save as c17f.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "ff.h"
/* --------------------------------- */
int main(void)
{
double ax,bx,nx, ny, nz, M;
 
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);

 printf(" v : x,y   -> %s  \n", veq); 
 printf(" u : x,y   -> %s\n\n", ueq);

 printf(" t : x     -> %s  \n", teq); 
 printf(" s : x     -> %s\n\n", seq);

 ax = 0; 
 bx = 1; 
 nx = 2*50; ny = 2*50; nz = 2*50;

 M =    trapezoid_dzdydx(   f,
                            ax, bx, nx,
                            sy, ty, ny,
                            uz, vz, nz);

 printf(" With the trapezoidal's rule.\n\n\n");
 printf("        (%+.1f  (%s     (%s    \n", bx, teq, veq);
 printf(" A = int(   int(    int(      %s  dz dy dx = %.6f\n",feq, M);
 printf("        (%+.1f  (%s     (%s\n\n\n", ax, seq, ueq);

 stop();

 return 0;
}
/* --------------------------------- */