code: def major_axis_endpoints(self):
          
  
          i = np.argmax(self.axlens)  # find the major axis
          v = self.paxes[:, i]  # vector from center to major axis endpoint
  
          return self.ctr - v, self.ctr + v
docs: Return the endpoints of the major axis.