Конвертация Fortran-77 в Fortran-90

      SUBROUTINE ASG0R(A,B,X,N,P)
      implicit none
      DIMENSION A(N,N),B(N),X(N)
      INTEGER N,P,F,L,M,I,J,K
      REAL A,B,X,S,R,W
      IF(N.GT.1) GO TO 15
      IF(N.LE.0) GO TO 10
      X(N)=B(N)/A(N,N)
      GO TO 10
   15 IF(P.NE.1) GO TO 4
      F=0
    1 F=F+1
      S=1./A(F,F)
      L=F+1
      DO 2 I=L,N,1
    2 A(I,F)=-(S*A(I,F))
      DO 3 J=L,N,1
      DO 3 I=L,N,1
    3 A(I,J)=A(I,J)+A(F,J)*A(I,F)
      IF(L-N) 1,4,4
    4 DO 5 I=1,N
    5 X(I)=B(I)
      M=N-1
      DO 6 J=1,M
      W=X(J)
      K=J+1
      DO 6 I=K,N
    6 X(I)=X(I)+A(I,J)*W
      X(N)=X(N)/A(N,N)
      IF(N.EQ.1) GO TO 10
      K=N
      DO 9 I=1,M
      R=0.0
      L=K-1
      DO 8 J=K,N
    8 R=R+A(L,J)*X(J)
      K=N-I
    9 X(K)=(X(K)-R)/A(K,K)
   10 RETURN
      END

!*==ASG0R.spg  processed by SPAG 6.72Dc at 23:43 on 17 Aug 2019
      SUBROUTINE ASG0R(A,B,X,N,P)
      IMPLICIT NONE
!*--ASG0R4
      DIMENSION A(N,N) , B(N) , X(N)
      INTEGER N , P , f , l , m , i , j , k
      REAL A , B , X , s , r , w
      IF ( N>1 ) THEN
         IF ( P==1 ) THEN
            f = 0
 20         f = f + 1
            s = 1./A(f,f)
            l = f + 1
            DO i = l , N , 1
               A(i,f) = -(s*A(i,f))
            ENDDO
            DO j = l , N , 1
               DO i = l , N , 1
                  A(i,j) = A(i,j) + A(f,j)*A(i,f)
               ENDDO
            ENDDO
            IF ( l<N ) GOTO 20
         ENDIF
         DO i = 1 , N
            X(i) = B(i)
         ENDDO
         m = N - 1
         DO j = 1 , m
            w = X(j)
            k = j + 1
            DO i = k , N
               X(i) = X(i) + A(i,j)*w
            ENDDO
         ENDDO
         X(N) = X(N)/A(N,N)
         IF ( N/=1 ) THEN
            k = N
            DO i = 1 , m
               r = 0.0
               l = k - 1
               DO j = k , N
                  r = r + A(l,j)*X(j)
               ENDDO
               k = N - i
               X(k) = (X(k)-r)/A(k,k)
            ENDDO
         ENDIF
      ELSEIF ( N>0 ) THEN
         X(N) = B(N)/A(N,N)
      ENDIF
      END