/* ---------------------------------- */
/* save as c17c.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "fc.h"
/* --------------------------------- */
int main(void)
{
double az,bz,nz,
       nx,ny,M;
 
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);
 
 printf(" v : x,z   -> %s  \n", veq);
 printf(" u : x,z   -> %s\n\n", ueq);

 printf(" t : z     -> %s  \n", teq);
 printf(" s : z     -> %s\n\n", seq);
 
 az = 0;
 bz = 1;
 nz = 2*50;nx = 2*50;ny = 2*50;

 M =  trapezoid_dydxdz(     f,
                            az, bz, nz,
                            sx, tx, nx,
                            uy, vy, ny);

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

 stop();

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