code: def point_displ(pt1, pt2):
      
  
      #Imports
      import numpy as np
  
      # Make iterable
      if not np.iterable(pt1):
          pt1 = np.float64(np.array([pt1]))
      else:
          pt1 = np.float64(np.array(pt1).squeeze())
      ## end if
      if not np.iterable(pt2):
          pt2 = np.float64(np.array([pt2]))
      else:
          pt2 = np.float64(np.array(pt2).squeeze())
      ## end if
  
      # Calculate the displacement vector and return
      displ = np.matrix(np.subtract(pt2, pt1)).reshape(3,1)
      return displ
docs: Calculate the displacement vector between two n-D points.
  
      pt1 - pt2
  
      .. todo:: Complete point_disp docstring