- 快速注册天涯账号
已有天涯账号？登陆
- 这里是匿名提问所提的问题，您需要登录才能参与回答。
- "天涯问答"是天涯社区旗下的问题分享平台。在这里您可以提问，回答感兴趣的问题，分享知识和经历，无论您在何时何地上线都可以访问，此平台完全免费，而且注册非常简单。
谁能给出一个求矩阵对角化的FORTRAN源程序啊？
谁能给出一个求矩阵对角化的FORTRAN源程序啊？
10-03-19 匿名提问 发布
-
ty_133797018 赞成
Program Main
[code]use PreciseTimeIntegration
Real(8), Dimension(8, 8) :: K,M,C,INVK
Real(8), Dimension(8, 4) :: R,U,V,A
Real(8) dt
dt=1
R=0
C=0
C(7,7)=5
C(7,8)=-5
C(8,7)=-5
C(8,8)=5
DO i=1,8
M(i,i)=8
END DO
DO i=1,7
K(i,i+1)=-4
k(i+1,i)=-4
END DO
DO i=2,7
k(i,i)=8
END DO
K(1,1)=4
K(8,8)=4
! WRITE (*,*) M
! WRITE (*,*) K
! WRITE (*,*) C
! INVK=INV(K,8)
! STOP
U=0
U(8,1)=10
V=0
A=0
call PreciseTIM (M, C, k, R, 8, 4, dt, U, V, A)
write (*,*) U
! write (*,*) V
! write (*,*) A
end[/code]◆◆
ty_133797018
17-10-18 |