ガウス-ジョルダン法

• 線形の連立方程式の解法
  • こちらのつづき
• 掃き出し法に近いもので、Gauss-Jordan法
• 拡大係数行列を変形させてみる

N<-3
A<-matrix(runif(N*(N+1)),N,N+1)

C<-A
for(i in 1:N){
C[i,]<-C[i,]/C[i,i]
	if(i != N){
	for(j in (i+1):N){
	C[j,]<-C[j,]-C[j,i]*C[i,]
	}}
}

for(i in N:1){
	if(i != 1){
	for(j in 1:(i-1)){
	C[j,]<-C[j,]-C[j,i]*C[i,]
	}}
}

C