docs: def molecular_orbital(coords, mocoeffs, gbasis):
      '''Return a molecular orbital given the nuclei coordinates, as well as
         molecular orbital coefficients and basis set specification as given by the cclib library.
  
         The molecular orbital is represented as a function that takes x, y, z coordinates (in a vectorized fashion)
         and returns a real number.
  
      '''
      
      # Making a closure
      def f(x, y, z, coords=coords, mocoeffs=mocoeffs, gbasis=gbasis):
          # The other functions take nanometers
          return sum(c * bf(x * 10, y*10, z*10) for c, bf in zip(mocoeffs, getbfs(coords * 10, gbasis)))
  
      return f