C Iwannou Spyros, Kontopodhs Dhmhtrhs PROGRAM EROTISI 3 REAL A(3,3),R(3,3) INTEGER I,N,NDIM NDIM=3 N=3 DATA ((A(I,J),I=1,3),J=1,3)/5.,-1.5,2.5,-3.5,1.,-0.5,0.,3.,2.5/ 1 FORMAT(3F10.5) WRITE(6,*) 'Pinakas A' DO 2 I=1,N WRITE(6,1) (A(I,J),J=1,N) 2 CONTINUE WRITE(6,*) 1H CALL INVERT(A,R,N,NDIM) WRITE(6,*) 'Antistrofos tou A' DO 4 I=1,N WRITE(6,1) (R(I,J),J=1,N) 4 CONTINUE CALL INVERT(R,A,N,NDIM) WRITE(6,*) 'Antistrofos tou antistrofou tou A' DO 5 I=1,N WRITE(6,1) (A(I,J),J=1,N) 5 CONTINUE STOP END SUBROUTINE INVERT(A,R,N,NDIM) INTEGER N,IPVT(N),I,J,NDIM REAL R(NDIM,N),A(NDIM,N),WORK(N),COND,T(N,N),M(N) DO 10 J=1,N DO 10 I=1,N T(I,J)=A(I,J) 10 CONTINUE CALL DECOMP(NDIM,N,A,COND,IPVT,WORK) IF ((COND+1).EQ.COND) THEN WRITE(6,*) 'PINAKAS MH ANTISTREPSIMOS,PROBLHMA STHN KLHSH INVERT' RETURN ENDIF DO 60 J=1,N DO 50 I=1,N IF (I .EQ. J) THEN M(I)=1 ELSE M(I)=0 ENDIF 50 CONTINUE CALL SOLVE(NDIM,N,A,M,IPVT) DO 55 I=1,N R(I,J)=M(I) 55 CONTINUE 60 CONTINUE DO 70 J=1,N DO 70 I=1,N C Epanaferei ton A pou prohgoumenws eixe krathsei ston pinaka T A(I,J)=T(I,J) 70 CONTINUE RETURN END