2个回答
展开全部
以下是求5阶方阵的逆矩阵,fortran77 程序。供参考。
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
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询