def initialize_axes(self, boundary=<NUM_LIT>):
bs = [(self.vals[:, i].max()-self.vals[:, i].min())*boundary for i in range(<NUM_LIT:3>)]<EOL>self.x_lim = np.array([self.vals[:, <NUM_LIT:0>].min()-bs[<NUM_LIT:0>], self.vals[:, <NUM_LIT:0>].max()+bs[<NUM_LIT:0>]])<EOL>self.y_lim = np.array([self.vals[:, <NUM_LIT:1>].min()-bs[<NUM_LIT:1>], self.vals[:, <NUM_LIT:1>].max()+bs[<NUM_LIT:1>]])<EOL>self.z_lim = np.array([self.vals[:, <NUM_LIT:2>].min()-bs[<NUM_LIT:2>], self.vals[:, <NUM_LIT:2>].max()+bs[<NUM_LIT:2>]])<EOL>
Set up the axes with the right limits and scaling.