function J = calcJ( m, ra, ri ) 
  J = m/2* ( ra^2 + ri^2 );
 end 