fortran语言矩阵求逆
Posted
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了fortran语言矩阵求逆相关的知识，希望对你有一定的参考价值。
利用fartran语言做一个套用 子程序 矩阵求逆，谁会，帮帮忙啊！！！
! aa为原矩阵，b为存放aa的逆矩阵，n为矩阵aa的维数subroutine nizhen(aa,b,n)
integer n,i,j,k
real:: aa(n,n),b(n,n),a(n,n)
a=aa
do i=1,n
b(i,i)=1
enddo
do i=1,n
b(i,:)=b(i,:)/a(i,i)
a(i,i:n)=a(i,i:n)/a(i,i)
do j=i+1,n
do 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
enddo
do i=n,1,-1
do j=i-1,1,-1
do k=1,n
b(j,k)=b(j,k)-b(i,k)*a(j,i)
enddo
enddo
enddo
end 参考技术A 没有必要自己编吧。
1 调用IMSL/MKL等现成的数学库，比如.i.就可以直接求逆。
2 直接copy徐士良那本<<Fortran常用算法集>>上的程序
以上是关于fortran语言矩阵求逆的主要内容，如果未能解决你的问题，请参考以下文章