main
| Line | |
|---|
| 1 | PROGRAM MXMDRIVER
|
|---|
| 2 |
|
|---|
| 3 | INTEGER N
|
|---|
| 4 | PARAMETER (N=4)
|
|---|
| 5 | REAL A(N,N), B(N,N), C(N,N)
|
|---|
| 6 | INTEGER I,J
|
|---|
| 7 |
|
|---|
| 8 | DO J=1,N
|
|---|
| 9 | DO I=1,N
|
|---|
| 10 | A(I,J) = 1.0*(I+J)
|
|---|
| 11 | B(I,J) = 1.0*(I-J)
|
|---|
| 12 | ENDDO
|
|---|
| 13 | ENDDO
|
|---|
| 14 |
|
|---|
| 15 | C DO I=1,N
|
|---|
| 16 | C DO J=1,N
|
|---|
| 17 | C C(I,J) = 0.0
|
|---|
| 18 | C DO K=1,N
|
|---|
| 19 | C C(I,J) = C(I,J) + A(I,K)*B(K,J)
|
|---|
| 20 | C ENDDO
|
|---|
| 21 | C ENDDO
|
|---|
| 22 | C ENDDO
|
|---|
| 23 |
|
|---|
| 24 | CALL MXM(A,N,B,N,C,N)
|
|---|
| 25 |
|
|---|
| 26 | DO J=1,N
|
|---|
| 27 | DO I=1,N
|
|---|
| 28 | PRINT *,"C(",I,",",J,") = ",C(I,J)
|
|---|
| 29 | ENDDO
|
|---|
| 30 | ENDDO
|
|---|
| 31 |
|
|---|
| 32 | END
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.
