/* ---------------------------------- */
/* save as c18c.c                   */
/* ---------------------------------- */
#include "x_hfile.h"
#include       "fc.h"
/* ---------------------------------- */
int main(void)
{
double az =  -1;
double bz =   1.;
int nz    = 2*100;
int ny    = 2*100;
double m= 0;

/* ---------------------------------- */
 clrscrn();
 printf(" Let S be the part of the graph of x = %s.  \n\n", keq);
 printf(" If F(x,y,z) = %si %sj %sk, find the flux of F through S\n\n\n",
          Meq,Neq,Peq);
 printf(" Consider g(x,y,z) = x - (%s)\n\n",keq);
 printf("          n        = grad(g(x,y,z)) / ||grad(g(x,y,z))||\n\n\n");

 m = flux_dydz( M,N,P,
                k, u,v,ny,
                  az,bz,nz,
                  H);

 printf(" The flux of F through S is \n\n");
 printf("     //                  \n");
 printf("    ||                   \n");
 printf("    || F.n  dS = %.3f\n",m);
 printf("    ||                   \n");
 printf("   //                    \n");
 printf("   S                 \n\n\n");
 stop();

/* ---------------------------------- */
  clrscrn();
 printf("    Formula : \n\n");
 printf("     //                  \n");
 printf("    ||                   \n");
 printf("    || F.n  dS = %.3f\n",m);
 printf("    ||                   \n");
 printf("   //                    \n");
 printf("   S                 \n\n\n");

 printf("     //\n");
 printf("    ||              (n)\n");
 printf("    || \n");
 printf("    ||     F  .  (+i-k_yj-k_zk)    dS = %.3f\n",m);
 printf("    ||            ------------    \n");
 printf("    ||          [1+k_y^2+k_z^2]^1/2\n");
 printf("   //  \n");
 printf("   S \n\n\n");
 stop();

/* ---------------------------------- */
 clrscrn();
 printf("    Formula :\n\n");
 printf("     //\n");
 printf("    || \n");
 printf("    ||     F  .  (+i-k_yj-k_zk)    dS = %.3f\n",m);
 printf("    ||            ------------   \n");
 printf("    ||           [1+k_y^2+k_z^2]^1/2\n");
 printf("    || \n");
 printf("   //  \n");
 printf("   S \n\n\n");

 printf("     //\n");
 printf("    ||    (F)       (n)                  (dS)\n");
 printf("    || \n");
 printf("    ||     F  .  (+i-k_yj-k_zk)   [1+k_y^2+k_z^2]^1/2 dA = %.3f\n",m);
 printf("    ||            -----------     \n");
 printf("    ||           [1+k_y^2+k_z^2]^1/2\n");
 printf("   //  \n    ");
 printf("   Rxy \n\n\n");
 stop();

/* ---------------------------------- */
 clrscrn();
 printf("     / b   / v(z)\n");
 printf("    |     |      \n");
 printf("    |     |     F.(+i-k_yj-k_zk)    [1+k_y^2+k_z^2]^1/2 dy dz = %.3f\n",m);
 printf("    |     |        -----------     \n");
 printf("    |     |      [1+k_y^2+k_z^2]^1/2\n");
 printf("    |     |      \n");
 printf("   /  a  /   u(z)\n\n\n");

 printf(" With.\n\n\n");
 printf(" F : (x,y,z)-> %si %sj %sk \n\n",Meq,Neq,Peq);
 printf(" k :   (y,z)-> %s    \n\n", keq);
 printf(" u :     (z)-> %s      \n", ueq);
 printf(" v :     (z)-> %s    \n\n", veq);
 printf(" a = %+.1f   b = %+.1f \n",az,bz);
 stop();

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