C
C
C     
      COMMON /DDD/Z

      REAL A,B,C,Z
      DIMENSION A(4),B(4),C(4)
      INTEGER J
      DIMENSION J(2)
      
C      A(1)=10.0
C      A(2)=1.22
C      A(3)=23.0
C      A(4)=1.5

      Z=1.7070
      Z=1.557
      CALL M$INIT(A)
      CALL M$INIT(B)
      CALL M$INIT(C)
 10   FORMAT (9HMATRIX A:)
 11   FORMAT (9HMATRIX B:)
 12   FORMAT (9HMATRIX C:)
      

C      A(1)=1.0
C      A(2)=2.0
C      A(3)=3.0
C      A(4)=4.0

      WRITE (1,10)
      CALL PMAT(A)
      CALL M$ROT(A,Z)
      CALL PMAT (A)
      PAUSE 4
C      A(1)=1.0
C      A(2)=2.0
C      A(3)=3.0
C      A(4)=4.0
      
      B(1)=1.7
      B(4)=1.7
C      B(2)=6.0
C      B(3)=7.0
C      B(4)=-8.0
      
      CALL M$MUL(C,A,B)
      
      WRITE (1,10)
      CALL PMAT(A)
      WRITE (1,11)
      CALL PMAT(B)
      WRITE (1,12)
      CALL PMAT(C)

      J(1)=4
      J(2)=10
      
 13   FORMAT (8HPOINT J:)
      WRITE (1,13)
      CALL PPOINT(J)

      CALL M$APLI(C,J)

 14   FORMAT (30HPOINT J AFTER TRANSFORMATION:   )
      WRITE (1,14)
      CALL PPOINT(J)


 22   FORMAT (20HPROGRAMMENDE            )

     
      WRITE (1,22)
 
      PAUSE 2
 99   GO TO 99
      END

CCCCCCCCCCCCCCCCCCCCCCC

      SUBROUTINE PMAT (MATRIX)

      DIMENSION MATRIX(4)
      REAL MATRIX
      REAL M1,M2,M3,M4
      M1=MATRIX(0)
      M2=MATRIX(1)
      M3=MATRIX(2)
      M4=MATRIX(3)

      WRITE (1,343) MATRIX
      WRITE (1,344)

 343  FORMAT(2H[ ,F18.10,2X,F18.10,2H ])
 344  FORMAT(X)

      RETURN
      END

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

      SUBROUTINE PPOINT(POINT)
      DIMENSION POINT(2)
      INTEGER POINT

      WRITE (1,200) POINT
 200  FORMAT(2H( ,I4,2H, I4,2H ))
      RETURN
      END


CCCCCCCCCCCCCCCCCCC
$0