matriz inversa de programación fortran

El siguiente es un programa fortran77 para encontrar la matriz inversa de una matriz cuadrada de quinto orden. Como referencia.

c Encuentre la matriz inversa

c En el programa principal, suponga dos matrices a(5,5), b(5,5)

c a ( 5,5)——Almacenamiento de la matriz cuadrada de quinto orden

c b(5,5)——Almacenamiento de la matriz unitaria

c En la subrutina, establezca una matriz c(5, 10), esta matriz está empalmada por las matrices a y b

c Para realizar el empalme y división de matrices, utilice declaraciones comunes

real?a(5,5) ,b(5,5)

común?/x/a,b

do?10?i=1,5

do?10?j =1,5

si?(es decir,eq.j)?entonces

b(i,j)=1

más

b(i,j)=0

endif

10 continuar

llamada?inversa

escribir(*,20) b

escribir(*,*)

escribir(*,20)((b(i,j),j=1,5),i=1,5)

20 formato(1x,5f10.3)

c lectura(*,*)?

fin

c?*** *********************************

¿subrutina?inversa

real?c(5,10)

común?/x/c

do?10?k=1,5

do?20?j= 10,k,-1

20 c(k,j)=c(k,j)/c(k,k)

hacer?40?i=1,5

si?(i.ne.k)?entonces

hacer?30?j=10,k,-1

30 c(i,j )=c(i,j)-c(i,k)*c(k,j)

endif

40 continuar

10 continuar p>

fin

c?******************************

bloque?datos

real?a(5,5),b(5,5)

común?/x/a,b

datos ?a/3,1,0,2,10,-2,0,1,3,1,9,3,1,0,1, 1 1,0,1,1,0,1,2 ,0 ,2,10/

fin