//moment of inertia of shaded region
//refer fig. 9.44
//The figure is divided into simple geometry
IAB=((80*80^3)/12)+((%pi*80^4)/128)-((%pi*40^4)/64)  //mm^4
printf("IAB=%.2d mm^4",IAB)