docs: def potbeslsho(x, y, z1, z2, labda, order, ilap, naq):
      """
      Parameters
      ---------- 
         x,y: Point where potential is computed
         z1: Complex begin point of line-sink
         z2: Complex end point of line-sink
         labda(naq): labda's (zero for first labda if Laplace)
         order: Order of the line-sink
         ilap: equals 1 when first value is Laplace line-sink and first labda equals zero
         naq: Number of aquifers
         rv(naq): Array to store return value (must be pre-allocated)
      Returns
      --------
         rv(naq): Potentials. First spot is Laplace value if ilap=1
      """
  
      rv = np.zeros(naq)
  
      # lstype = 1 means line-sink
      lstype = 1
  
      # Radius of convergence
      Rconv = 7.0
  
      # if (ilap==1) :
      #    istart = 1
      # else:
      #    istart = 0
      zin, z1in, z2in, Lin, z, zplus1, zmin1 = prepare_z(x, y, z1, z2)
      # Laplace linesink
      if ilap == 1:
          power = order + 1
          pcor = np.complex(0.0, 0.0)
          for n in range(1, int((power + 1) / 2) + 1):
              pcor = pcor + z ** (power - 2 * n + 1) / (2 * n - 1)
          pcor = 2.0 * pcor
          comega = (
              z ** power * np.log((zmin1) / (zplus1))
              + pcor
              - np.log(zmin1)
              + (-1.0) ** power * np.log(zplus1)
          )
          comega = -comega * Lin / (4.0 * np.pi * power)
          rv[0] = np.real(comega)
      # N-1 leakage factors
      for i in range(ilap, naq):
          # Check whether entire linesink is outside radius of convergence
          # Outside if |z-zc|>L/2+7lab, and thus |Z|>1+7lab*2/L, or |zeta|>1/biglab+7 (zeta is called z here)
          biglab = 2.0 * labda[i] / Lin
          z = (2.0 * zin - (z1in + z2in)) / (z2in - z1in) / biglab
  
          if abs(z) < (Rconv + 1.0 / biglab):
              pot = IntegralF(zin, z1in, z2in, Lin, labda[i], order, Rconv, lstype)
              rv[i] = -Lin / 2.0 * pot
          else:
              rv[i] = 0.0
      return rv