inversión de matriz de lenguaje fortran

! aa es la matriz original, b es la matriz inversa que almacena aa, n es la dimensión de la matriz aa

subrutina nizhen(aa, b, n)

entero n, i, j, k

real:: aa(n, n), b(n, n), a(n, n)

a=aa

hacer i=1, n

b(i, i)=1

enddo

hacer i=1, n< / p>

b(i,:)=b(i,:)/a(i,i)

a(i,i:n)=a(i,i:n) / a(i, i)

hacer j=i 1, n

hacer k=1, n

b(j, k)=b( j ,k)-b(i,k)*a(j,i)

enddo

a(j,i:n)=a(j,i:n) - a(i,i:n)*a(j,i)

enddo

do i=n, 1, -1

hacer j=i-1, 1, -1

hacer k=1, n

b(j, k)=b(j, k )-b (i,k)*a(j,i)

enddo

end