※ 引述《roussas (QQ)》之銘言： : : http://www.na.ntu.edu.tw/chinese/course/online/num_anal/chap2.htm : : 這裡有演算法 不知道有沒有幫助 : : 不然~~去圖書館借 : : 數值方法入門-使用C語言- 全華科技圖書 : : 一共有兩本 : : 裡面有C寫的程式跟講解 Jacobi Method的作法 : : 星期三我借的時候還有一本 一本我借走了 : : 希望對你有幫助 : http://www.cyut.edu.tw/~ckhung/b/la/gje.shtml : 這裡也有高斯jordan的做法... 如果是用 Gauss-Seidel 法 必須要考慮到數列收斂的問題 因為用這方法的前提是可以找到一初始值使得n至無窮時x會收斂到真正的值 所以就算方程組有解 用此法可能會成數列發散而找不到解 如果掉換一下各方程式的位置 重新寫一組矩陣 或許就可收斂了 或是改變初始值 ex: x+2y+4z=1 3y+z=1 2x+y+3z=1 (x,y,z)=(1/3,1/3,0) 用Gauss-Seidel法用不出來 但將第一列第三列交換數列就收斂了........ --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.113.185.114