IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION BM(3,8),BT(8,3),D(3,3),SON(8,8),CARP(8,3),BON(8,8)
     *,N(140),SR(141,41),SRM(140,140),YUK(140),DP(140)
      OPEN(2,FILE='PLAK.DAT',STATUS='OLD')
      OPEN(3,FILE='PLAK.OUT',FORM='FORMATTED')
      READ(2,*)NP,ND
      READ(2,*)E,V,A,B,H
      DO 20 PSAY=-1/SQRT(3),1/SQRT(3),2/SQRT(3)
      DO 20 ETA=-1/SQRT(3),1/SQRT(3),2/SQRT(3)
      BM(1,1)=-(1.-ETA)/(2.*A)
      BM(1,3)=-BM(1,1)
      BM(1,5)=(1.+ETA)/(2.*A)
      BM(1,7)=-(1.+ETA)/(2.*A)
      BM(2,2)=-(1.-PSAY)/(2.*B)
      BM(2,4)=-(1.+PSAY)/(2.*B)
      BM(2,6)=-BM(2,4)
      BM(2,8)=-BM(2,2)
      BM(3,1)=BM(2,2)
      BM(3,2)=BM(1,1)
      BM(3,3)=BM(2,4)
      BM(3,4)=BM(1,3)
      BM(3,5)=BM(2,6)
      BM(3,6)=BM(1,5)
      BM(3,7)=BM(2,8)
      BM(3,8)=BM(1,7)
      RJ=(A*B)/4.
      DO 10 I=1,8
      DO 10 J=1,3
   10 BT(I,J)=BM(J,I)
      RK=E/(1-V*V)
      RJ=(A*B)/4.
      D(1,1)=1*E/(1-V*V)
      D(1,2)=V*E/(1-V*V)
      D(2,1)=V*E/(1-V*V)
      D(2,2)=1*E/(1-V*V)
      D(3,3)=E*(1-V)/(2.*(1-V*V))
      CALL MATCARP(3,3,D,8,3,BT,CARP)
      CALL MATCARP(3,8,BM,8,3,CARP,BON)
      DO 20 I=1,8
      DO 20 J=1,8
   20 SON(I,J)=SON(I,J)+RJ*H*BON(I,J)
      CALL BIALT(E,V,A,B,ND,NP,SON,N,SR,SRM,YUK,DP)
      CLOSE(3)
      CLOSE(2)
      STOP
      END

      SUBROUTINE BIALT(E,V,A,B,ND,NP,SON,N,SR,SRM,YUK,DP)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION SON(8,8),N(ND),SR(ND+1,ND+1),SRM(ND,ND)
     *,YUK(ND),DP(ND)
C ====================================================
C ===== SISTEM RIJITLIK MATRISININ OLUSTURULMASI =====
C ====================================================
      DO 50 INP=1,NP
      READ(2,*)(N(I),I=1,8)
      DO 40 I=1,8
   40 N(I)=N(I)+1
      DO 50 I=1,8
      DO 50 J=1,8
   45 FORMAT(I2,'-)',10(1PE10.2))
   50 SR(N(I),N(J))=SR(N(I),N(J))+SON(I,J)
      DO 60 I=1,ND
      DO 60 J=1,ND
   60 SRM(I,J)=SR(I+1,J+1)
   80 FORMAT(I2,'-)',10(1PE10.2))
      WRITE(3,*)'SISTEM RIJITLIK MATRISI'
      DO 90 I=1,ND
   90 WRITE(3,80)I,(SRM(I,L),L=1,ND)
C ====================================================
C ==== SISTEM RIJITLIK MATRISININ TERSININ HESABI ====
C ====================================================
      DO 110 I=1,ND
      DO 120 J=1,ND
      DO 120 K=1,ND
      IF(J.EQ.I.OR.K.EQ.I)GO TO 120
      SRM(J,K)=SRM(J,K)-SRM(J,I)*SRM(I,K)/SRM(I,I)
  120 CONTINUE
      SRM(I,I)=-1.0/SRM(I,I)
      DO 130 J=1,ND
      IF(J.EQ.I)GO TO 130
      SRM(J,I)=SRM(J,I)*SRM(I,I)
  130 CONTINUE
      DO 140 K=1,ND
      IF(K.EQ.I)GO TO 140
      SRM(I,K)=SRM(I,K)*SRM(I,I)
  140 CONTINUE
  110 CONTINUE
      DO 150 I=1,ND
      DO 150 J=1,ND
  150 SRM(I,J)=-SRM(I,J)
      WRITE(3,*)'SISTEM RIJITLIK MATRISININ TERSI'
      DO 190 I=1,ND
  190 WRITE(3,80)I,(SRM(I,L),L=1,ND)
C ====================================================
C ==== SISTEME ETKIYEN YUKLERIN KUTUKTEN OKUNMASI ====
C ====================================================
      READ(2,*)(YUK(I),I=1,ND)
      WRITE(3,*)'SISTEME ETKIYEN YUKLER'
      DO 200 I=1,ND
  200 WRITE(3,85)I,YUK(I)
   85 FORMAT(I2,'-)',1PE10.2)
C ====================================================
C ============== DEFORMASYONLARIN HESABI =============
C ====================================================
      CALL MATCARP(ND,1,YUK,ND,ND,SRM,DP)
      WRITE(3,*)'DEFORMASYONLAR'
      DO 210 I=1,ND
  210 WRITE(3,85)I,DP(I)
C ====================================================
C ================== GERILME HESABI ==================
C ====================================================
      WRITE(3,*)'GERILMELER'
      DO 220 I=1,ND-1,2
      PS=E/(1-V*V)
      SIGMAX=PS*(DP(I)/A-V*DP(I+1)/B)
      SIGMAY=PS*(V*DP(I)/A-DP(I+1)/B)
      WRITE(3,85)I,SIGMAX
  220 WRITE(3,85)I+1,SIGMAY
      RETURN
      END
C =======================================================
C ======= IKI MATRISIN CARPIMINI YAPAN ALT PROGRAM ======
C =======================================================
      SUBROUTINE MATCARP(ID,JD,D,IB,JB,BT,CARP)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION BT(IB,JB),D(ID,JD),CARP(IB,JD)
      DO 100 I=1,IB
      DO 100 J=1,JD
      CARP(I,J)=0.
      DO 100 K=1,JB
  100 CARP(I,J)=CARP(I,J)+BT(I,K)*D(K,J)
      RETURN
      END