function measureHumerusContactPoint(obj)
  humerus = obj.container.shoulder.humerus;
  humeralHeadContactPointFinder = SphereContactPoint(...
    Sphere(humerus.center, humerus.radius));

  % muscles humeral head insertion is estimated to be the most lateral point of
  % the humeral head fitted sphere 
  MLAxis = obj.container.shoulder.scapula.coordSys.ML;
  muscleInsertion = humerus.center + humerus.radius * MLAxis;
  humeralHeadContactPointFinder.setAnchorPoint(muscleInsertion);

  humeralHeadContactPointFinder.setPolarPoint(obj.centroid);

  obj.forceApplicationPoint = humeralHeadContactPointFinder.getContactPoint();
end