clear
//Given 
//Entire area - hallow area
l_e = 60.0     //mm - length of the entire area
b_e = 40       //mm - width of the entire area
l_h = 30       //mm - length of the hallow area
b_h = 20       //mm - width of the hallow area
A_e = l_e*b_e  //mm2 - The entire area
A_h = -l_h*b_h //mm2 - The hallow area '-' because its hallow
A_re = A_e + A_h //mm2 resultant area
y_e = l_e/2      //  mm com from bottom 
y_h = 20+l_h/2   //mm com from bottom 
y_com = (A_e*y_e + A_h*y_h)/A_re 
//moment of inertia caliculatins - bh3/12 +ad2
I_e = b_e*(l_e**3)/12 + A_e*((y_e-y_com)**2)  //Parallel axis theorm
I_h =  b_h*(l_h**3)/12 - A_h*((y_h-y_com)**2) //Parallel axis theorm
I_total = I_e - I_h
printf("\n The moment of inertia of total system is  %e mm^4",I_total)