Rで連立一次方程式の解を求める。

Rで連立一次方程式の解を求める。
solve関数を使える。
http://cse.naro.affrc.go.jp/takezawa/r-tips/r/20.html

a <- matrix(c(0,1,2,3,4,5,6,7,9), 3,3)     
b<-matrix(c(1,0,-2))                     
solve(a,b)  
                                

#      3y + 6z =  1
#  x + 4y + 7z =  0
# 2x + 5y + 9z = -2
>|sh|
          [,1]
[1,] -2.333333
[2,]  2.333333
[3,] -1.000000

ガウス=サイデル法
http://homepage3.nifty.com/gakuyu/suti/renritu/hanpuku/index.html

a<-matrix(c(15,2,6,2,34,5,3,5,38),3,3)
b<-matrix(c(10,12,21))
kurikaesi<-7
x<-c(1,1,1)
diag<-matrix(a)

for(k in 1:kurikaesi){

	for(i in 1:length(x)){
		x[i]<-(b[i] - sum( a[i,][-i]*x[-i] ) )/a[i,i]
	}      


}

x
solve(a,b)