快速注册天涯账号

已有天涯账号?登陆

这里是匿名提问所提的问题,您需要登录才能参与回答。
"天涯问答"是天涯社区旗下的问题分享平台。在这里您可以提问,回答感兴趣的问题,分享知识和经历,无论您在何时何地上线都可以访问,此平台完全免费,而且注册非常简单。

谁能给出一个求矩阵对角化的FORTRAN源程序啊?

谁能给出一个求矩阵对角化的FORTRAN源程序啊?
10-03-19  匿名提问 发布
  • 1

    王者归来wp

    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]

    10-04-11 | 1条评论 | 打赏

    ty_133797018

    看不懂,也不会用,能再说一下吗

    17-10-18 | 回复

意见反馈 帮助