ce1f6adb |
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 | |
8eb88117 |
11 | C A(1)=10.0 |
12 | C A(2)=1.22 |
13 | C A(3)=23.0 |
14 | C A(4)=1.5 |
ce1f6adb |
15 | |
8eb88117 |
16 | Z=1.7070 |
17 | Z=1.557 |
ce1f6adb |
18 | CALL M$INIT(A) |
19 | CALL M$INIT(B) |
20 | CALL M$INIT(C) |
ce1f6adb |
21 | 10 FORMAT (9HMATRIX A:) |
22 | 11 FORMAT (9HMATRIX B:) |
23 | 12 FORMAT (9HMATRIX C:) |
24 | |
8eb88117 |
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 | |
ce1f6adb |
31 | WRITE (1,10) |
8eb88117 |
32 | CALL PMAT(A) |
33 | CALL M$ROT(A,Z) |
ce1f6adb |
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 | |
8eb88117 |
94 | 343 FORMAT(2H[ ,F18.10,2X,F18.10,2H ]) |
ce1f6adb |
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 |