1 subroutine gauss(RO,AP,MT,M,N,*)
3 c CALCULATES (RO**(-1))*AP BY GAUSS ELIMINATION
4 c RO IS A SQUARE MATRIX
5 c THE CALCULATED PRODUCT IS STORED IN AP
6 c ABNORMAL EXIT IF RO IS SINGULAR
8 integer MT, M, N, M1,I,J,IM,
10 double precision RO(MT,M),AP(MT,N),X,RM,PR,
14 if(dabs(X).le.1.0D-13) return 1
26 if(DABS(RO(J,I)).LE.RM) goto 2
40 if(dabs(X).le.1.0E-13) return 1
49 8 AP(J,K)=AP(J,K)-Y*AP(I,K)
51 9 RO(J,K)=RO(J,K)-Y*RO(I,K)
55 if(dabs(X).le.1.0E-13) return 1
65 15 X=X-AP(K,J)*RO(MI,K)