File: SOLVEN.FT of Tape: Various/ETH/eth11-2

(Source file text) 

---

C	SOLVEN.FOR - SAMPLE PROGRAM USING MATRICES TO TEST:
C  MINV		- MATRIX INVERSION
C  GMPRD	- GENERAL MATRIX (MODE 0) MULTIPLICATION
C  SIMQ		- SOLUTION OF SIMULTANEOUS LINEAR EQUATIONS
C
	DIMENSION A(3,3),X1(3),X2(3),B(3),C(3,3),D(10,10)
	DATA A/2.,1.,3.,3.,5.,-1.,1.,2.,1./,B/0.,6.,2./
C  OUTPUT CHANNEL = IOUT
	IOUT = 2
C
	DO 1 I=1,3
	DO 1 J=1,3
1	C(I,J)=A(I,J)
C	COPY THE ARRAY A INTO C
C
C	MINV
	CALL MINV(A,3,DETA,X1,X2)
	WRITE(IOUT,100) ((C(I,J),J=1,3),B(I),I=1,3)
100	FORMAT(/' USE OF MINV'//' (',3F8.3,' ) ( X )   (',F8.3,')'/
	1' (',3F8.3,' ) ( Y ) = (',F8.3,')'/' (',3F8.3,
	2' ) ( Z )   (',F8.3,')'//)
	WRITE(IOUT,101) ((A(I,J),J=1,3),I=1,3)
101	FORMAT(' INVERSE OF COEFFICIENT MATRIX'//
	13(' (',3F8.3,' )'/))
C
C	GMPRD
	CALL GMPRD(A,B,X1,3,3,1)
	WRITE(IOUT,102) X1
102	FORMAT(/' MULTIPLYING INVERSE BY CONSTANT VECTOR:'/
	1' X =',F8.3,/' Y =',F8.3,/' Z =',F8.3,//)
C
C	SIMQ
	CALL SIMQ(C,B,3,0)
	WRITE(IOUT,103) B
103	FORMAT(/' USE OF SIMQ FOR SOLUTION:'/
	1' X ='F8.3/' Y ='F8.3/' Z ='F8.3//)
	STOP
	END

---