求矩阵行列式Fortran程序program ex01
implicit none
integer i,j,k,p,q
integer::nn=1
integer,parameter::N=3
real,dimension(N,N)::matrix
real::det=1.0
real::temp
real::E
data((matrix(i,j),i=1,3),j=1,3)/1,2,3,5,2,7,4,3,1/
do j=1,N-1
if(matrix(j,j)==0)then
do p=j+1,N
if(matr...
program ex01
implicit none
integer i,j,k,p,q
integer::nn=1
integer,parameter::N=3
real,dimension(N,N)::matrix
real::det=1.0
real::temp
real::E
data((matrix(i,j),i=1,3),j=1,3)/1,2,3,5,2,7,4,3,1/
do j=1,N-1
if(matrix(j,j)==0)then
do p=j+1,N
if(matrix(p,j)/=0)then
do q=1,N
temp=matrix(p,q)
matrix(p,q)=matrix(j,q)
matrix(j,q)=temp
end do
nn=nn*(-1)
end if
end do
end if
do i=j+1,N
E=matrix(j,j)
do k=j+1,N
matrix(i,k)=matrix(i,k)-matrix(i,j)*matrix(j,k)/E end do
matrix(i,j)=0
end do
end do
do i=1,N
det=det*matrix(i,i)
end do
det=det*nn
write(*,*) det
do i=1,N
do j=1,N
write(*,*)matrix(i,j)
end do
end do
pause
end program
本文档为【求矩阵行列式Fortran程序】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。