function output = estimateInsertion(obj)
  humerus = obj.container.shoulder.humerus;

  scapulaCoordSys = humerus.shoulder.scapula.coordSys;
  humerusSuperior = humerus.center + humerus.radius * scapulaCoordSys.IS;
  humerusAnterior = humerus.center + humerus.radius * scapulaCoordSys.PA;
  humerusPosterior = humerus.center - humerus.radius * scapulaCoordSys.PA;
  humerusLateral = humerus.center + humerus.radius * scapulaCoordSys.ML;

  switch obj.getFullName()
    case "subscapularis"
      muscleInsertion = humerusAnterior;
    case "supraspinatus"
      muscleInsertion = humerusSuperior;
    case "infraspinatus"
      muscleInsertion = humerusPosterior;
    case "teres_minor"
      muscleInsertion = humerusPosterior;
    otherwise
      muscleInsertion = humerusLateral;
  end

  output = muscleInsertion;
end