Directory of image this file is from
This file as a plain text file

	SUBROUTINE MINV
C MATRIX INVERSION - DIMENSION N<=16
C MATRIX A : A11,A12,A13,...,ANN IN COMMON
C FORTRAN II
	COMMON N,A,B
	DIMENSION A(272),B(256)
	NN=N*N
	DO 1020 J=1,N
	DO 1000 I=1,N
	A(NN+I)=0.
1000	CONTINUE
	A(NN+J)=1.
	PV=A(J)
	K=NN+J
	DO 1001 L=J,K,N
	A(L)=A(L)/PV
1001	CONTINUE
	DO 1005 M=1,N
	IF(M-J)1002,1005,1002
1002	M1=M
	M2=M1+NN
	M3=J
	RW=A(M1)
	DO 1003 K1=M1,M2,N
	A(K1)=A(K1)-RW*A(M3)
	M3=M3+N
1003	CONTINUE
1005	CONTINUE
	DO 1010 I=1,NN
	A(I)=A(I+N)
1010	CONTINUE
1020	CONTINUE
	RETURN
	END