線型連立方程式を解く 以前はこちらやこちら Gauss-Seidel 法 , は の対角成分からなる対角行列 対角成分が大きくないといけないらしい N<-3 A<-matrix(runif(N^2),N,N) diag(A)<-diag(A)*10 #対角成分が大きくないといけない x<-runif(N) b<-runif(N) D<-ma…