/* ---------------------------------- */
/* save as c17a.c                     */
/* ---------------------------------- */
#include "x_hfile.h"
#include       "fa.h"
/* ---------------------------------- */
int main(void)
{
double  a = 0;
double  b = 9;  

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(" y = sqrt(x), x = 9, y = 0,  \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 :   (x)-> %s      \n", veq);
 printf(" u :   (x)-> %s    \n\n", ueq);
 printf(" b :   %+.1f   \n", b);
 printf(" a :   %+.1f \n\n",   a);
 stop();

/* ---------------------------------- */
 clrscrn();

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

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

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


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