fortran语言矩阵求逆
来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/08/19 06:44:03
fortran语言矩阵求逆
利用fartran语言做一个套用 子程序 矩阵求逆,谁会,帮帮忙啊!
利用fartran语言做一个套用 子程序 矩阵求逆,谁会,帮帮忙啊!
![fortran语言矩阵求逆](/uploads/image/z/5901866-26-6.jpg?t=fortran%E8%AF%AD%E8%A8%80%E7%9F%A9%E9%98%B5%E6%B1%82%E9%80%86)
!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
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