Question

求助：利用fortran写一个矩阵求逆的程序

Answer 1

c	求逆矩阵
c	在主程序中，设两个数组a（5，5），b(5,5)
c	a(5,5)－－存放5阶方阵
c	b(5,5)－－存放单位阵
c	在子程序中，设一个数组c(5,10)，该数组是a、b阵拼接起来的
c 要实现数组的拼接和拆分，用公用语句实现(正是由于二维数组按列存储才行）

	real a(5,5),b(5,5)
	common /x/a,b
	do 10 i=1,5
	do 10 j=1,5
	if (i.eq.j) then
		b(i,j)=1
	else
		b(i,j)=0
	endif
10	continue
	call inverse
	write(*,20)b
	write(*,*)
	write(*,20)((b(i,j),j=1,5),i=1,5)
20	format(1x,5f10.3)
c	read(*,*) 
	end

	subroutine inverse
	real c(5,10)
	common /x/c
	do 10 k=1,5
		do 20 j=10,k,-1
20		c(k,j)=c(k,j)/c(k,k)
		do 40 i=1,5
		if (i.ne.k) then
			do 30 j=10,k,-1
30			c(i,j)=c(i,j)-c(i,k)*c(k,j)
		endif
40		continue
10	continue
	end

	block data
	real a(5,5),b(5,5)
	common /x/a,b
	data 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/
	end

Answer 2

c	求逆矩阵
c	在主程序中，设两个数组a（5，5），b(5,5)
c	a(5,5)－－存放5阶方阵
c	b(5,5)－－存放单位阵
c	在子程序中，设一个数组c(5,10)，该数组是a、b阵拼接起来的
c	要实现数组的拼接和拆分，用公用语句实现
	real a(5,5),b(5,5)
	common /x/a,b
	do 10 i=1,5
	do 10 j=1,5
	if (i.eq.j) then
		b(i,j)=1
	else
		b(i,j)=0
	endif
10	continue
	call inverse
	write(*,20)b
	write(*,*)
	write(*,20)((b(i,j),j=1,5),i=1,5)
20	format(1x,5f10.3)
c	read(*,*) 
	end
c     ************************************
	subroutine inverse
	real c(5,10)
	common /x/c
	do 10 k=1,5
		do 20 j=10,k,-1
20		c(k,j)=c(k,j)/c(k,k)
		do 40 i=1,5
		if (i.ne.k) then
			do 30 j=10,k,-1
30			c(i,j)=c(i,j)-c(i,k)*c(k,j)
		endif
40		continue
10	continue
	end
c ******************************
	block data
	real a(5,5),b(5,5)
	common /x/a,b
	data 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/
	end