30 CONTINUE

END

*-------------------------------------------

C SUBPROGRAM-9

C SOLVE EQUATIONS BY GS METHOD SUBROUTINE BGSMT(NJ,NJ2,NW,KS,P)

DIMENSION P(NJ2)

REAL KS(NJ2,NW)

NJ1=NJ2-1

DO 20 K=1,NJ1

IF(NJ2.GT.K+NW-1) IM=K+NW-1

IF(NJ2.LE.K+NW-1) IM=NJ2

K1=K+1

DO 20 I=K1,IM

L=I-K+1

C=KS(K,L)/KS(K,1)

IW=NW-L+1

DO 10 J=1,IW

M=J+I-K

10 KS(I,J)=KS(I,J)-C*KS(K,M)

20 P(I)=P(I)-C*P(K)

P(NJ2)=P(NJ2)/KS(NJ2,1)

DO 40 I1=1,NJ1

I=NJ2-I1

IF(NW.GT.NJ2-I+1) JO=NJ2-I+1

IF(NW.LE.NJ2-I+1) JO=NW

DO 30 J=2,JO

K=J+I-1

30 P(I)=P(I)-KS(I,J)*P(K)

40 P(I)=P(I)/KS(I,1)

WRITE(6,50)

50 FORMAT(/1X,'NODAL DISPLACEMENTS'/3X, * 'NODE',5X,'X-DISP.',8X,'Y-DISP.')

DO 60 I=1,NJ

60 WRITE(6,70) I,P(2*I-1),P(2*I)

70 FORMAT(1X,I5,2E15.6)

END

*---------------------------------------------------

C SUBPROGRAM-10

C CALCULATE ELEMENT STRESS MATRIX SUBROUTINE SIGME(NE,NJ,NJ2,E,PR,LND,X,Y,P)