SUBROUTINE FUKUI(FK,KIND)
C
C     PROCEDURA WYZNACZA WEKTOR JADROWYCH FUNKCJI FUKUI'EGO.
C     TYMCZASOWO OGRANICZONE DO RHF I UHF.
C     JEST TO MODYFIKACJA PROCEDURY HELFEY Z PLIKU GRD1.SRC
C
C     !!!!!!!!!!!!!!!!!!!!!!!!!!!! UWAGA !!!!!!!!!!!!!!!!!!!!!!!!!!!!
C     KONIECZNIE PRZED UZYCIEM TEJ PROCEDURY WYWOLAJ PROCEDURE HRDNSS
C     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
C
C     KIND=                         ATAK:
C          1: FUKUI= F(0)-F(+)          ELEKTROFILOWY
C          2: FUKUI= F(-)-F(0)          NUKLEOFILOWY
C          3: FUKUI=(F(-)-F(+))/2       RODNIKOWY
C     GDZIE: F - SILY DZIAJACAJE NA JADRA ATOMOW (NAJCZESCIEJ: F(0)=0)
C            0 - WYJSCIOWY UKLAD
C            + - UKLAD POZBAWIONY ELEKTRONU
C            - - UKLAD Z DODATKOWYM ELEKTRONEM
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
      INTEGER MOINFO(8),FAHOMO,CAHOMO,FBHOMO,CBHOMO,
     *     FALUMO,CALUMO,FBLUMO,CBLUMO
      EQUIVALENCE (MOINFO(1),FAHOMO),(MOINFO(2),CAHOMO),
     *     (MOINFO(3),FBHOMO),(MOINFO(4),CBHOMO),
     *     (MOINFO(5),FALUMO),(MOINFO(6),CALUMO),
     *     (MOINFO(7),FBLUMO),(MOINFO(8),CBLUMO)
      COMMON /HOLUMO/ MOINFO
C
      LOGICAL IANDJ,OUT,DBG,NORM,DOUBLE,GOPARR,DSKWRK,MASWRK,NXT,SVDSKW
C
C  MAGIC NUMBERS: S(1)+P(3)+D(6)+F(10)+G(15)=35
C                 BASIS FUNCTION TYPES=5
C                 MAXIMUM BASIS TYPE SQUARED = (15)**2 = 225   <-- czy tu chodzi o kwadrat maksymalnej liczby funkcji na podpowloce??

=========================================================================


C     --- WLASNY SPOSOB GENEROWANIA NIBY MACIERZY GESTOSCI ---
C
C     ZAMIAST:
*                NDUM = IA(LOCI+I)+(LOCJ+J)  <----  cos zgubione po + 
*                DEN = DENAB(NDUM)
C     JEST:
              DEN=0.0D0
              DO INMOS=1,NMOS
                 DEN=DEN+X(LC+INMOS*NUM+LOCI+I)*
     *                X(LC+INMOS*NUM+LOCJ+J)
              ENDDO