请问如何用MATLAB解"久期方程"

我知道的Matlab解矩阵方法有：
A=[1 6 9;6 5 2;8 3 4]
b=[2;4;7]
1）x=inv(A)*b 【就是你用的方法】
2）x=A\b
3）x=A^(-1)*b
4）C=[A,b]
rref(C)
试试吧，也许好使。
还有一个网址，下载一些经典程序，如：
gmres_m.rar-gmres算法的一个特殊形式，用来求解大型稀疏矩阵方程,matlab
GMres.rar-数值算法中另一个经典算法gmres算法，用来求解大型矩阵方程问题。,matlab
arnoldi.rar-数值分析中经典的arnoldi算法，用来求解大型矩阵方程组的求解,matlab
http://www.programsalon.com/sitemap/sitemap153_1450.htm

补充：
2．利用矩阵的LU、QR和cholesky分解求方程组的解
（1）LU分解：
LU分解又称Gauss消去分解，可把任意方阵分解为下三角矩阵的基本变换形式（行交换）和上三角矩阵的乘积。即A=LU，L为下三角阵，U为上三角阵。
则：A*X=b 变成L*U*X=b
所以X=U\(L\b) 这样可以大大提高运算速度。
命令 [L，U]=lu (A)
例1-78 求方程组 的一个特解。
解：

>>A=[4 2 -1;3 -1 2;11 3 0];
>>B=[2 10 8]';
>>D=det(A)
>>[L,U]=lu(A)
>>X=U\(L\B)
显示结果如下：
D =
0