急に

連立方程式を解かないといけなくなって、
こういうとき、どうしようかと考えた。

ここで、やっぱ、Rだよね〜

Ax=b
なるのを解くとして、

> a <- matrix(c(8,2,0,1,2,8,0,1,0,0,8,1,-1,-1,-1,0),4,4)
> a
[,1] [,2] [,3] [,4]
[1,] 8 2 0 -1
[2,] 2 8 0 -1
[3,] 0 0 8 -1
[4,] 1 1 1 0
> b <- matrix(c(5.4,3.4,1.2,1.0))
> b
[,1]
[1,] 5.4
[2,] 3.4
[3,] 1.2
[4,] 1.0
> solve(a,b)
[,1]
[1,] 0.5974359
[2,] 0.2641026
[3,] 0.1384615
[4,] -0.0923077

で終わり。
やっぱ、早いわ。

お手軽でした＾＾