C		Iwannou Spyros, Kontopodis Dhmhtrhs

		PROGRAM EROTISI 2
		DOUBLE PRECISION H(3,3),WORK(3),COND
		INTEGER I,J,N,NDIM,IPVT(3)
		N=3
		NDIM=3
		DO 10 J=1,N
		DO 10 I=1,N
		H(I,J)=1.0/FLOAT(I+J-1)
10		CONTINUE
11		FORMAT(3F8.5)

		DO 15 I=1,N
		WRITE(6,11) (H(I,J),J=1,N)
15		CONTINUE
		CALL DECOMP(NDIM,N,H,COND,IPVT,WORK)
20		FORMAT('COND=',D20.15)
		WRITE(6,20) COND
		STOP
		END