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> 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