/* ---------------------------------- */
/* save as c17e.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "fe.h"
/* --------------------------------- */
int main(void)
{
double ay,by,n,ny,M;
 
 clrscrn();
 printf(" f : (x,y)-> %s\n\n\n", feq);
 printf(" v :   (y)-> %s    \n", veq);
 printf(" u :   (y)-> %s  \n\n", ueq);


 ay = 0;   by = 1;
 n = 2*5; ny = 2*5;
 
  M = simpson_dxdy(f, u,v,n, ay,by,ny);

 printf(" With the simpson's rule.\n\n");
 printf("    (%.3f  (%s\n", by, veq);
 printf(" int(   int(       %s  dx dy = %.6f\n", feq, M);
 printf("    (%.3f  (%s\n\n\n", ay, ueq);

 stop();

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