/* ---------------------------------- */
/* save as c17b.c                     */
/* ---------------------------------- */
#include "x_hfile.h"
#include       "fb.h"
/* ---------------------------------- */
int main(void)
{
double  a = -2;
double  b =  2; 

double nx = 2*20;
double ny = 2*20;  

double Ix =  0; 
double Iy =  0; 
double I0 =  0; 
  
/* ---------------------------------- */
 clrscrn();

 printf(" Find the polar moment of inertia of the lamina\n");
 printf(" that has the shape of the region bounded by the   \n");
 printf(" graphs of the given equations.\n\n");

 printf(" x = y**2, x = 4,  \n\n");
 printf(" and f(x,y)  (the area mass density at (x,y)) \n\n\n");
 printf(" f : (x,y)-> %s    \n\n", feq);
 printf(" These give :      \n\n");
 printf(" v :   (y)-> %s      \n", veq);
 printf(" u :   (y)-> %s    \n\n", ueq);
 printf(" b :   %+.1f   \n", b);
 printf(" a :   %+.1f \n\n",   a);
 stop();

/* ---------------------------------- */
 clrscrn();
 
   Ix = simpson_dxdy(h,u,v,nx,a,b,ny);
 printf(" Compute the moment of inertia Ix.\n");
 printf("         (b     (v(y)\n");
 printf(" Ix = int(   int(    y**2  f(x,y)  dx dy = %.1f\n", Ix);
 printf("         (a     (u(y)\n\n");

   Iy = simpson_dxdy(g,u,v,nx,a,b,ny);
 printf(" Compute the moment of inertia Iy.\n");
 printf("         (b     (v(y)\n");
 printf(" Iy = int(   int(    x**2  f(x,y)  dx dy = %.1f\n", Iy);
 printf("         (a     (u(y)\n\n");

 printf("   The polar moment of inertia : ");
 printf("     I0 = Ix + Iy = %.1f\n\n",Ix+Iy);

   I0 = simpson_dxdy(gplsh,u,v,nx,a,b,ny);
 printf(" or directly.\n\n");
 printf("         (b     (v(y)\n");
 printf(" I0 = int(   int(     (x**2 + y**2)  f(x,y)  dx dy = %.1f\n", I0);
 printf("         (a     (u(y)\n\n");
 stop();

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