1.fortran����Դ��
2.用fortran软件编制程序,矩矩阵计算矩阵的阵源行列式。要求使用子程序进行模块化编程
3.在fortran中矩阵C是矩矩阵(3,6),矩阵B1,阵源外卖系统的源码B2,矩矩阵B3都是阵源教育网站源码招生系统(3,2),我想让C的矩矩阵第1,2列取B1,3,阵源4列取B2,5,矩矩阵6列取B3
fortran����Դ��
参考代码:
program test
implicit none
integer :: arr(5,阵源3), arr_sum(5),itmp(3)
integer :: i,j,maxv_line(1)
print*, '输入一个5行3列整数矩阵'
do i = 1,5
read(*,*) arr(i,1:3)
arr_sum(i) = sum(arr(i,1:3))
end do
print*, '交换前的矩阵'
do i = 1,5
print*, arr(i,1:3), 'sum = ',sum(arr(i,:))
end do
maxv_line = maxloc(arr_sum)
itmp(:) = arr(maxv_line(1),:)
arr(maxv_line(1),:) = arr(1,:)
arr(1,:) = itmp(:)
print*, '交换后的矩阵'
do i = 1,5
print*, arr(i,1:3), 'sum = ',sum(arr(i,:))
end do
stop
end program
运行结果
输入
输入一个5行3列整数矩阵
1 2 3
2 3 0
3 8 1
0 1 3
5 1 2
输出
交换前的矩阵
1 2 3 sum = 6
2 3 0 sum = 5
3 8 1 sum = #和6交换
0 1 3 sum = 4
5 1 2 sum = 8
交换后的矩阵
3 8 1 sum =
2 3 0 sum = 5
1 2 3 sum = 6
0 1 3 sum = 4
5 1 2 sum = 8
用fortran软件编制程序,计算矩阵的矩矩阵行列式。要求使用子程序进行模块化编程
求矩阵行列式是阵源一个复杂的过程。有很多很多算法来做,矩矩阵qq浏览器查看源码但是阵源各有适用性。有的矩矩阵不适合病态矩阵等等。
以下是会议室预订系统 源码一个简单的全选主元高斯消去法。
摘自徐世良的《Fortran常用算法集》
Program Main
Implicit None
Real(8) :: rm(3,3) = reshape( (/1,2,4,5,7,3,,5,7/) , (/3,3/) )
Real(8) :: rDet
call BSDet( rm , 3 , rDet )
write(*,*) rDet
End Program Main
SUBROUTINE BSDET(A,N,DET)
DIMENSION A(N,N)
DOUBLE PRECISION A,DET,F,D,Q
F=1.0
DET=1.0
DO K=1,N-1
Q=0.0
DO I=K,N
DO J=K,N
IF (ABS(A(I,J)).GT.Q) THEN
Q=ABS(A(I,J))
IS=I
JS=J
END IF
CONTINUE
IF (Q+1.0.EQ.1.0) THEN
DET=0.0
RETURN
END IF
IF (IS.NE.K) THEN
F=-F
DO J=K,N
D=A(K,J)
A(K,J)=A(IS,J)
A(IS,J)=D
CONTINUE
END IF
IF (JS.NE.K) THEN
F=-F
DO I=K,N
D=A(I,JS)
A(I,JS)=A(I,K)
A(I,K)=D
CONTINUE
END IF
DET=DET*A(K,K)
DO I=K+1,N
D=A(I,K)/A(K,K)
DO J=K+1,N
A(I,J)=A(I,J)-D*A(K,J)
CONTINUE
CONTINUE
DET=F*DET*A(N,N)
RETURN
END
在fortran中矩阵C是(3,6),矩阵B1,B2,淘宝客网站源码知乎B3都是(3,2),我想让C的第1,2列取B1,3,4列取B2,5,6列取B3
① 代码{F}:program array_assignmentinteger, dimension (3, 6) :: c
integer, dimension (3, 2) :: b1, b2, b3
data b1/1,2,3,4,5,6/, b2/,,,,,/, b3/,,,,,/
c(:,1:2) = b1
c(:,3:4) = b2
c(:,5:6) = b3
do, i=1,3
write(*,'(6i5)') (c(i,j), j = 1, 6)
enddo
end program array_assignment
② 运行:
1 42 5
3 6
2024-11-30 11:41
2024-11-30 10:32
2024-11-30 10:17
2024-11-30 09:43
2024-11-30 09:42
2024-11-30 09:09