fortran语言矩阵求逆利用fartran语言做一个套用 子程序 矩阵求逆,谁会,帮帮忙啊!

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/05 20:31:30
fortran语言矩阵求逆利用fartran语言做一个套用 子程序   矩阵求逆,谁会,帮帮忙啊!

fortran语言矩阵求逆利用fartran语言做一个套用 子程序 矩阵求逆,谁会,帮帮忙啊!
fortran语言矩阵求逆
利用fartran语言做一个套用 子程序 矩阵求逆,谁会,帮帮忙啊!

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