| 1 | C |
| 2 | C |
| 3 | C |
| 4 | COMMON /DDD/Z |
| 5 | |
| 6 | REAL A,B,C,Z |
| 7 | DIMENSION A(4),B(4),C(4) |
| 8 | INTEGER J |
| 9 | DIMENSION J(2) |
| 10 | |
| 11 | C A(1)=10.0 |
| 12 | C A(2)=1.22 |
| 13 | C A(3)=23.0 |
| 14 | C A(4)=1.5 |
| 15 | |
| 16 | Z=1.7070 |
| 17 | Z=1.557 |
| 18 | CALL M$INIT(A) |
| 19 | CALL M$INIT(B) |
| 20 | CALL M$INIT(C) |
| 21 | 10 FORMAT (9HMATRIX A:) |
| 22 | 11 FORMAT (9HMATRIX B:) |
| 23 | 12 FORMAT (9HMATRIX C:) |
| 24 | |
| 25 | |
| 26 | C A(1)=1.0 |
| 27 | C A(2)=2.0 |
| 28 | C A(3)=3.0 |
| 29 | C A(4)=4.0 |
| 30 | |
| 31 | WRITE (1,10) |
| 32 | CALL PMAT(A) |
| 33 | CALL M$ROT(A,Z) |
| 34 | CALL PMAT (A) |
| 35 | PAUSE 4 |
| 36 | C A(1)=1.0 |
| 37 | C A(2)=2.0 |
| 38 | C A(3)=3.0 |
| 39 | C A(4)=4.0 |
| 40 | |
| 41 | B(1)=1.7 |
| 42 | B(4)=1.7 |
| 43 | C B(2)=6.0 |
| 44 | C B(3)=7.0 |
| 45 | C B(4)=-8.0 |
| 46 | |
| 47 | CALL M$MUL(C,A,B) |
| 48 | |
| 49 | WRITE (1,10) |
| 50 | CALL PMAT(A) |
| 51 | WRITE (1,11) |
| 52 | CALL PMAT(B) |
| 53 | WRITE (1,12) |
| 54 | CALL PMAT(C) |
| 55 | |
| 56 | J(1)=4 |
| 57 | J(2)=10 |
| 58 | |
| 59 | 13 FORMAT (8HPOINT J:) |
| 60 | WRITE (1,13) |
| 61 | CALL PPOINT(J) |
| 62 | |
| 63 | CALL M$APLI(C,J) |
| 64 | |
| 65 | 14 FORMAT (30HPOINT J AFTER TRANSFORMATION: ) |
| 66 | WRITE (1,14) |
| 67 | CALL PPOINT(J) |
| 68 | |
| 69 | |
| 70 | 22 FORMAT (20HPROGRAMMENDE ) |
| 71 | |
| 72 | |
| 73 | WRITE (1,22) |
| 74 | |
| 75 | PAUSE 2 |
| 76 | 99 GO TO 99 |
| 77 | END |
| 78 | |
| 79 | CCCCCCCCCCCCCCCCCCCCCCC |
| 80 | |
| 81 | SUBROUTINE PMAT (MATRIX) |
| 82 | |
| 83 | DIMENSION MATRIX(4) |
| 84 | REAL MATRIX |
| 85 | REAL M1,M2,M3,M4 |
| 86 | M1=MATRIX(0) |
| 87 | M2=MATRIX(1) |
| 88 | M3=MATRIX(2) |
| 89 | M4=MATRIX(3) |
| 90 | |
| 91 | WRITE (1,343) MATRIX |
| 92 | WRITE (1,344) |
| 93 | |
| 94 | 343 FORMAT(2H[ ,F18.10,2X,F18.10,2H ]) |
| 95 | 344 FORMAT(X) |
| 96 | |
| 97 | RETURN |
| 98 | END |
| 99 | |
| 100 | CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
| 101 | |
| 102 | SUBROUTINE PPOINT(POINT) |
| 103 | DIMENSION POINT(2) |
| 104 | INTEGER POINT |
| 105 | |
| 106 | WRITE (1,200) POINT |
| 107 | 200 FORMAT(2H( ,I4,2H, I4,2H )) |
| 108 | RETURN |
| 109 | END |
| 110 | |
| 111 | |
| 112 | CCCCCCCCCCCCCCCCCCC |
| 113 | $0 |