[問題] 請問矩陣對角化

[sinlin (新生活!~)]

我要做一個3x3的矩陣對角化
利用"Numerical Recipes in fortran"第二版
第460頁的jacobi這個subroutine下去代
照理說應該可以求出eigenvalue和eigenvactor
可是我現在算出來的結果都只有eigenvalue是對的,eigenvactor是錯的
我直接寫一個小程式就只有run這個subroutine不做多餘的事,答案還是錯
從書上抄上去的地方我也重覆檢查了n次
請問有可能會是書上提供的程式有問題嗎??這本應該是聖經吧 囧

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.109.103.226
推 YCTzeng:矩陣對角化這個我很熟，會不會你給的矩陣不能對角化?
→ YCTzeng:如果用Numerical Recipes的話還可以試試TRED2加上TQLI
→ YCTzeng:或是你輸入的參數不對呢?導致JACOBI只計算本徵值?
→ sinlin:我給的是實數對稱矩陣,這是jacobi要求的條件,參數就我的理
→ sinlin:解應該是沒有設不對,再不行也只好用TRED2加TQLI去對照了
推 YCTzeng:可以給我那個3x3的矩陣嗎?我也來試試看。