fork download
copy 
  1.          PROGRAM TRANSPOSE
  2. PARAMETER(M=3,N=6)
  3. DIMENSION A(M,N),B(M,M)
  4.  
  5. DATA B/1,2,4,2,1,-2,-3,0,5/
  6. PRINT*," MATRIX A"
  7. PRINT 20,((B(I,J),J=1,M),I=1,M)
  8. 20  FORMAT(3X,3F8.2/)
  9. PRINT*," TRANSPOSE OF THE MATRIX A"
  10. PRINT 20,((B(I,J),I=1,M),J=1,M)
  11. DO I=1,M
  12. DO J=1,N
  13. IF(J.LE.M) THEN
  14. !A(I,J)=AI(I,J)
  15. A(I,J)=B(I,J)
  16. ELSE
  17. IF(J.EQ.I+M) THEN
  18. A(I,J)=1.0
  19. ELSE
  20. A(I,J)=0.0
  21. ENDIF
  22. ENDIF
  23. ENDDO
  24. ENDDO
  25. PRINT*," AUGMENTED MATRIX"
  26. PRINT 10,((A(I,J),J=1,N),I=1,M)
  27. 10  FORMAT(3X,6F8.2/)
  28. DO I=1,M
  29. IF(A(I,I).NE.0.0) THEN
  30. P=A(I,I)
  31. ELSE
  32. PRINT*," PIVOT ELEMENT IS ZERO"
  33. STOP
  34. ENDIF
  35. DO JJ=1,N
  36. A(I,JJ)=A(I,JJ)/P
  37. ENDDO
  38. DO 3 J=1,M
  39. IF(J.EQ.I) GOTO 3
  40. COM=A(J,I)
  41. DO K=1,N
  42. A(J,K)=A(J,K)-A(I,K)*COM
  43. ENDDO
  44. 3  CONTINUE
  45. ENDDO
  46. PRINT*," AUGMENTED ECHELON MATRIX"
  47. PRINT 11,((A(I,J),J=1,N),I=1,M)
  48. 11  FORMAT(3X,6F8.2/)
  49. DO I=1,M
  50. DO J=M+1,N
  51. K=J-M
  52. B(I,K)=A(I,J)
  53. ENDDO
  54. ENDDO
  55. PRINT*," INVERSE OF MATRIX A"
  56. PRINT 30,((B(I,J),J=1,M),I=1,M)
  57. 30  FORMAT(3X,3F8.2/)
  58. END PROGRAM
Success #stdin #stdout 0s 5280KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
  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:
11 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!