fork download 
         PROGRAM TRANSPOSE
PARAMETER(M=3,N=6)
DIMENSION A(M,N),B(M,M)

DATA B/1,2,4,2,1,-2,-3,0,5/
PRINT*," MATRIX A"
PRINT 20,((B(I,J),J=1,M),I=1,M)
20  FORMAT(3X,3F8.2/)
PRINT*," TRANSPOSE OF THE MATRIX A"
PRINT 20,((B(I,J),I=1,M),J=1,M)
DO I=1,M
DO J=1,N
IF(J.LE.M) THEN
!A(I,J)=AI(I,J)
A(I,J)=B(I,J)
ELSE
IF(J.EQ.I+M) THEN
A(I,J)=1.0
ELSE
A(I,J)=0.0
ENDIF
ENDIF
ENDDO
ENDDO
PRINT*," AUGMENTED MATRIX"
PRINT 10,((A(I,J),J=1,N),I=1,M)
10  FORMAT(3X,6F8.2/)
DO I=1,M
IF(A(I,I).NE.0.0) THEN
P=A(I,I)
ELSE
PRINT*," PIVOT ELEMENT IS ZERO"
STOP
ENDIF
DO JJ=1,N
A(I,JJ)=A(I,JJ)/P
ENDDO
DO 3 J=1,M
IF(J.EQ.I) GOTO 3
COM=A(J,I)
DO K=1,N
A(J,K)=A(J,K)-A(I,K)*COM
ENDDO
3  CONTINUE
ENDDO
PRINT*," AUGMENTED ECHELON MATRIX"
PRINT 11,((A(I,J),J=1,N),I=1,M)
11  FORMAT(3X,6F8.2/)
DO I=1,M
DO J=M+1,N
K=J-M
B(I,K)=A(I,J)
ENDDO
ENDDO
PRINT*," INVERSE OF MATRIX A"
PRINT 30,((B(I,J),J=1,M),I=1,M)
30  FORMAT(3X,3F8.2/)
END PROGRAM
Success #stdin #stdout 0s 5280KB
comments (?)
stdin 
Standard input is empty
stdout 
  MATRIX A
       1.00    2.00   -3.00

       2.00    1.00    0.00

       4.00   -2.00    5.00

  TRANSPOSE OF THE MATRIX A
       1.00    2.00    4.00

       2.00    1.00   -2.00

      -3.00    0.00    5.00

  AUGMENTED MATRIX
       1.00    2.00   -3.00    1.00    0.00    0.00

       2.00    1.00    0.00    0.00    1.00    0.00

       4.00   -2.00    5.00    0.00    0.00    1.00

  AUGMENTED ECHELON MATRIX
       1.00    0.00    0.00    0.56   -0.44    0.33

      -0.00    1.00    0.00   -1.11    1.89   -0.67

      -0.00   -0.00    1.00   -0.89    1.11   -0.33

  INVERSE OF MATRIX A
       0.56   -0.44    0.33

      -1.11    1.89   -0.67

      -0.89    1.11   -0.33
https://ideone.com/A8dxGW
language:
Fortran (gfortran 8.3)
created:
23 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!