code: def squarelim():
      
      fig = gcf()
      xmin, xmax = fig.xlim
      ymin, ymax = fig.ylim
      zmin, zmax = fig.zlim
      width = max([abs(xmax - xmin), abs(ymax - ymin), abs(zmax - zmin)])
      xc = (xmin + xmax) / 2
      yc = (ymin + ymax) / 2
      zc = (zmin + zmax) / 2
      xlim(xc - width / 2, xc + width / 2)
      ylim(yc - width / 2, yc + width / 2)
      zlim(zc - width / 2, zc + width / 2)
docs: Set all axes with equal aspect ratio, such that the space is 'square'.