/* ---------------------------------- */
/* save as c18a.c                   */
/* --------------------------------- */
#include "x_hfile.h"
#include       "fa.h"
/* --------------------------------- */
int main(void)
{
double ax =  -1;
double bx =   1.;
int    nx = 2*50;

int ny = 2*50;
int nz = 2*50;

double m = 0;
/* --------------------------------- */ 
 clrscrn();
 printf(" Use the divergence theorem to find,\n\n");

 printf(" the flux of F through S.\n\n");

 printf("     //            ///         \n");
 printf("    ||            |||          \n");
 printf("    || F.n dS  =  ||| div F dV \n");
 printf("    ||            |||          \n");
 printf("   //            ///           \n");
 printf("   S                   Q        \n\n\n");

 printf("      If F = Mi + Nj + Pk                    \n\n\n");

 printf("    ///             ///                      \n");
 printf("   |||             |||                       \n");
 printf("   ||| div F dV =  |||  M_x + N_y + P_z dV   \n");
 printf("   |||             |||                       \n");
 printf("  ///             ///                        \n");
 printf("   Q               Q                     \n\n\n");
 stop();

/* --------------------------------- */
 clrscrn();
 
 m = flux_dzdydx(
 
                         M,N,P,
                         ax,bx,nx,
                         s, t, ny,
                         u,v,nz,
                           H);

 printf("     / b   / t(x)  / v(x, y)                       \n");
 printf("    |     |       |                                \n");
 printf("    |     |       |     M_x N_y P_z dzdydx = %.3f  \n",m);
 printf("    |     |       |                                \n");
 printf("   /  a  /   s(x) /  u(x, y)                   \n\n\n");

 printf(" With.\n\n\n");
 printf(" F : (x,y,z)-> %si %sj %sk \n\n",Meq,Neq,Peq);
 printf(" u :   (x,y)-> %s      \n", ueq);
 printf(" v :   (x,y)-> %s    \n\n", veq);
 printf(" s :     (x)-> %s      \n", seq);
 printf(" t :     (x)-> %s    \n\n", teq);
 printf(" ax = %+.1f   bx = %+.1f\n\n",ax,bx);
 printf(" dV = dz dy dx        \n\n");
 stop();

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