吕西林的《建筑结构抗震设计理论与实例》中多自由度体系振动里的程序fkmc里各数组的含义

C SHARE-BENDING MODEL THEN FIRST STOREY'S MASS IS M(1)
      SUBROUTINE FKMC(N,SB,TK,TC,H,BK,C,BM,BM0,BC0,PM,ALP0,ALP1,Q)
      REAL TK(N), BK(N,N),C(N,N),BM(N,N),BM0(N),PM(N),TC(N),H(N),
     @BC0(N),Q(N,N)
      BK=0.
      BM=0.
      C=0.
      IF(SB.GT.0.1)THEN
       DO 17 I=1, N
       DO 17 J=1, N
       AA=0.
       Q(I,J)=(H(I)*H(J)*H(1)-0.5*(H(I) H(J))*H(1)**2 H(1)**3/3.)/TC(1)
       DO 16 K=1, J-1
       AA=H(I)*H(J)*(H(K 1)-H(K))-0.5*(H(I)-H(J))*(H(K 1)**2-H(K)**2)
     @(H(K 1)**3-H(K)**3)/3.
16    Q(I,J)=Q(I,J) AA/TC(K 1)
17    Q(J,I)= Q(I,J)
      IF(SB.LT.0.9)THEN
       DO 19 I=1, N
       DO 19 J=1, I
         AA=0.
         DO 18 K=1,J
18        AA=AA 1./TK(K)
         Q(I,J)=Q(I,J) AA
19     Q(J,I)= Q(I,J)
      ENDIF
      DO 20 I=1,N
20    BK(I, I)=1.
       CALL GSN(Q,BK,N,N,1.E-11,KEY)
      ELSE
        CALL FKC(N,BK,TK)
        DO 21 I=1,N
21      WRITE(*,*)(BK(I,J),J=1,N)
      ENDIF
      DO 30 I=1,N
      BM(I,I)=BM0(I)
      PM(I)=BM0(I)
30    CONTINUE
      IF(ALP0.LT.1.E-3) THEN
       CALL FKC(N,C,BC0)
      ELSE
       DO 40 I=1, N
       DO 40 J=1, N
40     C(I,J)=ALP0*BM(I,J) ALP1*BK(I,J)
      ENDIF
      RETURN
      END