Question

求 算一个方阵的特征值的Fortran程序代码。

rt，给我个例子就行了，网上说用lapack函数库。
我的东西很简单，所以不想细抠lapack的手册，主要就是给我看一段能直接运行的Fortran代码，简要说说参数代表的东西，就行了。 比如给你一个4阶方阵，把他的4个特征值求出来就够了。
比如下面这个三阶的矩阵A：用Fortran编，代码怎么写，然后输出三个特征值是1,2,4.？

谢谢。
矩阵A是一般形式，不一定是对称的。





Answer 1

http://fcode.cn/resource_ebook-10-1.html
这本书上有几种求特征值的算法和代码。

追问

我按照你给的这个的第三章后面的部分编的，可以转化成上海森堡矩阵，但是求特征值时就出现了这个

意思好像是 数组A（i，j）的i j 等于0了？

怎么回事，跟书上对了好几遍了，没有错误。

追答

检查88行15行是否有循环导致数组越界

追问

我也知道是这行超了，但是也没法改啊，按原文这么写确实会 使得 L值低于1的。

我也觉得奇怪，下面截图中的40下面的一行就是第88行，按照每次执行这个代码的话，矩阵元 A(L,L-1)中的L-1 就会减到0，那么就超出L最小为1的这个下边界。但是单纯改L.GT.1为L.GT.2也是不行的。这就是原文的截图。