求一般实矩征特征值与特征向量

发信人: pink (hard days), 信区: Fortran
标  题: Re: [转载] 请问谁有求一般实矩征特征值与特征向量fo
发信站: 南京大学小百合站 (Tue Aug 10 21:48:42 2004)

lapack
file dgeev.f  dgeev.f plus dependencies
prec double
for  Computes the eigenvalues and left and right eigenvectors of 
     a general matrix.

【 在 leavy (☆星河☆) 的大作中提到: 】
:       subroutine jacb2(ndd,mm,Ematr,Evect,Evalu)
:       implicit double precision(a-h,o-z)
:       double precision Ematr(ndd,ndd),Evect(ndd,ndd),Evalu(ndd)
:       integer p,q
:       eps2=1.0d-12
:       Ncyc=0
:       Maxcyc=1000000
:       do 400 i=1,mm
:          Evalu(i)=0.0d0
:          Evect(i,i)=1.0d0
:          do 410 j=1,mm
:             if (i.ne.j) Evect(i,j)=0.0d0
: 410      continue
: 400   continue
: 510   fm=0d0
:       do 600 i=2,mm
:          do 610 j=1,i-1
:             if (dabs(Ematr(i,j)).gt.fm) then
:                fm=dabs(Ematr(i,j))
:                p=i
: (以下引言省略 ... ...)

--
人生最不幸处，是偶一失言而祸不及；偶一失谋而事幸成；偶一恣行而小利；
后乃视为故常，而恬不为意。则莫大之患由此生矣。

※ 来源:．南京大学小百合站 bbs.nju.edu.cn．[FROM: 172.16.66.140]